40 auto&
fld = tfld.ref();
44 const edge& verts = pointToVerts_[i];
50 label idx = verts.
first();
69 label idx = verts.
second();
70 if (idx == verts.
first())
91 const scalar d = s1-s0;
94 const scalar
s = (
iso_-s0)/d;
95 val =
s*v1+(1.0-
s)*v0;
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.
T & first() noexcept
The first element of the list, position [0].
const T & second() const noexcept
Return second element, which is also the last element.
void size(const label n)
Older name for setAddressableSize.
An edge is a list of two point labels. The functionality it provides supports the discretisation on a...
static autoPtr< isoSurfaceBase > New(const isoSurfaceParams ¶ms, const volScalarField &cellValues, const scalarField &pointValues, const scalar iso, const bitSet &ignoreCells=bitSet())
Create for specified algorithm type.
const scalarField & pVals_
Point values.
const scalar iso_
Iso value.
const polyMesh & mesh_
Reference to mesh.
const scalarField & cVals_
Cell values.
tmp< Field< Type > > interpolateTemplate(const Field< Type > &cellData, const Field< Type > &pointData) const
Interpolates cellData and pointData fields.
label nPoints() const noexcept
Number of mesh points.
A class for managing temporary objects.
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)
#define forAll(list, i)
Loop across all elements in list.