38 if (nPoints_ != sourceFld.
size())
41 <<
"Number of source points = " << nPoints_
42 <<
" number of values = " << sourceFld.
size()
59 fld[i] = sourceFld[verts[0]];
65 w[0]*sourceFld[verts[0]]
66 + w[1]*sourceFld[verts[1]];
72 w[0]*sourceFld[verts[0]]
73 + w[1]*sourceFld[verts[1]]
74 + w[2]*sourceFld[verts[2]];
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))
Generic templated field type.
A 1D vector of objects of type <T> with a fixed length <N>.
void size(const label n)
Older name for setAddressableSize.
bool interpolate() const noexcept
Same as isPointData()
A class for managing temporary objects.
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
errorManipArg< error, int > exit(error &err, const int errNo=1)
#define forAll(list, i)
Loop across all elements in list.