isoSurface Class Reference

A surface formed by the iso value. After "Regularised Marching Tetrahedra: improved iso-surface extraction", G.M. Treece, R.W. Prager and A.H. Gee. More...

Inheritance diagram for isoSurface:
[legend]
Collaboration diagram for isoSurface:
[legend]

Public Types

enum  filterType
 Filtering type. More...
 
- Public Types inherited from isoSurfaceBase
enum  algorithmType : uint8_t { ALGO_POINT, ALGO_CELL, ALGO_TOPO }
 The algorithm types. More...
 
enum  filterType : uint8_t {
  NONE = 0, CELL, DIAGCELL, PARTIAL = CELL,
  FULL = DIAGCELL
}
 The filtering (regularization) to apply. More...
 
- Public Types inherited from MeshedSurface< face >
typedef face FaceType
 Face type used. More...
 

Public Member Functions

 TypeName ("isoSurface")
 Runtime type information. More...
 
 isoSurface (const volScalarField &cellValues, const scalarField &pointValues, const scalar iso, const isoSurfaceBase::filterType filter, const boundBox &bounds=boundBox::invertedBox, const scalar mergeTol=1e-6)
 Construct from cell values and point values. More...
 
 isoSurface (const volScalarField &cellValues, const scalarField &pointValues, const scalar iso, const bool regularise, const boundBox &bounds=boundBox::invertedBox, const scalar mergeTol=1e-6)
 Construct from cell values and point values. More...
 
template<class Type >
tmp< Field< Type > > interpolate (const GeometricField< Type, fvPatchField, volMesh > &cCoords, const Field< Type > &pCoords) const
 Interpolates cCoords, pCoords. More...
 
template<class Type >
Foam::tmp< Foam::SlicedGeometricField< Type, Foam::fvPatchField, Foam::slicedFvPatchField, Foam::volMesh > > adaptPatchFields (const GeometricField< Type, fvPatchField, volMesh > &fld) const
 
template<class Type >
Foam::label generateFaceTriPoints (const volScalarField &cVals, const scalarField &pVals, const GeometricField< Type, fvPatchField, volMesh > &cCoords, const Field< Type > &pCoords, const DynamicList< Type > &snappedPoints, const labelList &snappedCc, const labelList &snappedPoint, const label facei, const scalar neiVal, const Type &neiPt, const bool hasNeiSnap, const Type &neiSnapPt, DynamicList< Type > &triPoints, DynamicList< label > &triMeshCells) const
 
template<class Type >
Foam::tmp< Foam::Field< Type > > interpolate (const label nPoints, const labelList &triPointMergeMap, const labelList &interpolatedPoints, const List< FixedList< label, 3 >> &interpolatedOldPoints, const List< FixedList< scalar, 3 >> &interpolationWeights, const DynamicList< Type > &unmergedValues)
 
template<class Type >
Foam::tmp< Foam::Field< Type > > interpolate (const GeometricField< Type, fvPatchField, volMesh > &cCoords, const Field< Type > &pCoords) const
 
- Public Member Functions inherited from isoSurfaceBase
 isoSurfaceBase (const scalar iso, const boundBox &bounds=boundBox::invertedBox)
 Construct with iso value. More...
 
scalar isoValue () const
 The iso-value associated with the surface. More...
 
const labelListmeshCells () const
 For each face, the original cell in mesh. More...
 
labelListmeshCells ()
 For each face, the original cell in mesh. More...
 
- Public Member Functions inherited from MeshedSurface< face >
 MeshedSurface ()
 Construct null, an empty surface. More...
 
 MeshedSurface (const MeshedSurface &surf)
 Copy construct. More...
 
 MeshedSurface (const UnsortedMeshedSurface< face > &surf)
 Copy construct from an UnsortedMeshedSurface. More...
 
 MeshedSurface (MeshedSurface &&surf)
 Move construct. More...
 
 MeshedSurface (UnsortedMeshedSurface< face > &&surf)
 Move construct from an UnsortedMeshedSurface. More...
 
 MeshedSurface (const pointField &pointLst, const UList< face > &faceLst, const UList< surfZone > &zoneLst)
 Copy construct from components (points, faces, zones). More...
 
 MeshedSurface (pointField &&pointLst, List< face > &&faceLst, const UList< surfZone > &zoneLst)
 Move construct from components (points, faces). More...
 
 MeshedSurface (const pointField &pointLst, const UList< face > &faceLst, const labelUList &zoneSizes=labelUList(), const UList< word > &zoneNames=UList< word >())
 Copy construct from components (points, faces). More...
 
 MeshedSurface (pointField &&pointLst, List< face > &&faceLst, const labelUList &zoneSizes=labelUList(), const UList< word > &zoneNames=UList< word >())
 Move construct from components (points, faces). More...
 
 MeshedSurface (const polyBoundaryMesh &bMesh, const bool globalPoints=false)
 Construct from a boundary mesh with local points/faces. More...
 
 MeshedSurface (const surfMesh &mesh)
 Construct from a surfMesh. More...
 
 MeshedSurface (const fileName &name)
 Construct from file name (uses extension to determine type) More...
 
 MeshedSurface (const fileName &name, const word &ext)
 Construct from file name (uses extension to determine type) More...
 
 MeshedSurface (Istream &is)
 Construct from Istream. More...
 
 MeshedSurface (const Time &runTime, const word &surfName=word::null)
 Construct from database. More...
 
 MeshedSurface (const MeshedSurface< face > &surf)
 
 MeshedSurface (MeshedSurface< face > &&surf)
 
