Merge commit 'bundle/olesenm' into home

ReleaseNotes-1.5

@@ -2,7 +2,7 @@
 #
 #+TITLE:          *OpenFOAM release notes for version 1.5*
 #+AUTHOR:                      OpenCFD Ltd.
-#+DATE:                       26 August 2008
+#+DATE:                       07 October 2008
 #+LINK:                  http://www.opencfd.co.uk
 #+OPTIONS: author:nil ^:{}
 
@@ -34,7 +34,7 @@
     and $WM_COMPILER_ARCH environment variables in the OpenFOAM-1.5/etc/bashrc
     (or cshrc) file.
   + The location of the installation of the compiler is controlled by the
-    $WM_COMPILER_INST environment variable in the OpenFOAM-1.5/etc/settings.sh
+    $compilerInstall environment variable in the OpenFOAM-1.5/etc/settings.sh
     (or settings.csh) file.
 
 * Developments to solvers (applications)

ReleaseNotes-1.5.html

@@ -49,7 +49,7 @@ quires the installation of Qt 4.3.? for ParaView-3 to run.
     $WM_COMPILER and $WM_COMPILER_ARCH environment variables in the
     OpenFOAM-1.5/etc/bashrc (or cshrc) file.
   * The location of the installation of the compiler is controlled by the
-    $WM_COMPILER_INST environment variable in the OpenFOAM-1.5/etc/settings.sh
+    $compilerInstall environment variable in the OpenFOAM-1.5/etc/settings.sh
     (or settings.csh) file.
 
 4 Developments to solvers (applications)
@@ -202,7 +202,7 @@ To aid common monitoring and post-processing activities.
   * Non-blocking, non-buffered, parallel transfers with potential scaling
     benefits for larger number of processors.
 
-Date: 26 August 2008
+Date: 07 October 2008
 
 HTML generated by org-mode 6.06b in emacs 23
 

applications/solvers/compressible/rhoCentralFoam/BCs/U/maxwellSlipUFvPatchVectorField.C

@@ -157,9 +157,8 @@ void maxwellSlipUFvPatchVectorField::updateCoeffs()
 
     if(thermalCreep_)
     {
-        const GeometricField<scalar, fvPatchField, volMesh>& vsfT =
-            this->db().objectRegistry::
-            lookupObject<GeometricField<scalar, fvPatchField, volMesh> >("T");
+        const volScalarField& vsfT =
+            this->db().objectRegistry::lookupObject<volScalarField>("T");
         label patchi = this->patch().index();
         const fvPatchScalarField& pT = vsfT.boundaryField()[patchi];
         Field<vector> gradpT = fvc::grad(vsfT)().boundaryField()[patchi];

applications/solvers/incompressible/boundaryFoam/boundaryFoam.C

@@ -26,8 +26,8 @@ Application
     boundaryFoam
 
 Description
-    Steady-state solver for 1D turbulent flow, typically to generate boundary 
-    layer conditions at an inlet, for use in a simulation. 
+    Steady-state solver for 1D turbulent flow, typically to generate boundary
+    layer conditions at an inlet, for use in a simulation.
 
     Boundary layer code to calculate the U, k and epsilon distributions.
     Used to create inlet boundary conditions for experimental comparisons
@@ -82,11 +82,14 @@ int main(int argc, char *argv[])
         U += (Ubar - UbarStar);
         gradP += (Ubar - UbarStar)/(1.0/UEqn.A())().weightedAverage(mesh.V());
 
+        label id = y.size() - 1;
+
         scalar wallShearStress =
-            flowDirection & turbulence->R()()[0] & wallNormal;
+            flowDirection & turbulence->R()()[id] & wallNormal;
 
         scalar yplusWall
-            = ::sqrt(mag(wallShearStress))*y[0]/laminarTransport.nu()()[0];
+//            = ::sqrt(mag(wallShearStress))*y[id]/laminarTransport.nu()()[id];
+            = ::sqrt(mag(wallShearStress))*y[id]/turbulence->nuEff()()[id];
 
         Info<< "Uncorrected Ubar = " << (flowDirection & UbarStar.value())<< tab
             << "pressure gradient = " << (flowDirection & gradP.value()) << tab

applications/solvers/incompressible/oodles/oodles.C

@@ -66,6 +66,9 @@ int main(int argc, char *argv[])
           + sgsModel->divDevBeff(U)
         );
 
