def test_write_3x3_cube(verbose=0):
"""
Read sample openfoam mesh from aggelos and try to write in Ido's hex format
"""
xmf_base=os.path.join(os.path.dirname(os.path.abspath(__file__)),
'hex_cube_3x3')
h5_base = xmf_base
mesh_info = readMeshXdmf(xmf_base,h5_base,verbose=0)
writeHexMesh(mesh_info,'hexMesh_3x3',index_base=0)
def test_3x3_cube(verbose=0):
"""
Read sample openfoam mesh from aggelos and check that the basic information is correct
"""
xmf_base=os.path.join(os.path.dirname(os.path.abspath(__file__)),
'hex_cube_3x3')
h5_base = xmf_base
mesh_info = readMeshXdmf(xmf_base,h5_base,verbose=0)
eq(mesh_info.nElements_global,27)
eq(mesh_info.nNodes_global,64)
eq(mesh_info.nElements_global,mesh_info.nElements_owned)
eq(mesh_info.nNodes_global,mesh_info.nNodes_owned)
eq(mesh_info.nodeArray.shape,(mesh_info.nNodes_owned,3))
eq(mesh_info.elementNodesArray.shape,(mesh_info.nElements_owned,8))
eq(mesh_info.elementTopologyName,'Hexahedron')
eq(len(mesh_info.nodeMaterialTypes),mesh_info.nNodes_owned)
eq(len(mesh_info.elementMaterialTypes),mesh_info.nElements_owned)