def main(): m = TetgenReader() base = 'tetgen_test_flipped.1' m.read_tetgen(base + '.node', base + '.smesh', base + '.ele', dimension_flag=3) m.write_nastran(base + '.bdf') m2 = TetgenReader() base = 'tetgen_test_flipped.1' m2.read_tetgen(base + '.node', base + '.smesh', base + '.ele', dimension_flag=2) ntris, four = m2.tris.shape #boundary_nodes = unique(m2.tris) #nboundary_nodes, = boundary_nodes.shape nboundary_nodes, three = m2.nodes.shape assert isinstance(nboundary_nodes, int), nboundary_nodes ntets, four = m.tets.shape nnodes, three = m.nodes.shape m.header = { 'inew': -1, 'nElements' : ntets, 'nPoints' : nnodes, 'nBoundPts' : nboundary_nodes, 'nViscPts' : nnodes, 'nViscElem' : ntets, 'tc' : 0.0, # d } write_usm3d_volume(m, base)
def main(): m = TetgenReader() base = 'tetgen_test_flipped.1' m.read_tetgen(base + '.node', base + '.smesh', base + '.ele', dimension_flag=3) m.write_nastran(base + '.bdf') m2 = TetgenReader() base = 'tetgen_test_flipped.1' m2.read_tetgen(base + '.node', base + '.smesh', base + '.ele', dimension_flag=2) ntris, four = m2.tris.shape #boundary_nodes = unique(m2.tris) #nboundary_nodes, = boundary_nodes.shape nboundary_nodes, three = m2.nodes.shape assert isinstance(nboundary_nodes, int), nboundary_nodes ntets, four = m.tets.shape nnodes, three = m.nodes.shape m.header = { 'inew': -1, 'nElements': ntets, 'nPoints': nnodes, 'nBoundPts': nboundary_nodes, 'nViscPts': nnodes, 'nViscElem': ntets, 'tc': 0.0, # d } write_usm3d_volume(m, base)
def tetgen_to_usm3d(base): m = Tetgen() m.read_tetgen(base + '.node', base + '.smesh', base + '.ele', dimension_flag=3) m.write_nastran(base + '.bdf') m2 = Tetgen() m2.read_tetgen(base + '.node', base + '.smesh', base + '.ele', dimension_flag=2) ntris = m2.tris.shape[0] nboundary_nodes = m2.nodes.shape[0] if not isinstance(nboundary_nodes, int): msg = 'nboundary_nodes=%r; type=%s must be an integer' % ( nboundary_nodes, type(nboundary_nodes)) raise TypeError(msg) ntets = m.tets.shape[0] nnodes = m.nodes.shape[0] #assert isinstance(nnodes, int), type(nnodes) m.header = { 'inew': -1, 'nElements': ntets, 'nPoints': nnodes, 'nBoundPts': nboundary_nodes, 'nViscPts': nnodes, 'nViscElem': ntets, 'tc': 0.0, # d } write_usm3d_volume(m, base)