*--* ex_triangulate_polygon
*--* Header Begin
*--* LAGriT Example Input file
* Carl Gable
* gable -at- lanl -dot- gov
*
*--*
*----------------------------------------------------
*--* Triangulate a polygon and then use refine/rivara to refine
*--* mesh to desired element size. Refine different parts of the
*--* mesh to different resolution.
*--*
*--* Also smooth and reconnect the mesh finishing with a call to
*--* recon/1 to eliminate all obtuse angles on boundaries.
*--*
*--*----------------------------------------------------
*--* Header End
*--* ex_triangulate_polygon
*--*
* Define the resolution of the point distribution.

define / REFINE1 / 0.05
define / REFINE2 / 0.025
define / REFINE3 / 0.0125

cmo / create / cmotri1 / / / tri

read / avs / polygon.inp / cmotri1
triangulate/counterclockwise
cmo / setatt / cmotri1 / imt / 1 0 0 / 1
cmo / setatt / cmotri1 / itetclr / 1 0 0 / 1
resetpts / itp
quality
dump / gmv / output_tri1.gmv /cmotri1
dump / avs / output_tri1.inp / cmotri1

refine/rivara///edge/1 0 0/REFINE1///inclusive
quality
dump / gmv / output_tri2.gmv / cmotri1
dump / avs / output_tri2.inp / cmotri1

smooth/position/esug/ 1 0 0

pset / p_gt / attribute / xic / 1 0 0 / 0.3 / gt
refine/rivara///edge/ pset get p_gt /REFINE2///inclusive
quality
dump / gmv / output_tri3.gmv / cmotri1
dump / avs / output_tri3.inp / cmotri1

pset / p_gt / attribute / xic / 1 0 0 / 0.6 / gt
refine/rivara///edge/ pset get p_gt /REFINE3///inclusive
quality
dump / gmv / output_tri4.gmv / cmotri1
dump / avs / output_tri4.inp / cmotri1

smooth/position/esug/ 1 0 0
smooth/position/esug/ 1 0 0
recon / 0
smooth/position/esug/ 1 0 0
smooth/position/esug/ 1 0 0
recon / 0
smooth/position/esug/ 1 0 0
smooth/position/esug/ 1 0 0
recon / 0
smooth/position/esug/ 1 0 0
recon / 0
smooth/position/esug/ 1 0 0
recon / 0
smooth/position/esug/ 1 0 0
recon / 0
smooth/position/esug/ 1 0 0
recon / 0
quality
*
dump / gmv / output_tri5.gmv / cmotri1
dump / avs / output_tri5.inp / cmotri1
recon / 1
quality
dump / gmv / output_tri6.gmv / cmotri1
dump / avs / output_tri6.inp / cmotri1
*
cmo/status
cmo/printatt//-all-/minmax
quality
finish