33template<
class Type,
template<
class>
class PatchField,
class GeoMesh>
34void Foam::rawTopoChangerFvMesh::setUnmappedValues
36 GeometricField<Type, PatchField, GeoMesh>&
fld,
37 const bitSet& mappedFace,
38 const GeometricField<Type, PatchField, GeoMesh>& baseFld
45 PatchField<Type>& fvp =
const_cast<PatchField<Type>&
>
47 fld.boundaryField()[patchi]
50 const label start = fvp.patch().start();
53 if (!mappedFace[start+i])
59 fvp[i] = baseFld.boundaryField()[patchi][i];
66template<
class Type,
template<
class>
class PatchField,
class GeoMesh>
67void Foam::rawTopoChangerFvMesh::zeroUnmappedValues
69 const bitSet& mappedFace
72 typedef GeometricField<Type, PatchField, GeoMesh> FieldType;
80 FieldType&
fld = lookupObjectRef<FieldType>(fldNames[i]);
98 dimensioned<Type>(
fld.dimensions(),
Zero)
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))
List< word > names(const UPtrList< T > &list, const UnaryMatchPredicate &matcher)
List< word > wordList
A List of words.
static constexpr const zero Zero
Global zero (0)
#define forAll(list, i)
Loop across all elements in list.