Commit 43ec00e5 authored by Mark Olesen's avatar Mark Olesen
Browse files

COMP: avoid ambiguous construct from tmp - test/ utils/

- missed these on the first sweep
parent a494ec35
......@@ -129,8 +129,8 @@ int main(int argc, char *argv[])
scalar eps = ::Foam::exp(-scalar(i + 1));
scalar x = xStart;
scalarField y = yStart;
scalarField dydx = dyStart;
scalarField y(yStart);
scalarField dydx(dyStart);
scalarField yScale(ode.nEqns(), 1.0);
scalar hEst = 0.6;
......@@ -148,7 +148,7 @@ int main(int argc, char *argv[])
scalar x = xStart;
scalar xEnd = x + 1.0;
scalarField y = yStart;
scalarField y(yStart);
scalarField yEnd(ode.nEqns());
yEnd[0] = ::Foam::j0(xEnd);
......
......@@ -41,18 +41,17 @@ int main(int argc, char *argv[])
# include "createTime.H"
# include "createMesh.H"
volScalarField fx = pow(mesh.C().component(vector::X), 2);
volScalarField fx(pow(mesh.C().component(vector::X), 2));
fx.write();
volScalarField gradx4 = fvc::grad(fx)().component(vector::X);
volScalarField gradx4(fvc::grad(fx)().component(vector::X));
gradx4.write();
//volVectorField curlC = fvc::curl(1.0*mesh.C());
//volVectorField curlC(fvc::curl(1.0*mesh.C()));
//curlC.write();
/*
surfaceScalarField xf = mesh.Cf().component(vector::X);
surfaceScalarField xf4 = pow(xf, 4);
surfaceScalarField xf(mesh.Cf().component(vector::X));
surfaceScalarField xf4(pow(xf, 4));
for (int i=1; i<xf4.size()-1; i++)
{
......
......@@ -47,10 +47,10 @@ int main()
x[i] = -3 + 0.06*i;
}
scalarField b = 0.5*(1.0 + erf(x));
scalarField c = 1.0 - b;
scalarField gradb = (1/::sqrt(constant::mathematical::pi))*exp(-sqr(x));
scalarField lapb = -2*x*gradb;
scalarField b(0.5*(1.0 + erf(x)));
scalarField c(1.0 - b);
scalarField gradb((1/::sqrt(constant::mathematical::pi))*exp(-sqr(x)));
scalarField lapb(-2*x*gradb);
r = lapb*b*c/(gradb*gradb);
......
......@@ -102,7 +102,7 @@ int main(int argc, char *argv[])
for (int corr=0; corr<nCorr; corr++)
{
volScalarField rAU = 1.0/UEqn.A();
volScalarField rAU(1.0/UEqn.A());
U = rAU*UEqn.H();
phi = (fvc::interpolate(U) & mesh.Sf())
......
/*
volTensorField gradU = fvc::grad(U);
volSymmTensorField D = symm(fvc::grad(U));
volTensorField Dprim = symm(fvc::grad(U - UMean));
volTensorField gradU(fvc::grad(U));
volSymmTensorField D(symm(fvc::grad(U)));
volTensorField Dprim(symm(fvc::grad(U - UMean)));
volScalarField prod = -((U - UMean)*(U - UMean)) && D;
volScalarField prod(-((U - UMean)*(U - UMean)) && D);
*/
/*
......
......@@ -112,7 +112,7 @@ int main(int argc, char *argv[])
label nInternalFaces = mesh.nInternalFaces();
vectorField unitAreas = mesh.faceAreas();
vectorField unitAreas(mesh.faceAreas());
unitAreas /= mag(unitAreas);
const polyPatchList& patches = mesh.boundaryMesh();
......
......@@ -209,7 +209,7 @@ int main(int argc, char *argv[])
// write bounding box corners
if (args.optionFound("blockMesh"))
{
pointField cornerPts = boundBox(surf.points()).corners();
pointField cornerPts(boundBox(surf.points()).corners());
Info<<"// blockMeshDict info" << nl;
......
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