31 template<
class Po
intField>
37 return (1.0/3.0)*(ps[0] + ps[1] + ps[2]);
59 const vector a = (pNext -
p)^(pAvg -
p);
76 const vector a = (pNext -
p)^(pAvg -
p);
79 const scalar an = a & sumAHat;
90 return (1.0/3.0)*sumAnc/sumAn;
99 template<
class Po
intField>
105 return 0.5*((ps[1] - ps[0])^(ps[2] - ps[0]));
127 const vector a = (pNext -
p)^(pAvg -
p);
166 pAvg += ps[operator[](
pi)];
167 fldAvg +=
fld[operator[](
pi)];
178 const point&
p = ps[operator[](
pi)];
179 const point& pNext = ps[operator[](fcIndex(
pi))];
181 const vector a = (pNext -
p)^(pAvg -
p);
193 const point&
p = ps[operator[](
pi)];
194 const point& pNext = ps[operator[](fcIndex(
pi))];
196 const vector a = (pNext -
p)^(pAvg -
p);
199 +
fld[operator[](fcIndex(
pi))]
202 const scalar an = a & sumAHat;
213 return (1.0/3.0)*sumAnf/sumAn;
#define forAll(list, i)
Loop across all elements in list.
void size(const label)
Override size to be inconsistent with allocated storage.
label fcIndex(const label i) const
Return the forward circular index, i.e. the next index.
static vector centre(const PointField &ps)
Return centre point given face points.
Type average(const pointField &, const Field< Type > &) const
Calculate average value at centroid of face.
static vector area(const PointField &ps)
Return vector area given face points.
gmvFile<< "tracers "<< particles.size()<< nl;{ pointField positions(particles.size());label particlei=0;forAllConstIter(Cloud< passiveParticle >, particles, iter) { positions[particlei++]=iter().position(mesh);} for(i=0;i< pTraits< point >::nComponents;i++) { forAll(positions, particlei) { gmvFile<< component(positions[particlei], i)<< ' ';} gmvFile<< nl;}}forAll(lagrangianScalarNames, i){ const word &name=lagrangianScalarNames[i];IOField< scalar > fld(IOobject(name, runTime.name(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
const dimensionedScalar c
Speed of light in a vacuum.
Form normalised(const VectorSpace< Form, Cmpt, Ncmpts > &vs)