def main(): m = Tetgen() base = 'tetgen_test_flipped.1' m.read_tetgen(base + '.node', base + '.smesh', base + '.ele', dimension_flag=3) m.write_nastran(base + '.bdf') m2 = Tetgen() base = 'tetgen_test_flipped.1' m2.read_tetgen(base + '.node', base + '.smesh', base + '.ele', dimension_flag=2) ntris = m2.tris.shape[0] #boundary_nodes = unique(m2.tris) #nboundary_nodes, = boundary_nodes.shape nboundary_nodes = m2.nodes.shape[0] assert isinstance(nboundary_nodes, int), nboundary_nodes ntets = m.tets.shape[0] nnodes = m.nodes[0] 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 = Tetgen() base = 'tetgen_test_flipped.1' m.read_tetgen(base + '.node', base + '.smesh', base + '.ele', dimension_flag=3) m.write_nastran(base + '.bdf') m2 = Tetgen() base = 'tetgen_test_flipped.1' m2.read_tetgen(base + '.node', base + '.smesh', base + '.ele', dimension_flag=2) ntris = m2.tris.shape[0] #boundary_nodes = unique(m2.tris) #nboundary_nodes, = boundary_nodes.shape nboundary_nodes = m2.nodes.shape[0] assert isinstance(nboundary_nodes, int), nboundary_nodes ntets = m.tets.shape[0] nnodes = m.nodes[0] m.header = { 'inew': -1, 'nElements' : ntets, 'nPoints' : nnodes, 'nBoundPts' : nboundary_nodes, 'nViscPts' : nnodes, 'nViscElem' : ntets, 'tc' : 0.0, # d } write_usm3d_volume(m, base)