diff --git a/applications/test/momentOfInertia/Test-momentOfInertia.C b/applications/test/momentOfInertia/Test-momentOfInertia.C index 22327f1bdd..82427abcf9 100644 --- a/applications/test/momentOfInertia/Test-momentOfInertia.C +++ b/applications/test/momentOfInertia/Test-momentOfInertia.C @@ -198,7 +198,7 @@ int main(int argc, char *argv[]) { const label cellI = args.optionLookupOrDefault("cell", 0); - tensorField mI = momentOfInertia::meshInertia(mesh); + tensorField mI(momentOfInertia::meshInertia(mesh)); tensor& J = mI[cellI]; diff --git a/applications/test/pointField/Test-PointField.C b/applications/test/pointField/Test-PointField.C index 1282c7ee22..3fe316e6af 100644 --- a/applications/test/pointField/Test-PointField.C +++ b/applications/test/pointField/Test-PointField.C @@ -60,7 +60,7 @@ int main(int argc, char *argv[]) pointPatchVectorField::calculatedType() ); - pointVectorField V = U + 2*U; + pointVectorField V(U + 2*U); Info<< "End\n" << endl; diff --git a/tutorials/basic/potentialFoam/cylinder/system/controlDict b/tutorials/basic/potentialFoam/cylinder/system/controlDict index 10797a6065..874de15927 100644 --- a/tutorials/basic/potentialFoam/cylinder/system/controlDict +++ b/tutorials/basic/potentialFoam/cylinder/system/controlDict @@ -94,13 +94,15 @@ functions Info<< "\nEvaluating analytical solution" << endl; - volVectorField centres = UA.mesh().C(); - volScalarField magCentres = mag(centres); - volScalarField theta = acos((centres & vector(1,0,0))/magCentres); + const volVectorField& centres = UA.mesh().C(); + volScalarField magCentres(mag(centres)); + volScalarField theta(acos((centres & vector(1,0,0))/magCentres)); - volVectorField cs2theta = + volVectorField cs2theta + ( cos(2*theta)*vector(1,0,0) - + sin(2*theta)*vector(0,1,0); + + sin(2*theta)*vector(0,1,0) + ); UA = uInfX*(dimensionedVector(vector(1,0,0)) - pow((radius/magCentres),2)*cs2theta); @@ -109,8 +111,10 @@ functions UA.write(); volScalarField error("error", mag(U-UA)/mag(UA)); + Info<<"Writing relative error in U to " << error.objectPath() << endl; + error.write(); #}; }