Commit 5dc01ed6 authored by Mark Olesen's avatar Mark Olesen
Browse files

COMP: avoid ambiguous construct from tmp - utils/ error estimation

parent 54d02bd2
...@@ -112,7 +112,7 @@ int main(int argc, char *argv[]) ...@@ -112,7 +112,7 @@ int main(int argc, char *argv[])
); );
ee.residual()().write(); ee.residual()().write();
volScalarField e = ee.error(); volScalarField e(ee.error());
e.write(); e.write();
mag(e)().write(); mag(e)().write();
} }
......
...@@ -108,7 +108,7 @@ int main(int argc, char *argv[]) ...@@ -108,7 +108,7 @@ int main(int argc, char *argv[])
-fvc::grad(p) -fvc::grad(p)
); );
volVectorField e = ee.error(); volVectorField e(ee.error());
e.write(); e.write();
mag(e)().write(); mag(e)().write();
} }
......
...@@ -101,8 +101,8 @@ int main(int argc, char *argv[]) ...@@ -101,8 +101,8 @@ int main(int argc, char *argv[])
# include "createPhi.H" # include "createPhi.H"
volScalarField ek = 0.5*magSqr(U); volScalarField ek(0.5*magSqr(U));
volTensorField gradU = fvc::grad(U); volTensorField gradU(fvc::grad(U));
// Divergence of the error in U squared // Divergence of the error in U squared
......
...@@ -106,9 +106,9 @@ int main(int argc, char *argv[]) ...@@ -106,9 +106,9 @@ int main(int argc, char *argv[])
# include "createPhi.H" # include "createPhi.H"
volVectorField gradT = fvc::grad(T); volVectorField gradT(fvc::grad(T));
volScalarField TE = 0.5*sqr(T); volScalarField TE(0.5*sqr(T));
volScalarField L volScalarField L
( (
......
...@@ -57,13 +57,14 @@ div ...@@ -57,13 +57,14 @@ div
scalarField aNorm(vols.size(), 0.0); scalarField aNorm(vols.size(), 0.0);
// Get sign of flux // Get sign of flux
const surfaceScalarField signF = pos(flux); const surfaceScalarField signF(pos(flux));
// Calculate gradient of the solution // Calculate gradient of the solution
GeometricField GeometricField
< <
typename outerProduct<vector, Type>::type, fvPatchField, volMesh typename outerProduct<vector, Type>::type, fvPatchField, volMesh
> gradVf = fvc::grad(vf); >
gradVf(fvc::grad(vf));
// Internal faces // Internal faces
forAll(owner, faceI) forAll(owner, faceI)
......
...@@ -122,7 +122,7 @@ laplacian ...@@ -122,7 +122,7 @@ laplacian
const scalarField& vols = mesh.V(); const scalarField& vols = mesh.V();
const surfaceVectorField& Sf = mesh.Sf(); const surfaceVectorField& Sf = mesh.Sf();
const surfaceScalarField magSf = mesh.magSf(); const surfaceScalarField magSf(mesh.magSf());
const fvPatchList& patches = mesh.boundary(); const fvPatchList& patches = mesh.boundary();
const labelUList& owner = mesh.owner(); const labelUList& owner = mesh.owner();
const labelUList& neighbour = mesh.neighbour(); const labelUList& neighbour = mesh.neighbour();
...@@ -137,7 +137,8 @@ laplacian ...@@ -137,7 +137,8 @@ laplacian
GeometricField GeometricField
< <
typename outerProduct<vector, Type>::type, fvPatchField, volMesh typename outerProduct<vector, Type>::type, fvPatchField, volMesh
> gradVf = fvc::grad(vf); >
gradVf(fvc::grad(vf));
// Internal faces // Internal faces
forAll(owner, faceI) forAll(owner, faceI)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment