68 out[i] = in[i].second();
70 if (
mag(in[i].first()) > VSMALL)
72 out[i] /= in[i].
first();
87 out[i].second() = out[i].
first()*in[i];
98 x.first() +=
y.first();
99 x.second() +=
y.second();
111 getPoints(
fld, pfld);
115 pfld = vt.transformPosition(pfld);
119 pfld = vt.invTransformPosition(pfld);
122 setPoints(pfld,
fld);
135 operator()(vt, forward,
fld);
147 getPoints(
fld, pfld);
149 cpp.transformPosition(pfld);
151 setPoints(pfld,
fld);
164 <<
" does not correspond to the number of points in the mesh "
186 if (
fld.size() != meshPoints.
size())
189 <<
" does not correspond to the number of points supplied "
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))
A 1D array of objects of type <T>, where the size of the vector is known and used for subscript bound...
void setSize(const label n)
Alias for resize()
A 2-tuple for storing two objects of dissimilar types. The container is similar in purpose to std::pa...
A 1D vector of objects of type <T>, where the size of the vector is known and can be used for subscri...
T & first()
Return the first element of the list.
void size(const label n)
Older name for setAddressableSize.
The coupledPolyPatch is an abstract base class for patches that couple regions of the computational d...
friend Ostream & operator(Ostream &, const faMatrix< Type > &)
Mesh consisting of general polyhedral cells.
label nPoints() const noexcept
Number of mesh points.
Wrapper for position + weight to be used in e.g. averaging.
static void getPoints(const UList< weightedPosition > &in, List< point > &out)
Get points.
weightedPosition()
Construct null.
static void syncPoints(const polyMesh &mesh, List< weightedPosition > &)
Synchronisation for mesh point positions.
static void setPoints(const UList< point > &in, List< weightedPosition > &out)
Set points.
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
gmvFile<< "tracers "<< particles.size()<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().x()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().y()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().z()<< " ";}gmvFile<< nl;forAll(lagrangianScalarNames, i){ word name=lagrangianScalarNames[i];IOField< scalar > s(IOobject(name, runTime.timeName(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
errorManipArg< error, int > exit(error &err, const int errNo=1)
#define forAll(list, i)
Loop across all elements in list.
Combine HashSet operation. Equivalent to 'a |= b'.
A non-counting (dummy) refCount.