15 #include "../config/config.hpp"
16 #include "../general/stable3d.hpp"
17 #include "../general/globals.hpp"
23 #include "../fem/eltrans.hpp"
24 #include "../fem/coefficient.hpp"
25 #include "../general/zstr.hpp"
26 #ifdef MFEM_USE_ADIOS2
27 #include "../general/adios2stream.hpp"
36 class GeometricFactors;
37 class FaceGeometricFactors;
40 class FiniteElementSpace;
61 #ifdef MFEM_USE_ADIOS2
248 #ifdef MFEM_USE_MEMALLOC
303 void ReadMFEMMesh(std::istream &input,
int version,
int &curved);
312 int &curved,
int &read_gf,
bool &finalize_topo);
313 void ReadVTKMesh(std::istream &input,
int &curved,
int &read_gf,
314 bool &finalize_topo);
316 bool &finalize_topo,
const std::string &xml_prefix=
"");
317 void ReadNURBSMesh(std::istream &input,
int &curved,
int &read_gf);
318 void ReadInlineMesh(std::istream &input,
bool generate_edges =
false);
319 void ReadGmshMesh(std::istream &input,
int &curved,
int &read_gf);
321 #ifdef MFEM_USE_NETCDF
322 void ReadCubit(
const char *filename,
int &curved,
int &read_gf);
355 int *edge1,
int *edge2,
int *middle)
361 int *edge1,
int *edge2,
int *middle)
362 {
Bisection(i, v_to_v, edge1, edge2, middle); }
399 bool update_nodes =
true);
416 double threshold,
int nc_limit = 0,
420 const int *fine,
int nfine,
int op);
457 const FaceInfo &fi,
bool is_ghost);
485 int v0,
int v1,
int v2);
488 int v0,
int v1,
int v2,
int v3);
503 void InitMesh(
int Dim_,
int spaceDim_,
int NVert,
int NElem,
int NBdrElem);
508 void Loader(std::istream &input,
int generate_edges = 0,
509 std::string parse_tag =
"");
515 std::string section_delimiter =
"")
const;
523 double sx,
double sy,
double sz,
bool sfc_ordering);
532 bool generate_edges,
bool sfc_ordering);
535 void Make1D(
int n,
double sx = 1.0);
564 explicit Mesh(
const Mesh &mesh,
bool copy_nodes =
true);
585 int generate_edges = 0,
int refine = 1,
586 bool fix_orientation =
true);
597 int nx,
int ny,
Element::Type type,
bool generate_edges =
false,
598 double sx = 1.0,
double sy = 1.0,
bool sfc_ordering =
true);
606 double sx = 1.0,
double sy = 1.0,
double sz = 1.0,
607 bool sfc_ordering =
true);
671 const std::vector<Vector> &translations,
double tol = 1e-8)
const;
685 int *element_attributes,
int num_elements,
687 int *boundary_attributes,
int num_boundary_elements,
688 int dimension,
int space_dimension = -1);
692 Mesh(
int Dim_,
int NVert,
int NElem,
int NBdrElem = 0,
int spaceDim_ = -1)
694 if (spaceDim_ == -1) { spaceDim_ = Dim_; }
695 InitMesh(Dim_, spaceDim_, NVert, NElem, NBdrElem);
706 int AddVertex(
double x,
double y = 0.0,
double z = 0.0);
714 int AddTriangle(
int v1,
int v2,
int v3,
int attr = 1);
718 int AddQuad(
int v1,
int v2,
int v3,
int v4,
int attr = 1);
719 int AddQuad(
const int *vi,
int attr = 1);
721 int AddTet(
int v1,
int v2,
int v3,
int v4,
int attr = 1);
722 int AddTet(
const int *vi,
int attr = 1);
724 int AddWedge(
int v1,
int v2,
int v3,
int v4,
int v5,
int v6,
int attr = 1);
725 int AddWedge(
const int *vi,
int attr = 1);
727 int AddPyramid(
int v1,
int v2,
int v3,
int v4,
int v5,
int attr = 1);
730 int AddHex(
int v1,
int v2,
int v3,
int v4,
int v5,
int v6,
int v7,
int v8,
732 int AddHex(
const int *vi,
int attr = 1);
747 int AddBdrQuad(
int v1,
int v2,
int v3,
int v4,
int attr = 1);
756 bool fix_orientation =
true);
759 bool fix_orientation =
true);
762 bool fix_orientation =
true);
765 bool fix_orientation =
true);
768 bool fix_orientation =
true);
771 void FinalizeMesh(
int refine = 0,
bool fix_orientation =
true);
800 virtual void Finalize(
bool refine =
false,
bool fix_orientation =
false);
824 int iterations = 4,
int window = 4,
825 int period = 2,
int seed = 0,
826 bool verbose =
false,
double time_limit = 0);
842 double sx = 1.0,
double sy = 1.0,
double sz = 1.0,
843 bool sfc_ordering =
true)
845 Make3D(nx, ny, nz, type, sx, sy, sz, sfc_ordering);
852 double sx = 1.0,
double sy = 1.0,
bool sfc_ordering =
true)
854 Make2D(nx, ny, type, sx, sy, generate_edges, sfc_ordering);
860 explicit Mesh(
int n,
double sx = 1.0)
869 explicit Mesh(
const char *filename,
int generate_edges = 0,
int refine = 1,
870 bool fix_orientation =
true);
875 explicit Mesh(std::istream &input,
int generate_edges = 0,
int refine = 1,
876 bool fix_orientation =
true);
879 Mesh(
Mesh *mesh_array[],
int num_pieces);
883 Mesh(
Mesh *orig_mesh,
int ref_factor,
int ref_type);
891 virtual void Load(std::istream &input,
int generate_edges = 0,
892 int refine = 1,
bool fix_orientation =
true)
894 Loader(input, generate_edges);
952 virtual long long ReduceInt(
int value)
const {
return value; }
1035 bool zerocopy =
false);
1052 return faces[i]->GetGeometryType();
1057 return elements[i]->GetGeometryType();
1062 return boundary[i]->GetGeometryType();
1135 faces[i]->GetVertices(vert);
1270 if (
faces_info[FaceNo].Elem2No < 0) {
return NULL; }
1451 void GetFaceInfos (
int Face,
int *Inf1,
int *Inf2)
const;
1452 void GetFaceInfos (
int Face,
int *Inf1,
int *Inf2,
int *NCFace)
const;
1530 void GetNode(
int i,
double *coord)
const;
1531 void SetNode(
int i,
const double *coord);
1575 virtual void SetCurvature(
int order,
bool discont =
false,
int space_dim = -1,
1600 int nonconforming = -1,
int nc_limit = 0);
1605 int nonconforming = -1,
int nc_limit = 0);
1609 int nonconforming = -1,
int nc_limit = 0);
1613 double eps = 0.0,
int nonconforming = -1);
1618 int nonconforming = -1,
int nc_limit = 0);
1623 int nonconforming = -1,
int nc_limit = 0);
1631 int nc_limit = 0,
int op = 1);
1635 int nc_limit = 0,
int op = 1);
1648 void EnsureNCMesh(
bool simplices_nonconforming =
false);
1675 virtual void Save(
const char *fname,
int precision=16)
const;
1678 #ifdef MFEM_USE_ADIOS2
1689 void PrintVTK(std::ostream &os,
int ref,
int field_data=0);
1697 bool high_order_output=
false,
1698 int compression_level=0,
1699 bool bdr_elements=
false);
1701 virtual void PrintVTU(std::string fname,
1703 bool high_order_output=
false,
1704 int compression_level=0,
1710 bool high_order_output=
false,
1711 int compression_level=0);
1720 std::ostream &os,
int elem_attr = 0)
const;
1724 int interior_faces = 0);
1761 double &kappa_min,
double &kappa_max,
1817 void Swap(
Mesh& other,
bool non_geometry);
1830 std::ostream &
operator<<(std::ostream &
out,
const Mesh &mesh);
1959 Mesh *
Extrude1D(Mesh *mesh,
const int ny,
const double sy,
1960 const bool closed =
false);
1963 Mesh *
Extrude2D(Mesh *mesh,
const int nz,
const double sz);
1969 a = c; c =
b; b =
t;
1973 std::ostream&
operator<<(std::ostream& os,
const Mesh::FaceInformation& info);
static Mesh MakePeriodic(const Mesh &orig_mesh, const std::vector< int > &v2v)
Create a periodic mesh by identifying vertices of orig_mesh.
Abstract class for all finite elements.
const IntegrationRule * IntRule
MFEM_DEPRECATED Mesh(int n, double sx=1.0)
Deprecated: see MakeCartesian1D.
void Loader(std::istream &input, int generate_edges=0, std::string parse_tag="")
double GetGeckoElementOrdering(Array< int > &ordering, int iterations=4, int window=4, int period=2, int seed=0, bool verbose=false, double time_limit=0)
Geometry::Constants< Geometry::PYRAMID > pyr_t
void PrintSurfaces(const Table &Aface_face, std::ostream &out) const
Print set of disjoint surfaces:
virtual void Save(const char *fname, int precision=16) const
void GetPointMatrix(int i, DenseMatrix &pointmat) const
static const int vtk_quadratic_hex[27]
int * CartesianPartitioning(int nxyz[])
FaceInformation GetFaceInformation(int f) const
Class for an integration rule - an Array of IntegrationPoint.
const DenseMatrix * PointMatrix
int GetBdrAttribute(int i) const
Return the attribute of boundary element i.
const double * GetVertex(int i) const
Return pointer to vertex i's coordinates.
void ScaleElements(double sf)
Class for grid function - Vector with associated FE space.
Table * GetEdgeVertexTable() const
Returns the edge-to-vertex Table (3D)
void GetBdrElementEdges(int i, Array< int > &edges, Array< int > &cor) const
Return the indices and the orientations of all edges of bdr element i.
void ReadVTKMesh(std::istream &input, int &curved, int &read_gf, bool &finalize_topo)
void FreeElement(Element *E)
int AddQuad(int v1, int v2, int v3, int v4, int attr=1)
int CheckElementOrientation(bool fix_it=true)
Check (and optionally attempt to fix) the orientation of the elements.
static bool remove_unused_vertices
void SetVertices(const Vector &vert_coord)
static const int vtk_quadratic_tet[10]
void GetEdgeVertices(int i, Array< int > &vert) const
Returns the indices of the vertices of edge i.
void Make2D(int nx, int ny, Element::Type type, double sx, double sy, bool generate_edges, bool sfc_ordering)
bool FaceIsTrueInterior(int FaceNo) const
Vector J
Jacobians of the element transformations at all quadrature points.
Vector X
Mapped (physical) coordinates of all quadrature points.
virtual void Eval(Vector &V, ElementTransformation &T, const IntegrationPoint &ip)
Evaluate the vector coefficient in the element described by T at the point ip, storing the result in ...
static Mesh MakeSimplicial(const Mesh &orig_mesh)
Base class for vector Coefficients that optionally depend on time and space.
void UniformRefinement3D_base(Array< int > *f2qf=NULL, DSTable *v_to_v_p=NULL, bool update_nodes=true)
const GeometricFactors * GetGeometricFactors(const IntegrationRule &ir, const int flags, MemoryType d_mt=MemoryType::DEFAULT)
Return the mesh geometric factors corresponding to the given integration rule.
static void PrintElementsByGeometry(int dim, const Array< int > &num_elems_by_geom, std::ostream &out)
Auxiliary method used by PrintCharacteristics().
void AddHexAsTets(const int *vi, int attr=1)
void AddHexAsWedges(const int *vi, int attr=1)
virtual void PrintInfo(std::ostream &os=mfem::out)
In serial, this method calls PrintCharacteristics(). In parallel, additional information about the pa...
virtual void Eval(Vector &V, ElementTransformation &T, const IntegrationPoint &ip)=0
Evaluate the vector coefficient in the element described by T at the point ip, storing the result in ...
Array< Element * > boundary
Geometry::Constants< Geometry::SQUARE > quad_t
int * GeneratePartitioning(int nparts, int part_method=1)
void GetBdrElementFace(int i, int *f, int *o) const
Return the index and the orientation of the face of bdr element i. (3D)
CoarseFineTransformations CoarseFineTr
int GetNumGeometries(int dim) const
Return the number of geometries of the given dimension present in the mesh.
virtual bool HasBoundaryElements() const
Checks if the mesh has boundary elements.
void MoveVertices(const Vector &displacements)
int GetNBE() const
Returns number of boundary elements.
IsoparametricTransformation Transformation
void GetLocalPtToSegTransformation(IsoparametricTransformation &, int)
Used in GetFaceElementTransformations (...)
void GetCharacteristics(double &h_min, double &h_max, double &kappa_min, double &kappa_max, Vector *Vh=NULL, Vector *Vk=NULL)
static FiniteElement * GetTransformationFEforElementType(Element::Type)
FaceGeometricFactors(const Mesh *mesh, const IntegrationRule &ir, int flags, FaceType type, MemoryType d_mt=MemoryType::DEFAULT)
static Mesh MakeCartesian3D(int nx, int ny, int nz, Element::Type type, double sx=1.0, double sy=1.0, double sz=1.0, bool sfc_ordering=true)
void NewNodes(GridFunction &nodes, bool make_owner=false)
Replace the internal node GridFunction with the given GridFunction.
virtual void UniformRefinement2D()
Refine a mixed 2D mesh uniformly.
void SwapNodes(GridFunction *&nodes, int &own_nodes_)
Element::Type GetElementType(int i) const
Returns the type of element i.
void GetBoundingBox(Vector &min, Vector &max, int ref=2)
Returns the minimum and maximum corners of the mesh bounding box.
void ReadNetgen2DMesh(std::istream &input, int &curved)
void ShiftRight(int &a, int &b, int &c)
Element::Type GetBdrElementType(int i) const
Returns the type of boundary element i.
static Mesh MakeRefined(Mesh &orig_mesh, int ref_factor, int ref_type)
Create a refined (by any factor) version of orig_mesh.
int EulerNumber2D() const
Equals 1 - num_holes.
void CreateVTKMesh(const Vector &points, const Array< int > &cell_data, const Array< int > &cell_offsets, const Array< int > &cell_types, const Array< int > &cell_attributes, int &curved, int &read_gf, bool &finalize_topo)
int AddTriangle(int v1, int v2, int v3, int attr=1)
Vector detJ
Determinants of the Jacobians at all quadrature points.
int AddBdrPoint(int v, int attr=1)
void GetGeometries(int dim, Array< Geometry::Type > &el_geoms) const
Return all element geometries of the given dimension present in the mesh.
void GetFaceVertices(int i, Array< int > &vert) const
Returns the indices of the vertices of face i.
Data type dense matrix using column-major storage.
Mesh * Extrude1D(Mesh *mesh, const int ny, const double sy, const bool closed)
Extrude a 1D mesh.
void Transform(void(*f)(const Vector &, Vector &))
void ApplyLocalSlaveTransformation(FaceElementTransformations &FT, const FaceInfo &fi, bool is_ghost)
int GetBdrElementEdgeIndex(int i) const
void AverageVertices(const int *indexes, int n, int result)
Averages the vertices with given indexes and saves the result in vertices[result].
int GetNE() const
Returns number of elements.
const Element * GetFace(int i) const
void FinalizeWedgeMesh(int generate_edges=0, int refine=0, bool fix_orientation=true)
Finalize the construction of a wedge Mesh.
Element * GetElement(int i)
void GenerateBoundaryElements()
Geometry::Type GetFaceBaseGeometry(int i) const
void PrintCharacteristics(Vector *Vh=NULL, Vector *Vk=NULL, std::ostream &os=mfem::out)
Compute and print mesh characteristics such as number of vertices, number of elements, number of boundary elements, minimal and maximal element sizes, minimal and maximal element aspect ratios, etc.
Structure for storing mesh geometric factors: coordinates, Jacobians, and determinants of the Jacobia...
NodeExtrudeCoefficient(const int dim, const int n_, const double s_)
void SetMeshGen()
Determine the mesh generator bitmask meshgen, see MeshGenerator().
static int GetQuadOrientation(const int *base, const int *test)
Returns the orientation of "test" relative to "base".
void GetVertices(Vector &vert_coord) const
void ReadNetgen3DMesh(std::istream &input)
Data arrays will be written in ASCII format.
void ReadInlineMesh(std::istream &input, bool generate_edges=false)
void GetBdrElementAdjacentElement2(int bdr_el, int &el, int &info) const
For the given boundary element, bdr_el, return its adjacent element and its info, i...
FaceElementTransformations * GetInteriorFaceTransformations(int FaceNo)
void RemoveInternalBoundaries()
void DebugDump(std::ostream &out) const
Output an NCMesh-compatible debug dump.
Mesh * Extrude2D(Mesh *mesh, const int nz, const double sz)
Extrude a 2D mesh.
virtual FaceElementTransformations * GetFaceElementTransformations(int FaceNo, int mask=31)
Geometry::Type GetElementBaseGeometry(int i) const
void FinalizeHexMesh(int generate_edges=0, int refine=0, bool fix_orientation=true)
Finalize the construction of a hexahedral Mesh.
Geometry::Type GetBdrElementGeometry(int i) const
void GetVertexToVertexTable(DSTable &) const
void RedRefinement(int i, const DSTable &v_to_v, int *edge1, int *edge2, int *middle)
void MarkTriMeshForRefinement()
void ReadCubit(const char *filename, int &curved, int &read_gf)
bool HasGeometry(Geometry::Type geom) const
Return true iff the given geom is encountered in the mesh. Geometries of dimensions lower than Dimens...
void InitRefinementTransforms()
IsoparametricTransformation FaceTransformation
void UniformRefinement2D_base(bool update_nodes=true)
Array< FaceGeometricFactors * > face_geom_factors
Optional face geometric factors.
Array< NCFaceInfo > nc_faces_info
Element * ReadElement(std::istream &)
void KnotInsert(Array< KnotVector * > &kv)
void GetLocalTriToPyrTransformation(IsoparametricTransformation &loc, int i)
Used in GetFaceElementTransformations (...)
A parallel extension of the NCMesh class.
std::vector< int > CreatePeriodicVertexMapping(const std::vector< Vector > &translations, double tol=1e-8) const
Creates a mapping v2v from the vertex indices of the mesh such that coincident vertices under the giv...
Element * NewElement(int geom)
FaceElementTransformations FaceElemTr
void GetEdgeTransformation(int i, IsoparametricTransformation *EdTr)
const GridFunction * GetNodes() const
Element::Type GetFaceElementType(int Face) const
void SetVerticesFromNodes(const GridFunction *nodes)
Helper to set vertex coordinates given a high-order curvature function.
Structure for storing face geometric factors: coordinates, Jacobians, determinants of the Jacobians...
ElementTransformation * GetBdrElementTransformation(int i)
Returns the transformation defining the i-th boundary element.
int AddBdrTriangle(int v1, int v2, int v3, int attr=1)
void PrintTopo(std::ostream &out, const Array< int > &e_to_k) const
int GetNumFaces() const
Return the number of faces (3D), edges (2D) or vertices (1D).
void MakeRefined_(Mesh &orig_mesh, const Array< int > ref_factors, int ref_type)
Internal function used in Mesh::MakeRefined.
Geometry::Constants< Geometry::SEGMENT > seg_t
void GetElementJacobian(int i, DenseMatrix &J)
static int GetTetOrientation(const int *base, const int *test)
Returns the orientation of "test" relative to "base".
MFEM_DEPRECATED Mesh(int nx, int ny, int nz, Element::Type type, bool generate_edges=false, double sx=1.0, double sy=1.0, double sz=1.0, bool sfc_ordering=true)
Deprecated: see MakeCartesian3D.
int AddTri(const int *vi, int attr=1)
int AddVertex(double x, double y=0.0, double z=0.0)
bool Nonconforming() const
Vector J
Jacobians of the element transformations at all quadrature points.
void MoveNodes(const Vector &displacements)
void GetEdgeOrdering(DSTable &v_to_v, Array< int > &order)
void GetElementData(int geom, Array< int > &elem_vtx, Array< int > &attr) const
void PrintWithPartitioning(int *partitioning, std::ostream &os, int elem_attr=0) const
Prints the mesh with boundary elements given by the boundary of the subdomains, so that the boundary ...
NCFaceInfo(bool slave, int master, const DenseMatrix *pm)
void FinalizeTriMesh(int generate_edges=0, int refine=0, bool fix_orientation=true)
Finalize the construction of a triangular Mesh.
double f(const Vector &xvec)
double GetElementSize(ElementTransformation *T, int type=0)
Symmetric 3D Table stored as an array of rows each of which has a stack of column, floor, number nodes. The number of the node is assigned by counting the nodes from zero as they are pushed into the table. Diagonals of any kind are not allowed so the row, column and floor must all be different for each node. Only one node is stored for all 6 symmetric entries that are indexable by unique triplets of row, column, and floor.
void PrintVTK(std::ostream &os)
void FinalizeTetMesh(int generate_edges=0, int refine=0, bool fix_orientation=true)
Finalize the construction of a tetrahedral Mesh.
void AddBdrQuadAsTriangles(const int *vi, int attr=1)
void Make3D(int nx, int ny, int nz, Element::Type type, double sx, double sy, double sz, bool sfc_ordering)
void AddSegmentFaceElement(int lf, int gf, int el, int v0, int v1)
virtual void SetAttributes()
void EnsureNCMesh(bool simplices_nonconforming=false)
GeometryList(Mesh &mesh, int dim)
Construct a GeometryList of all geometries of dimension dim in mesh.
Vector detJ
Determinants of the Jacobians at all quadrature points.
void AddHexAsPyramids(const int *vi, int attr=1)
int FindCoarseElement(int i)
GeometryList(Mesh &mesh)
Construct a GeometryList of all element geometries in mesh.
int AddTet(int v1, int v2, int v3, int v4, int attr=1)
void GetFaceTransformation(int i, IsoparametricTransformation *FTr)
Returns the transformation defining the given face element in a user-defined variable.
void GetLocalSegToTriTransformation(IsoparametricTransformation &loc, int i)
void UniformRefinement(int i, const DSTable &, int *, int *, int *)
void AddQuadFaceElement(int lf, int gf, int el, int v0, int v1, int v2, int v3)
void CheckPartitioning(int *partitioning_)
void LoadPatchTopo(std::istream &input, Array< int > &edge_to_knot)
Read NURBS patch/macro-element mesh.
void ReadGmshMesh(std::istream &input, int &curved, int &read_gf)
STable3D * GetElementToFaceTable(int ret_ftbl=0)
virtual void SetCurvature(int order, bool discont=false, int space_dim=-1, int ordering=1)
int MeshGenerator()
Get the mesh generator/type.
Type
Constants for the classes derived from Element.
void GetBdrPointMatrix(int i, DenseMatrix &pointmat) const
Geometry::Type GetFaceGeometry(int i) const
void CheckDisplacements(const Vector &displacements, double &tmax)
This structure stores the low level information necessary to interpret the configuration of elements ...
VTKFormat
Data array format for VTK and VTU files.
int AddBdrSegment(int v1, int v2, int attr=1)
Geometry::Constants< Geometry::TETRAHEDRON > tet_t
int AddElement(Element *elem)
The parameter elem should be allocated using the NewElement() method.
void ReadLineMesh(std::istream &input)
virtual void Load(std::istream &input, int generate_edges=0, int refine=1, bool fix_orientation=true)
void SetNodesOwner(bool nodes_owner)
Set the mesh nodes ownership flag.
void SetLayer(const int l)
int GetElementToEdgeTable(Table &, Array< int > &)
void GetBdrElementAdjacentElement(int bdr_el, int &el, int &info) const
For the given boundary element, bdr_el, return its adjacent element and its info, i...
void GetHilbertElementOrdering(Array< int > &ordering)
void RandomRefinement(double prob, bool aniso=false, int nonconforming=-1, int nc_limit=0)
Refine each element with given probability. Uses GeneralRefinement.
const Element * GetElement(int i) const
void GetLocalQuadToHexTransformation(IsoparametricTransformation &loc, int i)
Used in GetFaceElementTransformations (...)
static const int vtk_quadratic_pyramid[13]
void GetElementCenter(int i, Vector ¢er)
IsoparametricTransformation BdrTransformation
void GetLocalFaceTransformation(int face_type, int elem_type, IsoparametricTransformation &Transf, int info)
A helper method that constructs a transformation from the reference space of a face to the reference ...
IsoparametricTransformation Transformation2
virtual void SetNodalFESpace(FiniteElementSpace *nfes)
Element * GetBdrElement(int i)
Vector normal
Normals at all quadrature points.
A class for non-conforming AMR. The class is not used directly by the user, rather it is an extension...
FaceElementTransformations * GetBdrFaceTransformations(int BdrElemNo)
Geometry::Constants< Geometry::PRISM > pri_t
virtual MFEM_DEPRECATED void ReorientTetMesh()
void ReadMFEMMesh(std::istream &input, int version, int &curved)
virtual void LocalRefinement(const Array< int > &marked_el, int type=3)
This function is not public anymore. Use GeneralRefinement instead.
Data type tetrahedron element.
virtual void MarkTetMeshForRefinement(DSTable &v_to_v)
Array< Triple< int, int, int > > tmp_vertex_parents
MFEM_DEPRECATED Mesh(int nx, int ny, Element::Type type, bool generate_edges=false, double sx=1.0, double sy=1.0, bool sfc_ordering=true)
Deprecated: see MakeCartesian2D.
List of mesh geometries stored as Array<Geometry::Type>.
void PrintVTU(std::ostream &os, int ref=1, VTKFormat format=VTKFormat::ASCII, bool high_order_output=false, int compression_level=0, bool bdr_elements=false)
int CheckBdrElementOrientation(bool fix_it=true)
Check the orientation of the boundary elements.
int SpaceDimension() const
GeometricFactors(const Mesh *mesh, const IntegrationRule &ir, int flags, MemoryType d_mt=MemoryType::DEFAULT)
void PrintElementsWithPartitioning(int *partitioning, std::ostream &out, int interior_faces=0)
bool FaceIsInterior(int FaceNo) const
Return true if the given face is interior.
void GetLocalQuadToWdgTransformation(IsoparametricTransformation &loc, int i)
Used in GetFaceElementTransformations (...)
void AddPointFaceElement(int lf, int gf, int el)
Used in GenerateFaces()
void GetFaceElements(int Face, int *Elem1, int *Elem2) const
void Make1D(int n, double sx=1.0)
Creates a 1D mesh for the interval [0,sx] divided into n equal intervals.
const Element *const * GetElementsArray() const
Array< int > bdr_attributes
A list of all unique boundary attributes used by the Mesh.
virtual void NonconformingRefinement(const Array< Refinement > &refinements, int nc_limit=0)
This function is not public anymore. Use GeneralRefinement instead.
virtual ~Mesh()
Destroys Mesh.
void FinalizeMesh(int refine=0, bool fix_orientation=true)
Finalize the construction of any type of Mesh.
void GetFaceInfos(int Face, int *Inf1, int *Inf2) const
int AddSegment(int v1, int v2, int attr=1)
Geometry::Constants< Geometry::CUBE > hex_t
Class FiniteElementSpace - responsible for providing FEM view of the mesh, mainly managing the set of...
Geometry::Type geom_buf[Geometry::NumGeom]
static Mesh LoadFromFile(const char *filename, int generate_edges=0, int refine=1, bool fix_orientation=true)
void RemoveUnusedVertices()
Remove unused vertices and rebuild mesh connectivity.
static void PrintElement(const Element *, std::ostream &)
virtual void NURBSUniformRefinement()
Refine NURBS mesh.
int AddHex(int v1, int v2, int v3, int v4, int v5, int v6, int v7, int v8, int attr=1)
const IntegrationRule * IntRule
MemoryType
Memory types supported by MFEM.
Vector X
Mapped (physical) coordinates of all quadrature points.
void FinalizeTopology(bool generate_bdr=true)
Finalize the construction of the secondary topology (connectivity) data of a Mesh.
static void GetElementArrayEdgeTable(const Array< Element * > &elem_array, const DSTable &v_to_v, Table &el_to_edge)
GridFunction * GetNodes()
Return a pointer to the internal node GridFunction (may be NULL).
void SetSize(int nsize)
Change the logical size of the array, keep existing entries.
void PrepareNodeReorder(DSTable **old_v_to_v, Table **old_elem_vert)
int AddBdrQuad(int v1, int v2, int v3, int v4, int attr=1)
bool IsSlaveFace(const FaceInfo &fi) const
static Mesh MakeCartesian2D(int nx, int ny, Element::Type type, bool generate_edges=false, double sx=1.0, double sy=1.0, bool sfc_ordering=true)
void InitMesh(int Dim_, int spaceDim_, int NVert, int NElem, int NBdrElem)
Begin construction of a mesh.
void SetBdrAttribute(int i, int attr)
Set the attribute of boundary element i.
void ReadXML_VTKMesh(std::istream &input, int &curved, int &read_gf, bool &finalize_topo, const std::string &xml_prefix="")
Geometry::Type GetElementGeometry(int i) const
void RefineAtVertex(const Vertex &vert, double eps=0.0, int nonconforming=-1)
Refine elements sharing the specified vertex. Uses GeneralRefinement.
static Mesh MakeCartesian1D(int n, double sx=1.0)
Array< Element * > elements
const Table & ElementToElementTable()
double GetLength(int i, int j) const
Return the length of the segment from node i to node j.
const CoarseFineTransformations & GetRefinementTransforms()
bool RefineByError(const Array< double > &elem_error, double threshold, int nonconforming=-1, int nc_limit=0)
Operation GetLastOperation() const
Return type of last modification of the mesh.
Table * GetFaceToElementTable() const
void GetLocalTriToTetTransformation(IsoparametricTransformation &loc, int i)
Used in GetFaceElementTransformations (...)
int EulerNumber() const
Equals 1 + num_holes - num_loops.
NURBSExtension * NURBSext
Optional NURBS mesh extension.
int GetNV() const
Returns number of vertices. Vertices are only at the corners of elements, where you would expect them...
void PrintBdrVTU(std::string fname, VTKFormat format=VTKFormat::ASCII, bool high_order_output=false, int compression_level=0)
void GetNode(int i, double *coord) const
const Table & ElementToFaceTable() const
Class for integration point with weight.
void Swap(Mesh &other, bool non_geometry)
Element * ReadElementWithoutAttr(std::istream &)
int GetNumFacesWithGhost() const
Return the number of faces (3D), edges (2D) or vertices (1D) including ghost faces.
void GetElementEdges(int i, Array< int > &edges, Array< int > &cor) const
Return the indices and the orientations of all edges of element i.
static const int vtk_quadratic_wedge[18]
bool OwnsNodes() const
Return the mesh nodes ownership flag.
void GetElementTransformation(int i, IsoparametricTransformation *ElTr)
void GetLocalSegToQuadTransformation(IsoparametricTransformation &loc, int i)
const FiniteElementSpace * GetNodalFESpace() const
virtual void Print(std::ostream &os=mfem::out) const
void FinalizeQuadMesh(int generate_edges=0, int refine=0, bool fix_orientation=true)
Finalize the construction of a quadrilateral Mesh.
void GetBdrElementData(int geom, Array< int > &bdr_elem_vtx, Array< int > &bdr_attr) const
virtual void Finalize(bool refine=false, bool fix_orientation=false)
Finalize the construction of a general Mesh.
bool DerefineByError(Array< double > &elem_error, double threshold, int nc_limit=0, int op=1)
void ReorderElements(const Array< int > &ordering, bool reorder_vertices=true)
void GetFaceEdges(int i, Array< int > &edges, Array< int > &o) const
void GetLocalQuadToPyrTransformation(IsoparametricTransformation &loc, int i)
Used in GetFaceElementTransformations (...)
long long GetGlobalNE() const
Return the total (global) number of elements.
void DegreeElevate(int rel_degree, int degree=16)
virtual ~NodeExtrudeCoefficient()
constexpr int dimension
This example only works in 3D. Kernels for 2D are not implemented.
void MakeSimplicial_(const Mesh &orig_mesh, int *vglobal)
Array< FaceInfo > faces_info
STable3D * GetFacesTable()
std::ostream & operator<<(std::ostream &os, SparseMatrix const &mat)
NCMesh * ncmesh
Optional nonconforming mesh extension.
Geometry::Type GetBdrElementBaseGeometry(int i) const
virtual void UniformRefinement3D()
Refine a mixed 3D mesh uniformly.
int AddBdrElement(Element *elem)
static void PrintElementWithoutAttr(const Element *, std::ostream &)
Mesh(int Dim_, int NVert, int NElem, int NBdrElem=0, int spaceDim_=-1)
Init constructor: begin the construction of a Mesh object.
int AddWedge(int v1, int v2, int v3, int v4, int v5, int v6, int attr=1)
virtual bool NonconformingDerefinement(Array< double > &elem_error, double threshold, int nc_limit=0, int op=1)
NC version of GeneralDerefinement.
Geometry::Type GetFaceGeometryType(int Face) const
void Clear()
Clear the contents of the Mesh.
int GetNEdges() const
Return the number of edges.
Geometry::Constants< Geometry::TRIANGLE > tri_t
double AggregateError(const Array< double > &elem_error, const int *fine, int nfine, int op)
Derefinement helper.
double * GetVertex(int i)
Return pointer to vertex i's coordinates.
int GetNFaces() const
Return the number of faces in a 3D mesh.
virtual int FindPoints(DenseMatrix &point_mat, Array< int > &elem_ids, Array< IntegrationPoint > &ips, bool warn=true, InverseElementTransformation *inv_trans=NULL)
Find the ids of the elements that contain the given points, and their corresponding reference coordin...
void SetAttribute(int i, int attr)
Set the attribute of element i.
Mesh & operator=(Mesh &&mesh)
Move assignment operstor.
virtual void PrintXG(std::ostream &os=mfem::out) const
Print the mesh to the given stream using Netgen/Truegrid format.
void ReadTrueGridMesh(std::istream &input)
const FaceGeometricFactors * GetFaceGeometricFactors(const IntegrationRule &ir, const int flags, FaceType type, MemoryType d_mt=MemoryType::DEFAULT)
Return the mesh geometric factors for the faces corresponding to the given integration rule...
Array< int > GetFaceToBdrElMap() const
virtual int GetNFbyType(FaceType type) const
Returns the number of faces according to the requested type, does not count master nonconforming face...
static int GetTriOrientation(const int *base, const int *test)
Returns the orientation of "test" relative to "base".
void SetNode(int i, const double *coord)
virtual long long ReduceInt(int value) const
Utility function: sum integers from all processors (Allreduce).
int GetBdrFace(int BdrElemNo) const
Return the local face index for the given boundary face.
void GetLocalTriToWdgTransformation(IsoparametricTransformation &loc, int i)
Used in GetFaceElementTransformations (...)
void GetElementColoring(Array< int > &colors, int el0=0)
void ChangeVertexDataOwnership(double *vertices, int len_vertices, bool zerocopy=false)
Set the internal Vertex array to point to the given vertices array without assuming ownership of the ...
void DeleteGeometricFactors()
Destroy all GeometricFactors stored by the Mesh.
IsoparametricTransformation EdgeTransformation
MemAlloc< Tetrahedron, 1024 > TetMemory
Table * GetFaceEdgeTable() const
Returns the face-to-edge Table (3D)
void DoNodeReorder(DSTable *old_v_to_v, Table *old_elem_vert)
void SetNodes(const Vector &node_coord)
void SetNodalGridFunction(GridFunction *nodes, bool make_owner=false)
void ReadNURBSMesh(std::istream &input, int &curved, int &read_gf)
void ScaleSubdomains(double sf)
int AddPyramid(int v1, int v2, int v3, int v4, int v5, int attr=1)
OutStream out(std::cout)
Global stream used by the library for standard output. Initially it uses the same std::streambuf as s...
void InitFromNCMesh(const NCMesh &ncmesh)
Initialize vertices/elements/boundary/tables from a nonconforming mesh.
void Bisection(int i, const DSTable &, int *, int *, int *)
Bisect a triangle: element with index i is bisected.
void GeneralRefinement(const Array< Refinement > &refinements, int nonconforming=-1, int nc_limit=0)
Class for parallel meshes.
Abstract data type element.
int GetAttribute(int i) const
Return the attribute of element i.
void AddTriangleFaceElement(int lf, int gf, int el, int v0, int v1, int v2)
const Table & ElementToEdgeTable() const
Array< GeometricFactors * > geom_factors
Optional geometric factors.
void BdrBisection(int i, const HashTable< Hashed2 > &)
Bisect a boundary triangle: boundary element with index i is bisected.
void GetBdrElementVertices(int i, Array< int > &v) const
Returns the indices of the vertices of boundary element i.
void GenerateNCFaceInfo()
void GetElementVertices(int i, Array< int > &v) const
Returns the indices of the vertices of element i.
double GetElementVolume(int i)
Array< int > attributes
A list of all unique element attributes used by the Mesh.
Table * GetVertexToElementTable()
The returned Table must be destroyed by the caller.
const Element * GetBdrElement(int i) const
void GetElementData(const Array< Element * > &elem_array, int geom, Array< int > &elem_vtx, Array< int > &attr) const
void UpdateNodes()
Update the nodes of a curved mesh after refinement.
void GetElementFaces(int i, Array< int > &faces, Array< int > &ori) const
Return the indices and the orientations of all faces of element i.
void GreenRefinement(int i, const DSTable &v_to_v, int *edge1, int *edge2, int *middle)
void Printer(std::ostream &out=mfem::out, std::string section_delimiter="") const
void NodesUpdated()
This function should be called after the mesh node coordinates have changed, e.g. after the mesh has ...
Class used to extrude the nodes of a mesh.
void AddVertexParents(int i, int p1, int p2)
Mark vertex i as nonconforming, with parent vertices p1 and p2.