Source code for Mordicus.Modules.Cemosis.IO.FeelppMeshReader

from Mordicus.Core.IO.MeshReaderBase import MeshReaderBase
import feelpp

def ReadMesh(meshFileName):
    """
    Functional API
    
    Readsthe mesh defined in the Feelpp mesh file "meshFileName" (.geo, .msh, .json)
    
    Parameters
    ----------
    meshFileName : str
        mesh file
        
    Returns
    -------
    FeelppMesh
        mesh of the model
    """
    reader = FeelppMeshReader(meshFileName=meshFileName)
    return reader.ReadMesh()

class FeelppMeshReader(MeshReaderBase):
    """
    Class containing a reader for Feelpp mesh file
    
    Attributes
    ----------
    meshFileName : str
        name of the mesh file (.geo, .msh, .json)
    dim : int
        dimension of the mesh
    gorder : int
        order of the mesh
    realdim : int
        real dimension of the mesh
    h : double
        hsize of the mesh (for geo)
    """

    def __init__(self, meshFileName, dim, gorder=1, realdim=None, h=0.1):
        """Feel++ Mesh Reader

        Args:
            meshFileName (_type_): mesh file name
            dim (int): topological dimension of the mesh. (1,2 or 3)
            gorder (int, optional): order of the mesh. Defaults to 1.
            realdim (int, optional): real dimension of the mesg. Defaults to dim.
            h (float, optional): mesh size. Defaults to 0.1.
        """        
        assert isinstance(meshFileName, str)
        assert feelpp.Environment.initialized

        super().__init__()
        self.meshFileName = meshFileName
        self.dim = dim
        self.gorder = gorder
        self.realdim = dim if realdim is None else realdim
        self.h = h

    def ReadMesh(self):
        """
        Read the mesh
        
        Returns
        -------
        FeelppMesh
            mesh of the model
        """
        from Mordicus.Modules.Cemosis.Containers.Meshes import FeelppUnstructuredMesh
        self.mesh = feelpp.load(feelpp.mesh(dim=self.dim, geo=self.gorder, realdim=self.realdim), self.meshFileName, self.h)
        if feelpp.Environment.isMasterRank():
            print("mesh ", self.mesh.dimension(), "D nelts:", self.mesh.numGlobalElements())
            print("mesh ", self.mesh.dimension(), "D nfaces:", self.mesh.numGlobalFaces())
            print("mesh ", self.mesh.dimension(), "D hmin:", self.mesh.hMin())
            print("mesh ", self.mesh.dimension(), "D havg:", self.mesh.hAverage())
            print("mesh ", self.mesh.dimension(), "D hmax:", self.mesh.hMax())
            print("mesh ", self.mesh.dimension(), "D measure:", self.mesh.measure())

        from Mordicus.Modules.Cemosis.Containers.Meshes import FeelppUnstructuredMesh as FUM
        return FUM.FeelppUnstructuredMesh(self.mesh)



if __name__ == "__main__":# pragma: no cover

    from Mordicus import RunTestFile
    RunTestFile(__file__)