*--* ex_polygon_multi *--* Header Begin *--* LAGriT Example Input file * Carl Gable * gable -at- lanl -dot- gov * *--* *--* Triangulate a set of polygons, *--* Stitch the polygons together, *--* Refine and smooth the final 2D mesh. *--* Use MASSAGE to de-refine the mesh. *--* *--*---------------------------------------------------- *--* Header End *--* ex_polygon_multi *--* * * Triangulate the first polygon * cmo / create / cmotri1 / / / tri read / avs / polygon1.inp / cmotri1 define / ICOLOR / 1 * * Call another LAGriT input file which has the triangulation commands * infile / lagrit.input_tri_macro dump / gmv / output_tri1.gmv /cmotri1 ********** * * Triangulate the next polygon * cmo / create / cmotri2 / / / tri read / avs / polygon2.inp / cmotri2 define / ICOLOR / 2 infile / lagrit.input_tri_macro dump / gmv / output_tri2.gmv / cmotri2 dump / avs / output_tri2.inp / cmotri2 **** * * Triangulate the next polygon * cmo / create / cmotri3 / / / tri read / avs / polygon3.inp / cmotri3 define / ICOLOR / 3 infile / lagrit.input_tri_macro dump / gmv / output_tri3.gmv / cmotri3 * * * Triangulate the next polygon * cmo / create / cmotri4 / / / tri read / avs / polygon4.inp / cmotri4 define / ICOLOR / 4 infile / lagrit.input_tri_macro dump / gmv / output_tri4.gmv / cmotri4 * * Merge all the triangulated polygons into a single MO * addmesh / merge / cmo_all / cmotri1 / cmotri2 addmesh / merge / cmo_all / cmo_all / cmotri3 addmesh / merge / cmo_all / cmo_all / cmotri4 * * Figure out where there are duplicate points and remove duplicates * filter / 1 0 0 * * Recompute the jtet element adjancy information. * geniee * * Recompute the itp array * resetpts / itp * * Create parent/child points at material interfaces * settets * * Reset the imt values of nodes to be the same as the element itetclr * resetpts / itetclr dump / gmv / output_all.gmv / cmo_all * * Refine all triangles using edge mid-point bisection * infile lagrit.input_refine_smooth dump / gmv / output_smooth1.gmv / cmo_all * * Refine again * infile lagrit.input_refine_smooth dump / gmv / output_smooth2.gmv / cmo_all * * Refine again * infile lagrit.input_refine_smooth dump / gmv / output_smooth3.gmv / cmo_all * * Refine again * infile lagrit.input_refine_smooth dump / gmv / output_smooth4.gmv / cmo_all * * Modify the mesh using the MASSAGE command * pset / p_1 / attribute / imt / 1 0 0 / 1 / eq massage / 1.e20 / 0.08 / 0.001 / pset get p_1 / strictmergelength quality dump / gmv / output_massage1.gmv / cmo_all / ascii pset / p_2 / attribute / imt / 1 0 0 / 2 / eq massage / 1.e20 / 0.04 / 0.001 / pset get p_2 / strictmergelength quality dump / gmv / output_massage2.gmv / cmo_all / ascii pset / p_3 / attribute / imt / 1 0 0 / 3 / eq massage / 1.e20 / 0.02 / 0.001 / pset get p_3 / strictmergelength quality dump / gmv / output_massage3.gmv / cmo_all / ascii pset / p_4 / attribute / imt / 1 0 0 / 4 / eq massage / 1.e20 / 0.01 / 0.001 / pset get p_4 / strictmergelength quality dump / gmv / output_massage4.gmv / cmo_all / ascii infile lagrit.input_smooth_recon quality dump / gmv / output_massage5.gmv / cmo_all / ascii finish