40 label triI = triFaces.
size();
47 return split(SPLITTRIANGLE,
points, triI, quadI, triFaces, quadFaces);
78 for (label pI=0; pI<
nPoints; pI++)
80 centrePoint += meshPoints[operator[](pI)];
81 cf +=
fld[operator[](pI)];
90 for (label pI=0; pI<
nPoints; pI++)
103 (meshPoints[
operator[](pI)] - centrePoint)
104 ^ (meshPoints[
operator[]((pI + 1) %
nPoints)] - centrePoint)
113 return sumAf/(3*sumA);
Info<< nl<< "Wrote faMesh in vtk format: "<< writer.output().name()<< nl;}{ vtk::lineWriter writer(aMesh.points(), aMesh.edges(), fileName(aMesh.mesh().time().globalPath()/"finiteArea-edges"));writer.writeGeometry();writer.beginCellData(4);writer.writeProcIDs();{ Field< scalar > fld(faMeshTools::flattenEdgeField(aMesh.magLe(), true))
virtual void average()
Calculate the average.
A 1D vector of objects of type <T> that resizes itself as necessary to accept the new objects.
void setSize(const label n)
Same as resize()
Generic templated field type.
A 1D vector of objects of type <T>, where the size of the vector is known and can be used for subscri...
void size(const label n)
Older name for setAddressableSize.
const DynamicList< triPoints > triangles() const
Const access to the triangulation.
label nTriangles() const
Number of triangles after splitting.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
static constexpr const zero Zero
Global zero (0)