bool read (const fileName &name, const word &ext)
 Read from file. Chooses reader based on explicit extension. More...
 
virtual bool read (const fileName &name)
 Read from file. Chooses reader based on detected extension. More...
 
virtual void write (const fileName &name, const dictionary &options=dictionary::null) const
 Generic write routine. Chooses writer based on extension. More...
 
void write (const Time &t, const word &surfName=word::null) const
 Write to database. More...
 
 ClassName ("MeshedSurface")
 Runtime type information. More...
 
 declareRunTimeSelectionTable (autoPtr, MeshedSurface, fileExtension,(const fileName &name),(name))
 
virtual ~MeshedSurface ()
 Destructor. More...
 
 declareMemberFunctionSelectionTable (void, UnsortedMeshedSurface, write, fileExtension,(const fileName &name, const MeshedSurface< face > &surf, const dictionary &options),(name, surf, options))
 
label size () const
 The surface size is the number of faces. More...
 
const List< face > & surfFaces () const
 Return const access to the faces. More...
 
const surfZoneListsurfZones () const
 Const access to the surface zones. More...
 
const vectorFieldSf () const
 Face area vectors (normals) More...
 
const scalarFieldmagSf () const
 Face area magnitudes. More...
 
const vectorFieldCf () const
 Face centres. More...
 
virtual void clear ()
 Clear all storage. More...
 
virtual void addZones (const UList< surfZone > &, const bool cullEmpty=false)
 Add surface zones. More...
 
virtual void addZones (const labelUList &sizes, const UList< word > &names, const bool cullEmpty=false)
 Add surface zones. More...
 
virtual void addZones (const labelUList &sizes, const bool cullEmpty=false)
 Add surface zones. More...
 
bool addZonesToFaces ()
 Propagate zone information on face regions. More...
 
bool addZonesToFaces ()
 Specialization for labelledTri. More...
 
bool addZonesToFaces ()
 
virtual void removeZones ()
 Remove surface zones. More...
 
virtual void movePoints (const pointField &newPoints)
 Move points. More...
 
virtual void scalePoints (const scalar scaleFactor)
 Scale points. A non-positive factor is ignored. More...
 
virtual void cleanup (const bool verbose)
 Remove invalid faces. More...
 
virtual bool stitchFaces (const scalar tol=SMALL, const bool verbose=false)
 
virtual bool checkFaces (const bool verbose=false)
 
virtual label nTriangles () const
 Count number of triangles. More...
 
virtual label nTriangles (List< label > &faceMap) const
 Count number of triangles, returning a face map of original ids. More...
 
virtual label triangulate ()
 Triangulate in-place, returning the number of triangles added. More...
 
virtual label triangulate (List< label > &faceMap)
 Triangulate in-place, returning the number of triangles added. More...
 
MeshedSurface subsetMesh (const BoolListType &include, labelList &pointMap, labelList &faceMap) const
 Return new surface. More...
 
MeshedSurface subsetMesh (const bitSet &include) const
 Return new surface. More...
 
MeshedSurface subsetMesh (const labelHashSet &include) const
 Return new surface. More...
 
Foam::MeshedSurface< facesubsetMesh (const BoolListType &include, labelList &pointMap, labelList &faceMap) const
 
void swap (MeshedSurface< face > &surf)
 Swap contents. More...
 
void transfer (pointField &pointLst, List< face > &faceLst)
 Transfer the components. More...
 
void transfer (MeshedSurface< face > &surf)
 Transfer the contents of the argument and annul the argument. More...
 
void transfer (UnsortedMeshedSurface< face > &surf)
 Transfer the contents of the argument and annul the argument. More...
 
autoPtr< MeshedSurface< face > > releaseGeom ()
 Release (clear) geometry and return for reuse. More...
 
void swapFaces (List< face > &faces)
 Swap the stored faces. More...
 
void swapPoints (pointField &points)
 Swap the stored points. More...
 
void swapZones (surfZoneList &zones)
 Swap the stored zones. More...
 