+        // Optionally ensure diagonal-dominance of the momentum matrix
+        UEqn.relax();
+
         if (momentumPredictor)
         {
             solve(UEqn == -fvc::grad(p));

applications/solvers/incompressible/simpleFoam/Make/options

@@ -7,5 +7,4 @@ EXE_LIBS = \
     -lincompressibleRASModels \
     -lincompressibleTransportModels \
     -lfiniteVolume \
-    -lmeshTools \
-    /* $(LIB_WM_OPTIONS_DIR)/libfbsdmalloc.o */
+    -lmeshTools

applications/test/List/ListLoop.s

-	.file	"ListLoop.C"
-	.section	.ctors,"aw",@progbits
-	.align 4
-	.long	_GLOBAL__I__Z4funcRN4Foam4ListIdEERKS1_RKNS0_IiEES7_
-	.text
-	.align 2
-	.p2align 4,,15
-	.type	_Z41__static_initialization_and_destruction_0ii, @function
-_Z41__static_initialization_and_destruction_0ii:
-.LFB2550:
-	pushl	%ebp
-.LCFI0:
-	movl	%esp, %ebp
-.LCFI1:
-	pushl	%ebx
-.LCFI2:
-	call	__i686.get_pc_thunk.bx
-	addl	$_GLOBAL_OFFSET_TABLE_, %ebx
-	subl	$20, %esp
-.LCFI3:
-	decl	%eax
-	je	.L7
-.L5:
-	addl	$20, %esp
-	popl	%ebx
-	leave
-	ret
-	.p2align 4,,7
-.L7:
-	cmpl	$65535, %edx
-	jne	.L5
-	leal	_ZSt8__ioinit@GOTOFF(%ebx), %eax
-	movl	%eax, (%esp)
-	call	_ZNSt8ios_base4InitC1Ev@PLT
-	movl	__dso_handle@GOT(%ebx), %eax
-	movl	$0, 4(%esp)
-	movl	%eax, 8(%esp)
-	leal	__tcf_0@GOTOFF(%ebx), %eax
-	movl	%eax, (%esp)
-	call	__cxa_atexit@PLT
-	addl	$20, %esp
-	popl	%ebx
-	leave
-	ret
-.LFE2550:
-	.size	_Z41__static_initialization_and_destruction_0ii, .-_Z41__static_initialization_and_destruction_0ii
-.globl __gxx_personality_v0
-	.align 2
-	.p2align 4,,15
-	.type	_GLOBAL__I__Z4funcRN4Foam4ListIdEERKS1_RKNS0_IiEES7_, @function
-_GLOBAL__I__Z4funcRN4Foam4ListIdEERKS1_RKNS0_IiEES7_:
-.LFB2552:
-	pushl	%ebp
-.LCFI4:
-	movl	$65535, %edx
-	movl	$1, %eax
-	movl	%esp, %ebp
-.LCFI5:
-	leave
-	jmp	_Z41__static_initialization_and_destruction_0ii
-.LFE2552:
-	.size	_GLOBAL__I__Z4funcRN4Foam4ListIdEERKS1_RKNS0_IiEES7_, .-_GLOBAL__I__Z4funcRN4Foam4ListIdEERKS1_RKNS0_IiEES7_
-	.align 2
-	.p2align 4,,15
-	.type	__tcf_0, @function
-__tcf_0:
-.LFB2551:
-	pushl	%ebp
-.LCFI6:
-	movl	%esp, %ebp
-.LCFI7:
-	pushl	%ebx
-.LCFI8:
-	call	__i686.get_pc_thunk.bx
-	addl	$_GLOBAL_OFFSET_TABLE_, %ebx
-	subl	$4, %esp
-.LCFI9:
-	leal	_ZSt8__ioinit@GOTOFF(%ebx), %eax
-	movl	%eax, (%esp)
-	call	_ZNSt8ios_base4InitD1Ev@PLT
-	addl	$4, %esp
-	popl	%ebx
-	leave
-	ret
-.LFE2551:
-	.size	__tcf_0, .-__tcf_0
-	.align 2
-	.p2align 4,,15
-.globl _Z4funcRN4Foam4ListIdEERKS1_RKNS0_IiEES7_
-	.type	_Z4funcRN4Foam4ListIdEERKS1_RKNS0_IiEES7_, @function
-_Z4funcRN4Foam4ListIdEERKS1_RKNS0_IiEES7_:
-.LFB2352:
-	pushl	%ebp
-.LCFI10:
-	movl	%esp, %ebp
-.LCFI11:
-	pushl	%edi
-.LCFI12:
-	pushl	%esi
-.LCFI13:
-	subl	$12, %esp
-.LCFI14:
-	movl	8(%ebp), %edx
-	movl	(%edx), %eax
-	testl	%eax, %eax
-	movl	%eax, -20(%ebp)
-	jle	.L16
-	movl	16(%ebp), %eax
-	movl	4(%edx), %edx
-	xorl	%ecx, %ecx
-	movl	4(%eax), %eax
-	movl	%edx, -12(%ebp)
-	movl	%eax, -16(%ebp)
-	movl	20(%ebp), %eax
-	movl	4(%eax), %edi
-	movl	12(%ebp), %eax
-	movl	4(%eax), %esi
-	.p2align 4,,7
-.L15:
-	movl	-16(%ebp), %edx
-	movl	(%edx,%ecx,4), %eax
-	movl	-12(%ebp), %edx
-	leal	(%edx,%eax,8), %eax
-	movl	(%edi,%ecx,4), %edx
-	incl	%ecx
-	cmpl	%ecx, -20(%ebp)
-	fldl	(%eax)
-	fsubl	(%esi,%edx,8)
-	fstpl	(%eax)
-	jne	.L15
-.L16:
-	addl	$12, %esp
-	popl	%esi
-	popl	%edi
-	leave
-	ret
-.LFE2352:
-	.size	_Z4funcRN4Foam4ListIdEERKS1_RKNS0_IiEES7_, .-_Z4funcRN4Foam4ListIdEERKS1_RKNS0_IiEES7_
-	.local	_ZSt8__ioinit
-	.comm	_ZSt8__ioinit,1,1
-	.section	.eh_frame,"a",@progbits
-.Lframe1:
-	.long	.LECIE1-.LSCIE1
-.LSCIE1:
-	.long	0x0
-	.byte	0x1
-	.string	"zPR"
-	.uleb128 0x1
-	.sleb128 -4
-	.byte	0x8
-	.uleb128 0x6
-	.byte	0x9b
-	.long	DW.ref.__gxx_personality_v0-.
-	.byte	0x1b
-	.byte	0xc
-	.uleb128 0x4
-	.uleb128 0x4
-	.byte	0x88
-	.uleb128 0x1
-	.align 4
-.LECIE1:
-.LSFDE1:
-	.long	.LEFDE1-.LASFDE1
-.LASFDE1:
-	.long	.LASFDE1-.Lframe1
-	.long	.LFB2550-.
-	.long	.LFE2550-.LFB2550
-	.uleb128 0x0
-	.byte	0x4
-	.long	.LCFI0-.LFB2550
-	.byte	0xe
-	.uleb128 0x8
-	.byte	0x85
-	.uleb128 0x2
-	.byte	0x4
-	.long	.LCFI1-.LCFI0
-	.byte	0xd
-	.uleb128 0x5
-	.byte	0x4
-	.long	.LCFI2-.LCFI1
-	.byte	0x83
-	.uleb128 0x3
-	.align 4
-.LEFDE1:
-.LSFDE5:
-	.long	.LEFDE5-.LASFDE5
-.LASFDE5:
-	.long	.LASFDE5-.Lframe1
-	.long	.LFB2551-.
-	.long	.LFE2551-.LFB2551
-	.uleb128 0x0
-	.byte	0x4
-	.long	.LCFI6-.LFB2551
-	.byte	0xe
-	.uleb128 0x8
-	.byte	0x85
-	.uleb128 0x2
-	.byte	0x4
-	.long	.LCFI7-.LCFI6
-	.byte	0xd
-	.uleb128 0x5
-	.byte	0x4
-	.long	.LCFI8-.LCFI7
-	.byte	0x83
-	.uleb128 0x3
-	.align 4
-.LEFDE5:
-	.hidden DW.ref.__gxx_personality_v0
-	.weak	DW.ref.__gxx_personality_v0
-	.section	.data.DW.ref.__gxx_personality_v0,"awG",@progbits,DW.ref.__gxx_personality_v0,comdat
-	.align 4
-	.type	DW.ref.__gxx_personality_v0, @object
-	.size	DW.ref.__gxx_personality_v0, 4
-DW.ref.__gxx_personality_v0:
-	.long	__gxx_personality_v0
-	.ident	"GCC: (GNU) 4.1.1"
-	.section	.text.__i686.get_pc_thunk.bx,"axG",@progbits,__i686.get_pc_thunk.bx,comdat
-.globl __i686.get_pc_thunk.bx
-	.hidden	__i686.get_pc_thunk.bx
-	.type	__i686.get_pc_thunk.bx, @function
-__i686.get_pc_thunk.bx:
-	movl	(%esp), %ebx
-	ret
-	.section	.note.GNU-stack,"",@progbits

applications/test/xfer/Make/files

+xferListTest.C
+
+/* EXE = $(FOAM_USER_APPBIN)/xferListTest */
+EXE = ./xferListTest

src/turbulenceModels/RAS/incompressible/kOmega/kOmegaWallViscosityI.H → applications/test/xfer/xferListTest.C

@@ -22,48 +22,66 @@ License
     along with OpenFOAM; if not, write to the Free Software Foundation,
     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
 
-Global
-    kOmegaWallViscosity
+Application
 
 Description
-    Calculate wall viscosity from wall-functions.
 
 \*---------------------------------------------------------------------------*/
 
-{
-    scalar Cmu25 = pow(Cmu_.value(), 0.25);
+#include "OSspecific.H"
+
+#include "IOstreams.H"
+#include "IStringStream.H"
+#include "labelList.H"
+
+using namespace Foam;
 
-    const fvPatchList& patches = mesh_.boundary();
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+//  Main program:
 
-    forAll(patches, patchi)
+int main(int argc, char *argv[])
+{
+    List<label> lstA(10);
+    List<label> lstC(IStringStream("(1 2 3 4)")());
+    
+    forAll(lstA, i)
     {
-        const fvPatch& curPatch = patches[patchi];
-
-        if (isType<wallFvPatch>(curPatch))
-        {
-            const scalarField& nuw = nu().boundaryField()[patchi];
-            scalarField& nutw = nut_.boundaryField()[patchi];
-
-            forAll(curPatch, facei)
-            {
-                label faceCelli = curPatch.faceCells()[facei];
-
-                scalar yPlus =
-                    Cmu25*y_[patchi][facei]*sqrt(k_[faceCelli])/nuw[facei];
-
-                if (yPlus > yPlusLam_)
-                {
-                    nutw[facei] =
-                         nuw[facei]
-                        *(yPlus*kappa_.value()/log(E_.value()*yPlus) - 1);
-                }
-                else
-                {
-                    nutw[facei] = 0.0;
-                }
-            }
-        }
+        lstA[i] = i;
     }
+        
+    Info<< "lstA: " << lstA << endl;
+    Info<< "lstC: " << lstC << endl;
+
+    xfer<List<label> > xA(lstA, true);
+    xfer<List<label> > xB;
+
+//    List<label> lstB( xferTmp(List<label>, lstC) );
+//    List<label> lstB( xfer<List<label> >(lstC) );
+    List<label> lstB( xA );
+
+    Info<< "xA: " << *xA << endl;
+    Info<< "xB: " << *xB << endl;
+    Info<< "lstA: " << lstA << endl;
+    Info<< "lstB: " << lstB << endl;
+    Info<< "lstC: " << lstC << endl;
+
+    xA = lstB;
+
+    Info<< "xA: " << *xA << endl;
+    Info<< "xB: " << *xB << endl;
+    Info<< "lstA: " << lstA << endl;
+    Info<< "lstB: " << lstB << endl;
+    Info<< "lstC: " << lstC << endl;
+
+    xB = xA;
+
+    Info<< "xA: " << *xA << endl;
+    Info<< "xB: " << *xB << endl;
+    Info<< "lstA: " << lstA << endl;
+    Info<< "lstB: " << lstB << endl;
+    Info<< "lstC: " << lstC << endl;
+
+    return 0;
 }
 
 

applications/utilities/mesh/advanced/collapseEdges/collapseEdges.C

@@ -462,6 +462,7 @@ int main(int argc, char *argv[])
 
 #   include "setRootCase.H"
 #   include "createTime.H"
+    runTime.functionObjects().off();
 #   include "createPolyMesh.H"
 
     scalar minLen(readScalar(IStringStream(args.additionalArgs()[0])()));

applications/utilities/mesh/advanced/combinePatchFaces/combinePatchFaces.C

@@ -439,6 +439,7 @@ int main(int argc, char *argv[])
 
 #   include "setRootCase.H"
 #   include "createTime.H"
+    runTime.functionObjects().off();
 #   include "createPolyMesh.H"
 
     scalar featureAngle(readScalar(IStringStream(args.additionalArgs()[0])()));

applications/utilities/mesh/advanced/modifyMesh/modifyMesh.C

@@ -332,6 +332,7 @@ int main(int argc, char *argv[])
 
 #   include "setRootCase.H"
 #   include "createTime.H"
+    runTime.functionObjects().off();
 #   include "createPolyMesh.H"
 
     bool overwrite = args.options().found("overwrite");

applications/utilities/mesh/advanced/refineHexMesh/refineHexMesh.C

@@ -56,6 +56,7 @@ int main(int argc, char *argv[])
     argList::validArgs.append("cellSet");
 #   include "setRootCase.H"
 #   include "createTime.H"
+    runTime.functionObjects().off();
 #   include "createMesh.H"
     pointMesh pMesh(mesh);
 

applications/utilities/mesh/advanced/refineWallLayer/refineWallLayer.C

@@ -54,6 +54,7 @@ int main(int argc, char *argv[])
 
 #   include "setRootCase.H"
 #   include "createTime.H"
+    runTime.functionObjects().off();
 #   include "createPolyMesh.H"
 
     word patchName(args.additionalArgs()[0]);

applications/utilities/mesh/advanced/removeFaces/removeFaces.C

@@ -53,6 +53,7 @@ int main(int argc, char *argv[])
 
 #   include "setRootCase.H"
 #   include "createTime.H"
+    runTime.functionObjects().off();
 #   include "createMesh.H"
 
     bool overwrite = args.options().found("overwrite");

applications/utilities/mesh/advanced/splitCells/splitCells.C

@@ -532,6 +532,7 @@ int main(int argc, char *argv[])
 
 #   include "setRootCase.H"
 #   include "createTime.H"
+    runTime.functionObjects().off();
 #   include "createPolyMesh.H"
 
     scalar featureAngle(readScalar(IStringStream(args.additionalArgs()[0])()));

applications/utilities/mesh/conversion/polyDualMesh/polyDualMeshApp.C

@@ -47,6 +47,7 @@ int main(int argc, char *argv[])
 
 #   include "setRootCase.H"
 #   include "createTime.H"
+    runTime.functionObjects().off();
 #   include "createPolyMesh.H"
 
     scalar featureAngle(readScalar(IStringStream(args.additionalArgs()[0])()));

applications/utilities/mesh/conversion/writeMeshObj/writeMeshObj.C

@@ -346,6 +346,7 @@ int main(int argc, char *argv[])
 #   include "addTimeOptions.H"
 #   include "setRootCase.H"
 #   include "createTime.H"
+    runTime.functionObjects().off();
 
     bool patchFaces = args.options().found("patchFaces");
     bool doCell = args.options().found("cell");

applications/utilities/mesh/generation/extrude2DMesh/doExtrude2DMesh.C

@@ -61,6 +61,7 @@ int main(int argc, char *argv[])
     argList::validOptions.insert("overwrite", "");
 #   include "setRootCase.H"
 #   include "createTime.H"
+    runTime.functionObjects().off();
 #   include "createPolyMesh.H"
 
     scalar thickness(readScalar(IStringStream(args.additionalArgs()[0])()));