nargs=1,
metavar=("N"),
type=int,
required=True,
help="number of initial refine steps",
)
parser.add_argument(
"--P",
nargs=1,
metavar=("P"),
type=int,
required=True,
help="number of and refine/smoothing steps",
)
args = parser.parse_args()
fid = open(args.input_file)
myFace = Face.from_ts(fid)
a = trimesh.Trimesh(vertices=myFace.vertex, faces=myFace.connect)
for i in range(args.N[0]):
a = a.subdivide()
for i in range(args.P[0]):
a = a.subdivide()
a = trimesh.smoothing.filter_taubin(a)
myFace = Face(a.vertices, a.faces)
basename, ext = os.path.splitext(args.input_file)
myFace.write(f"{basename}_refined_smooth_{args.N[0]}{ext}")
structured_grid.smooth(args.smooth[0])
structured_grid.generate_vertex()
structured_grid.generate_connect()
structured_grid.isolate_hole(args.hole)
if args.proj:
structured_grid.proj_vertex(args.proj[0])
if args.translate:
structured_grid.translate(args.translate)
basename, ext = os.path.splitext(args.output_file)
nsolid = max(structured_grid.solid_id) + 1
if nsolid == 1:
myFace = Face(structured_grid.vertex, structured_grid.connect)
if structured_grid.is_sparse:
myFace.reindex(structured_grid.vid_lookup)
myFace.write(f"{basename}{ext}")
else:
for sid in range(nsolid):
idtr = np.where(structured_grid.solid_id == sid)[0]
aVid = np.unique(structured_grid.connect[idtr, :].flatten())
myFace = Face(vertex=None, connect=structured_grid.connect[idtr, :])
if structured_grid.is_sparse:
myFace.reindex(structured_grid.vid_lookup)
myFace.write(f"{basename}{sid}{ext}",
structured_grid.vertex,
write_full_vertex_array=False)