void writeStats (Ostream &os) const
 
void operator= (const MeshedSurface< face > &surf)
 Copy assignment. More...
 
void operator= (MeshedSurface< face > &&surf)
 Move assignment. More...
 
 operator MeshedSurfaceProxy< face > () const
 Conversion operator to MeshedSurfaceProxy. More...
 

Friends

class isoSurfaceCell
 Declare friendship to share some functionality. More...
 
class isoSurfaceTopo
 

Additional Inherited Members

- Static Public Member Functions inherited from isoSurfaceBase
static filterType getFilterType (const dictionary &dict, const isoSurfaceBase::filterType deflt)
 Get 'regularise' as bool or enumeration. More...
 
- Static Public Member Functions inherited from MeshedSurface< face >
static void write (const fileName &name, const MeshedSurface< face > &surf, const dictionary &options=dictionary::null)
 Write to file, selecting writer based on its extension. More...
 
static void write (const fileName &name, const word &ext, const MeshedSurface< face > &surf, const dictionary &options=dictionary::null)
 Write to file, selecting writer based on the given extension. More...
 
static bool canRead (const fileName &name, bool verbose=false)
 Can we read this file format? More...
 
static bool canReadType (const word &ext, bool verbose=false)
 Can we read this file format? More...
 
static bool canWriteType (const word &ext, bool verbose=false)
 Can we write this file format? More...
 
static wordHashSet readTypes ()
 Known readable file-types. More...
 
static wordHashSet writeTypes ()
 Known writable file-types. More...
 
static autoPtr< MeshedSurfaceNew (const fileName &name, const word &ext)
 Select constructed from filename (explicit extension) More...
 
static autoPtr< MeshedSurfaceNew (const fileName &name)
 Select constructed from filename (implicit extension) More...
 
- Static Public Attributes inherited from isoSurfaceBase
static const Enum< algorithmTypealgorithmNames
 Names for the iso-surface algorithms. More...
 
static const Enum< filterTypefilterNames
 Names for the filtering types. More...
 
- Protected Types inherited from isoSurfaceBase
typedef meshedSurface MeshStorage
 
- Protected Member Functions inherited from MeshedSurface< face >
void transcribe (MeshedSurface< face > &surf)
 Transfer points/zones from 'face' to other other shapes. More...
 
void transcribe (MeshedSurface< face > &surf)
 
void transcribe (MeshedSurface< face > &surf)
 
void transcribe (MeshedSurface< face > &surf)
 
void checkZones ()
 Basic sanity check on zones. More...
 
pointFieldstoredPoints ()
 Non-const access to global points. More...
 
List< face > & storedFaces ()
 Non-const access to the faces. More...
 
surfZoneListstoredZones ()
 Non-const access to the zones. More...
 
void sortFacesAndStore (DynamicList< face > &unsortedFaces, DynamicList< label > &zoneIds, const bool sorted)
 Sort faces by zones and store sorted faces. More...
 
virtual void remapFaces (const labelUList &faceMap)
 Set new zones from faceMap. More...
 
- Protected Attributes inherited from isoSurfaceBase
const scalar iso_
 Iso value. More...
 
const boundBox bounds_
 Optional bounds. More...
 
labelList meshCells_
 For every face, the original cell in mesh. More...
 

Detailed Description

A surface formed by the iso value. After "Regularised Marching Tetrahedra: improved iso-surface extraction", G.M. Treece, R.W. Prager and A.H. Gee.

Note:

  • does tets without using cell centres/cell values. Not tested.
  • regularisation can create duplicate triangles/non-manifold surfaces. Current handling of those is bit ad-hoc for now and not perfect.
  • regularisation does not do boundary points so as to preserve the boundary perfectly.
  • uses geometric merge with fraction of bounding box as distance.
  • triangles can be between two cell centres so constant sampling does not make sense.
  • on empty patches behaves like zero gradient.
  • does not do 2D correctly, creates non-flat iso surface.
  • on processor boundaries might two overlapping (identical) triangles (one from either side)

The handling on coupled patches is a bit complex. All fields (values and coordinates) get rewritten so

  • empty patches get zerogradient (value) and facecentre (coordinates)
  • separated processor patch faces get interpolate (value) and facecentre (coordinates). (this is already the default for cyclics)
  • non-separated processor patch faces keep opposite value and cell centre

Now the triangle generation on non-separated processor patch faces can use the neighbouring value. Any separated processor face or cyclic face gets handled just like any boundary face.

Source files

Definition at line 89 of file isoSurface.H.

Member Enumeration Documentation

◆ filterType

enum filterType : uint8_t
strong

Filtering type.

Definition at line 87 of file isoSurfaceBase.H.

Constructor & Destructor Documentation

◆ isoSurface() [1/2]

isoSurface ( const volScalarField cellValues,
const scalarField pointValues,
const scalar  iso,
const isoSurfaceBase::filterType  filter,
const boundBox bounds = boundBox::invertedBox,
const scalar  mergeTol = 1e-6 
)

Construct from cell values and point values.

Uses boundaryField for boundary values. Holds reference to cellIsoVals and pointIsoVals.

Parameters
boundsoptional bounding box for trimming
mergeTolfraction of mesh bounding box for merging points

Definition at line 1339 of file isoSurface.C.

◆ isoSurface() [2/2]

isoSurface ( const volScalarField cellValues,
const scalarField pointValues,
const scalar  iso,
const bool  regularise,
const boundBox bounds = boundBox::invertedBox,
const scalar  mergeTol = 1e-6 
)

Construct from cell values and point values.

Uses boundaryField for boundary values. Holds reference to cellIsoVals and pointIsoVals.

Parameters
boundsoptional bounding box for trimming
mergeTolfraction of mesh bounding box for merging points

Definition at line 1361 of file isoSurface.C.

References triSurface::clearOut(), polyPatch::coupled(), Foam::expressions::patchExpr::debug, Foam::dimLength, Foam::dimless, Foam::endl(), f(), forAll, Foam::inplaceRenumber(), Foam::max(), Foam::min(), Foam::nl, IOobject::NO_READ, IOobject::NO_WRITE, patches, polyPatch::patchSlice(), fileName::path(), Foam::Pout, GeometricField< Type, PatchField, GeoMesh >::primitiveField(), GeometricField< Type, PatchField, GeoMesh >::primitiveFieldRef(), List< T >::set(), List< T >::setSize(), FixedList< T, N >::size(), polyPatch::start(), triSurface::swapPoints(), fvMesh::time(), Time::timeName(), List< T >::transfer(), triSurface::triFaceFaces(), triSurfaceTools::validTri(), Ostream::write(), and Foam::Zero.

Here is the call graph for this function:

Member Function Documentation

◆ TypeName()

TypeName ( "isoSurface"  )

Runtime type information.

◆ interpolate() [1/3]

tmp<Field<Type> > interpolate ( const GeometricField< Type, fvPatchField, volMesh > &  cCoords,
const Field< Type > &  pCoords 
) const

Interpolates cCoords, pCoords.

Uses the references to the original fields used to create the iso surface.

◆ adaptPatchFields()

Definition at line 47 of file isoSurfaceTemplates.C.

References f(), fld, forAll, mesh, Foam::New(), patches, fvPatchField< Type >::patchInternalField(), and fvPatchField< Type >::patchNeighbourField().

Here is the call graph for this function:

◆ generateFaceTriPoints()

Foam::label generateFaceTriPoints ( const volScalarField cVals,
const scalarField pVals,
const GeometricField< Type, fvPatchField, volMesh > &  cCoords,
const Field< Type > &  pCoords,
const DynamicList< Type > &  snappedPoints,
const labelList snappedCc,
const labelList snappedPoint,
const label  facei,
const scalar  neiVal,
const Type &  neiPt,
const bool  hasNeiSnap,
const Type &  neiSnapPt,
DynamicList< Type > &  triPoints,
DynamicList< label > &  triMeshCells 
) const

Definition at line 447 of file isoSurfaceTemplates.C.

References DynamicList< T, SizeMin >::append(), f(), forAll, FixedList< T, N >::size(), and Foam::Zero.

Here is the call graph for this function:

◆ interpolate() [2/3]

Foam::tmp<Foam::Field<Type> > interpolate ( const label  nPoints,
const labelList triPointMergeMap,
const labelList interpolatedPoints,
const List< FixedList< label, 3 >> &  interpolatedOldPoints,
const List< FixedList< scalar, 3 >> &  interpolationWeights,
const DynamicList< Type > &  unmergedValues 
)

Definition at line 748 of file isoSurfaceTemplates.C.

References forAll, Foam::New(), nPoints, Foam::HashTableOps::values(), and Foam::Zero.

Here is the call graph for this function:

◆ interpolate() [3/3]

Foam::tmp<Foam::Field<Type> > interpolate ( const GeometricField< Type, fvPatchField, volMesh > &  cCoords,
const Field< Type > &  pCoords 
) const

Definition at line 811 of file isoSurfaceTemplates.C.

References Foam::constant::physicoChemical::c2, Foam::interpolate(), and points.

Here is the call graph for this function:

Friends And Related Function Documentation

◆ isoSurfaceCell

friend class isoSurfaceCell
friend

Declare friendship to share some functionality.

Definition at line 385 of file isoSurface.H.

◆ isoSurfaceTopo

friend class isoSurfaceTopo
friend

Definition at line 386 of file isoSurface.H.


The documentation for this class was generated from the following files: