* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Jun 14 07:11:53 * * * Manual: http://lagrit.lanl.gov * * * * * * * * * * * * * * * * * * * * * * * * * * * * * -----oOo----- LaGriT V3 LACC-2012-084 LaGriT Copyright: This program was prepared by Los Alamos National Security, LLC at Los Alamos National Laboratory (LANL) under contract No. DE-AC52-06NA25396 with the U.S. Department of Energy (DOE). All rights in the program are reserved by the DOE and Los Alamos National Security, LLC. Permission is granted to the public to copy and use this software without charge, provided that this Notice and any statement of authorship are reproduced on all copies. Neither the U.S. Government nor LANS makes any warranty, express or implied, or assumes any liability or responsibility for the use of this software. -----oOo----- Output log file: outx3dgen Command log file: logx3dgen # Create 2D slice then translate to add each on Y axis # get point distribution from original hex mesh top read avs hex_top.inp mo geniee finish cmo/status/brief The current-mesh-object(CMO) is: mo 1 Mesh Object name: mo number of nodes = 3567 number of elements = 2240 dimensions geometry = 3 element type = hex dimensions topology = 3 8 nodes 6 faces 12 edges boundary flag = 16000000 status = active finish cmo printatt mo -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 3.000000000E+02 1.000000000E+03 7.000000000E+02 3567 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 3567 zic 5.500000000E+02 6.000000000E+02 5.000000000E+01 3567 # EXTRACT TRI SURFACE extract surfmesh 1,0,0 mos mo cmo/addatt/mos/itetclr0/vint/scalar/nelements////0 finish cmo/addatt/mos/itetclr1/vint/scalar/nelements////0 finish cmo/addatt/mos/idnode0/vint/scalar/nnodes////0 finish cmo/addatt/mos/idelem0/vint/scalar/nelements////0 finish cmo/addatt/mos/idelem1/vint/scalar/nelements////0 finish cmo/addatt/mos/facecol/vint/scalar/nelements////0 finish cmo/addatt/mos/idface0/vint/scalar/nelements////0 finish cmo/addatt/mos/idface1/vint/scalar/nelements////0 finish cmo/addatt/mos/ncon50/int/scalar/scalar/constant/permanent/x/0 finish cmo/addatt/mos/nconbnd/int/scalar/scalar/constant/permanent/x/0 finish cmo/addatt/mos/icontab/vint/scalar/ncon50/constant/permanent/x/0 finish resetpts/itp geniee finish finish cmo select mos settets normal assign_color_normal epsilon for zero: 1.0000000E-15 1392 out of 2512 element colors changed. dump gmv hex_top_surf.gmv mos cmo/modatt/-def-/-def-/ioflag/x finish cmo select mos eltset ebot itetclr eq 1 THE ELTSET ebot HAS 1120 ELEMENTS eltset eleft itetclr eq 5 THE ELTSET eleft HAS 80 ELEMENTS eltset ekeep union ebot eleft THE ELTSET ekeep HAS 1200 ELEMENTS eltset eduds not ekeep THE ELTSET eduds HAS 1312 ELEMENTS rmpoint element eltset,get,eduds Dudded 1243 points plus parents, and 1312 elements removed. rmpoint compress 1243 points removed and 0 elements removed. The new point count is 1271. The new element count is 1200. dump gmv tmp_ex_tri.gmv mos cmo/modatt/-def-/-def-/ioflag/x finish # EXTRACT LINE # reduce to a single set of points along x and z cmo setatt mos itetclr 1 1200 values reset for attribute itetclr resetpts itp geniee finish extract surfmesh 1,0,0 mol mos cmo/addatt/mol/itetclr0/vint/scalar/nelements////0 finish cmo/addatt/mol/itetclr1/vint/scalar/nelements////0 finish cmo/addatt/mol/idnode0/vint/scalar/nnodes////0 finish cmo/addatt/mol/idelem0/vint/scalar/nelements////0 finish cmo/addatt/mol/idelem1/vint/scalar/nelements////0 finish cmo/addatt/mol/facecol/vint/scalar/nelements////0 finish cmo/addatt/mol/idface0/vint/scalar/nelements////0 finish cmo/addatt/mol/idface1/vint/scalar/nelements////0 finish cmo/addatt/mol/idnode0/VINT/scalar/nnodes/linear/temporary/agltn CMO_ADDATT warning: attribute already exist: mol idnode0 finish cmo/addatt/mol/ncon50/int/scalar/scalar/constant/permanent/x/0 finish cmo/addatt/mol/nconbnd/int/scalar/scalar/constant/permanent/x/0 finish cmo/addatt/mol/icontab/vint/scalar/ncon50/constant/permanent/x/0 finish resetpts/itp geniee finish finish cmo delete mos Released Mesh Object: mos sort/mol/line_graph/ascending/ikey/ Default line_graph option: using "elements" cmo/addatt/mol/ikey/vint/scalar/nelements///gax/0 finish cmo/DELATT/mol/cid DELATT no action, attribute does not exist: cid finish cmo/addatt/mol/cid/VINT/scalar/nelements/linear/temporary finish cmo/DELATT/mol/ctype DELATT no action, attribute does not exist: ctype finish cmo/addatt/mol/ctype/VINT/scalar/nelements/linear/temporary finish cmo/DELATT/mol/loopid DELATT no action, attribute does not exist: loopid finish cmo/addatt/mol/loopid/VINT/scalar/nelements/linear/temporary finish SORT: order key written to attribute: ikey reorder/mol/ikey geniee finish dump gmv tmp_ex_line1.gmv mol cmo/modatt/-def-/-def-/ioflag/x finish pset/pyout/attribute yic/1,0,0/ne 0. THE PSET pyout HAS 109 POINTS pset/pxout/attribute xic/1,0,0/gt 700. THE PSET pxout HAS 63 POINTS pset/pout/union pyout pxout THE PSET pout HAS 121 POINTS rmpoint pset,get,pout Dudded 119 points plus parents, plus 120 elements removed. rmpoint compress 119 points removed and 0 elements removed. The new point count is 21. The new element count is 20. sort/mol/line_graph/ascending/ikey2/nodes cmo/addatt/mol/ikey2/vint/scalar/nnodes///gax/0 finish SORT: order key written to attribute: ikey2 reorder/mol/ikey2 dump gmv tmp_poly.gmv mol cmo/modatt/-def-/-def-/ioflag/x finish ### NOTE HARDWIRE REMOVAL OF EXTRA BITS ### cmo setatt mol itetclr/4,5,1/2 2 values reset for attribute itetclr rmmat 2 RMMAT DUDDED 0 POINTS: FLAGGED 2 ELEMENTS rmpoint compress 2 points removed and 2 elements removed. The new point count is 19. The new element count is 18. dump gmv top_poly.gmv mol cmo/modatt/-def-/-def-/ioflag/x finish dump avs2 top_poly_pts.inp mol 1 0 0 0 cmo/modatt/-def-/-def-/ioflag/x finish cmo delete mos Mesh Object release skipped, it does not exist: mos cmo delete mol Released Mesh Object: mol cmo delete mo Released Mesh Object: mo # Done with poly points for top triangle ############################################## # # cmo create motri///triplane read avs top_poly_pts.inp motri cmo/status/brief The current-mesh-object(CMO) is: motri 1 Mesh Object name: motri number of nodes = 19 number of elements = 0 dimensions geometry = 2 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = active finish # transpose y and z axis # since triangulate will not work in xz orientation cmo copyatt motri motri yic zic Warning: nelements = 0 in cmo motri Warning: nelements = 0 in cmo motri 19 copied from motri zic to -> motri yic cmo setatt motri zic 0. Warning: nelements = 0 19 values reset for attribute zic resetpts itp geniee geniee_cmo warning: no elements with faces so no jtet finish triangulate/counterclockwise filter/1,0,0 FILTER:Use internal epsilonl value = 0.895090418262E-10 FILTER:Dudding 0 duplicate points. finish rmpoint/compress 0 points removed and 0 elements removed. The new point count is 19. The new element count is 17. finish resetpts/itp geniee finish finish cmo/setatt/motri/imt/1 19 values reset for attribute imt1 cmo/setatt/motri/itetclr/1 17 values reset for attribute itetclr resetpts itp geniee finish dump/gmv/tmpxy_tri0.gmv/motri cmo/modatt/-def-/-def-/ioflag/x finish cmo copyatt motri motri zic yic 19 copied from motri yic to -> motri zic cmo setatt motri yic 0. 19 values reset for attribute yic resetpts itp geniee finish dump/gmv/tmpxz_tri0.gmv/motri cmo/modatt/-def-/-def-/ioflag/x finish massage/150./1.e-5/1.e-5/1,0,0/strictmergelength nodes in point set = 19 recon/0/0.10000000E-04/ RECON2D reconnection loop2d: 1 RECON: done. finish nodes in point set = 19 SGD: 0 total relaxations on point set of 19 indept. nodes. nodes in point set = 19 Face iteration, number of faces: 1 1 Face iteration, number of faces: 1 1 Face iteration, number of faces: 1 2 Face iteration, number of faces: 1 2 Face iteration, number of faces: 1 3 Face iteration, number of faces: 1 4 Face iteration, number of faces: 1 4 Face iteration, number of faces: 1 4 Face iteration, number of faces: 1 1 Face iteration, number of faces: 1 1 Face iteration, number of faces: 1 1 recon/0/0.10000000E-04/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 18 Recon2d: iteration number= 2 number of flips= 19 Recon2d: iteration number= 3 number of flips= 18 Recon2d: iteration number= 4 number of flips= 15 Recon2d: iteration number= 5 number of flips= 12 Recon2d: iteration number= 6 number of flips= 10 Recon2d: iteration number= 7 number of flips= 8 Recon2d: iteration number= 8 number of flips= 5 Recon2d: iteration number= 9 number of flips= 2 Recon2d: iteration number= 10 number of flips= 1 RECON: done. finish smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 43 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.3906250E+02 Max. area BEFORE smoothing= 0.1093750E+04 Iteration= 1; root mean square error= 0.3664630E+01 Iteration= 2; root mean square error= 0.2800341E+01 Iteration= 3; root mean square error= 0.2437500E+01 Iteration= 4; root mean square error= 0.2073102E+01 Iteration= 5; root mean square error= 0.1856393E+01 Iteration= 6; root mean square error= 0.1489907E+01 Iteration= 7; root mean square error= 0.1364325E+01 Iteration= 8; root mean square error= 0.1291844E+01 Iteration= 9; root mean square error= 0.8904487E+00 Iteration= 10; root mean square error= 0.5559862E+00 Iteration= 11; root mean square error= 0.3406355E+00 Iteration= 12; root mean square error= 0.2203083E+00 Iteration= 13; root mean square error= 0.1466666E+00 Iteration= 14; root mean square error= 0.9982798E-01 Iteration= 15; root mean square error= 0.6906332E-01 Iteration= 16; root mean square error= 0.4820909E-01 Iteration= 17; root mean square error= 0.3371242E-01 Min. area AFTER smoothing= 0.4683489E+02 Max. area AFTER smoothing= 0.3600475E+03 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 9 Recon2d: iteration number= 2 number of flips= 2 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 43 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.3306436E+02 Max. area BEFORE smoothing= 0.3600475E+03 Iteration= 1; root mean square error= 0.2236291E+01 Iteration= 2; root mean square error= 0.1463186E+01 Iteration= 3; root mean square error= 0.1143209E+01 Iteration= 4; root mean square error= 0.8358901E+00 Iteration= 5; root mean square error= 0.5780359E+00 Iteration= 6; root mean square error= 0.3833611E+00 Iteration= 7; root mean square error= 0.2542287E+00 Iteration= 8; root mean square error= 0.1677897E+00 Iteration= 9; root mean square error= 0.1107011E+00 Iteration= 10; root mean square error= 0.7329417E-01 Iteration= 11; root mean square error= 0.4875879E-01 Iteration= 12; root mean square error= 0.3252839E-01 Iteration= 13; root mean square error= 0.2170792E-01 Min. area AFTER smoothing= 0.4050915E+02 Max. area AFTER smoothing= 0.3927129E+03 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 3 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 43 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.2523021E+02 Max. area BEFORE smoothing= 0.3927129E+03 Iteration= 1; root mean square error= 0.1408345E+01 Iteration= 2; root mean square error= 0.7105313E+00 Iteration= 3; root mean square error= 0.4491723E+00 Iteration= 4; root mean square error= 0.3324129E+00 Iteration= 5; root mean square error= 0.2313966E+00 Iteration= 6; root mean square error= 0.1560290E+00 Iteration= 7; root mean square error= 0.1051469E+00 Iteration= 8; root mean square error= 0.7082071E-01 Iteration= 9; root mean square error= 0.4748951E-01 Iteration= 10; root mean square error= 0.3181415E-01 Iteration= 11; root mean square error= 0.2131859E-01 Iteration= 12; root mean square error= 0.1428936E-01 Iteration= 13; root mean square error= 0.9578319E-02 Min. area AFTER smoothing= 0.3422483E+02 Max. area AFTER smoothing= 0.3952063E+03 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 1 RECON: done. dump gmv tmp_tri01.gmv motri cmo/modatt/-def-/-def-/ioflag/x finish massage/100./1.e-5/1.e-5/1,0,0/strictmergelength nodes in point set = 43 recon/0/0.10000000E-04/ RECON2D reconnection loop2d: 1 RECON: done. finish nodes in point set = 43 SGD: 3 total relaxations on point set of 43 indept. nodes. rmpoint/compress 0 points removed and 0 elements removed. The new point count is 43. The new element count is 62. finish nodes in point set = 43 recon/0/0.10000000E-04/ RECON2D reconnection loop2d: 1 RECON: done. finish nodes in point set = 43 SGD: 0 total relaxations on point set of 43 indept. nodes. nodes in point set = 43 Face iteration, number of faces: 1 4 recon/0/0.10000000E-04/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 4 Recon2d: iteration number= 2 number of flips= 3 Recon2d: iteration number= 3 number of flips= 2 Recon2d: iteration number= 4 number of flips= 1 RECON: done. finish smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 47 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.3392444E+02 Max. area BEFORE smoothing= 0.5038008E+03 Iteration= 1; root mean square error= 0.1880656E+01 Iteration= 2; root mean square error= 0.1049119E+01 Iteration= 3; root mean square error= 0.6389711E+00 Iteration= 4; root mean square error= 0.2689473E+00 Iteration= 5; root mean square error= 0.1190588E+00 Iteration= 6; root mean square error= 0.5598946E-01 Iteration= 7; root mean square error= 0.2699164E-01 Iteration= 8; root mean square error= 0.1320330E-01 Min. area AFTER smoothing= 0.5383409E+02 Max. area AFTER smoothing= 0.3450371E+03 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 1 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 47 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.5383409E+02 Max. area BEFORE smoothing= 0.3450371E+03 Iteration= 1; root mean square error= 0.9276807E+00 Iteration= 2; root mean square error= 0.3468118E+00 Iteration= 3; root mean square error= 0.1626026E+00 Iteration= 4; root mean square error= 0.8861719E-01 Iteration= 5; root mean square error= 0.4660048E-01 Iteration= 6; root mean square error= 0.2537085E-01 Iteration= 7; root mean square error= 0.1435978E-01 Iteration= 8; root mean square error= 0.8286350E-02 Min. area AFTER smoothing= 0.5369288E+02 Max. area AFTER smoothing= 0.3451011E+03 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 1 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 47 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.5369288E+02 Max. area BEFORE smoothing= 0.3451011E+03 Iteration= 1; root mean square error= 0.7621089E+00 Iteration= 2; root mean square error= 0.1588434E+00 Iteration= 3; root mean square error= 0.6198920E-01 Iteration= 4; root mean square error= 0.3182651E-01 Iteration= 5; root mean square error= 0.1839822E-01 Iteration= 6; root mean square error= 0.1062443E-01 Iteration= 7; root mean square error= 0.6142976E-02 Min. area AFTER smoothing= 0.5369248E+02 Max. area AFTER smoothing= 0.3451691E+03 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 1 RECON: done. dump gmv tmp_tri02.gmv motri cmo/modatt/-def-/-def-/ioflag/x finish massage/50./1.e-5/1.e-5/1,0,0/strictmergelength nodes in point set = 47 recon/0/0.10000000E-04/ RECON2D reconnection loop2d: 1 RECON: done. finish nodes in point set = 47 SGD: 2 total relaxations on point set of 47 indept. nodes. rmpoint/compress 0 points removed and 0 elements removed. The new point count is 47. The new element count is 66. finish nodes in point set = 47 recon/0/0.10000000E-04/ RECON2D reconnection loop2d: 1 RECON: done. finish nodes in point set = 47 SGD: 0 total relaxations on point set of 47 indept. nodes. nodes in point set = 47 Face iteration, number of faces: 1 8 recon/0/0.10000000E-04/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 8 Recon2d: iteration number= 2 number of flips= 2 RECON: done. finish smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 55 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.3906250E+02 Max. area BEFORE smoothing= 0.2991819E+03 Iteration= 1; root mean square error= 0.1357812E+01 Iteration= 2; root mean square error= 0.3304494E+00 Iteration= 3; root mean square error= 0.1112627E+00 Iteration= 4; root mean square error= 0.5408884E-01 Iteration= 5; root mean square error= 0.2774065E-01 Iteration= 6; root mean square error= 0.1455994E-01 Iteration= 7; root mean square error= 0.7739773E-02 Min. area AFTER smoothing= 0.3906250E+02 Max. area AFTER smoothing= 0.2466099E+03 recon 0 RECON2D reconnection loop2d: 1 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 55 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.3906250E+02 Max. area BEFORE smoothing= 0.2466099E+03 Iteration= 1; root mean square error= 0.4144912E-02 Iteration= 2; root mean square error= 0.2222128E-02 Iteration= 3; root mean square error= 0.1187247E-02 Iteration= 4; root mean square error= 0.6316262E-03 Iteration= 5; root mean square error= 0.3347277E-03 Iteration= 6; root mean square error= 0.1768185E-03 Iteration= 7; root mean square error= 0.9316487E-04 Iteration= 8; root mean square error= 0.4898983E-04 Iteration= 9; root mean square error= 0.2572072E-04 Min. area AFTER smoothing= 0.3906250E+02 Max. area AFTER smoothing= 0.2466390E+03 recon 0 RECON2D reconnection loop2d: 1 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 55 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.3906250E+02 Max. area BEFORE smoothing= 0.2466390E+03 Iteration= 1; root mean square error= 0.1348768E-04 Iteration= 2; root mean square error= 0.7066227E-05 Iteration= 3; root mean square error= 0.3699369E-05 Iteration= 4; root mean square error= 0.1935662E-05 Iteration= 5; root mean square error= 0.1012391E-05 Iteration= 6; root mean square error= 0.5293314E-06 Iteration= 7; root mean square error= 0.2766940E-06 Iteration= 8; root mean square error= 0.1446073E-06 Iteration= 9; root mean square error= 0.7556453E-07 Min. area AFTER smoothing= 0.3906250E+02 Max. area AFTER smoothing= 0.2466391E+03 recon 0 RECON2D reconnection loop2d: 1 RECON: done. dump gmv tmp_tri03.gmv motri cmo/modatt/-def-/-def-/ioflag/x finish massage/25./12./1.e-5/1,0,0/strictmergelength nodes in point set = 55 rmpoint/element Dudded 0 points plus parents, and 6 elements removed. finish geniee finish Merged 3 nodes. rmpoint/element Dudded 0 points plus parents, and 2 elements removed. finish geniee finish Merged 1 nodes. recon/0/0.10000000E-04/ RECON2D reconnection loop2d: 1 RECON: done. finish nodes in point set = 51 SGD: 1 total relaxations on point set of 51 indept. nodes. rmpoint/compress 4 points removed and 0 elements removed. The new point count is 51. The new element count is 66. finish nodes in point set = 51 recon/0/0.10000000E-04/ RECON2D reconnection loop2d: 1 RECON: done. finish nodes in point set = 51 SGD: 0 total relaxations on point set of 51 indept. nodes. rmpoint/compress 0 points removed and 0 elements removed. The new point count is 51. The new element count is 66. finish nodes in point set = 51 recon/0/0.10000000E-04/ RECON2D reconnection loop2d: 1 RECON: done. finish nodes in point set = 51 SGD: 0 total relaxations on point set of 51 indept. nodes. nodes in point set = 51 Face iteration, number of faces: 1 15 Face iteration, number of faces: 1 6 Face iteration, number of faces: 1 1 recon/0/0.10000000E-04/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 2 RECON: done. finish smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 73 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.1953125E+02 Max. area BEFORE smoothing= 0.2450072E+03 Iteration= 1; root mean square error= 0.1073323E+01 Iteration= 2; root mean square error= 0.2378405E+00 Iteration= 3; root mean square error= 0.8803066E-01 Iteration= 4; root mean square error= 0.3851313E-01 Iteration= 5; root mean square error= 0.1896915E-01 Iteration= 6; root mean square error= 0.1009401E-01 Min. area AFTER smoothing= 0.1757812E+02 Max. area AFTER smoothing= 0.2289303E+03 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 2 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 73 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.1757812E+02 Max. area BEFORE smoothing= 0.2289303E+03 Iteration= 1; root mean square error= 0.6031523E+00 Iteration= 2; root mean square error= 0.1805008E+00 Iteration= 3; root mean square error= 0.4846648E-01 Iteration= 4; root mean square error= 0.1653041E-01 Iteration= 5; root mean square error= 0.6841464E-02 Iteration= 6; root mean square error= 0.3057351E-02 Min. area AFTER smoothing= 0.1674107E+02 Max. area AFTER smoothing= 0.2297997E+03 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 1 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 73 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.1674107E+02 Max. area BEFORE smoothing= 0.2297997E+03 Iteration= 1; root mean square error= 0.3661379E+00 Iteration= 2; root mean square error= 0.1234853E+00 Iteration= 3; root mean square error= 0.6680915E-01 Iteration= 4; root mean square error= 0.3654366E-01 Iteration= 5; root mean square error= 0.2004169E-01 Iteration= 6; root mean square error= 0.1101261E-01 Iteration= 7; root mean square error= 0.6058148E-02 Iteration= 8; root mean square error= 0.3330000E-02 Min. area AFTER smoothing= 0.1674107E+02 Max. area AFTER smoothing= 0.2217266E+03 recon 0 RECON2D reconnection loop2d: 1 RECON: done. dump gmv tmp_tri04.gmv motri cmo/modatt/-def-/-def-/ioflag/x finish # exclude join sides from node movements cmo addatt motri pmass VINT scalar nnodes infile mass_psave.mlgi pset/p1/attribute zic/1,0,0/eq 550. THE PSET p1 HAS 17 POINTS pset/p2/attribute xic/1,0,0/eq 300. THE PSET p2 HAS 3 POINTS pset/pstay/union p1 p2 THE PSET pstay HAS 19 POINTS cmo setatt motri pmass pset,get,pstay 1 19 values reset for attribute pmass pset/pmove/not pstay THE PSET pmove HAS 54 POINTS dump gmv tmp_pstay.gmv motri cmo/modatt/-def-/-def-/ioflag/x finish pset/p1/delete PSET p1 DELETED pset/p2/delete PSET p2 DELETED pset/pstay/delete PSET pstay DELETED finish dump gmv tmp_pstay.gmv motri cmo/modatt/-def-/-def-/ioflag/x finish cmo setatt motri itetclr 1 94 values reset for attribute itetclr cmo setatt motri imt 1 73 values reset for attribute imt1 resetpts itp geniee finish massage/25. 25. .1/pset,get,pmove/strictmergelength nodes in point set = 54 rmpoint/element Dudded 0 points plus parents, and 24 elements removed. finish geniee finish Merged 17 nodes. rmpoint/element Dudded 0 points plus parents, and 21 elements removed. finish geniee finish Merged 13 nodes. rmpoint/element Dudded 0 points plus parents, and 9 elements removed. finish geniee finish Merged 6 nodes. rmpoint/element Dudded 0 points plus parents, and 5 elements removed. finish geniee finish Merged 3 nodes. rmpoint/element Dudded 0 points plus parents, and 6 elements removed. finish geniee finish Merged 3 nodes. recon/0/0.10000000E+00/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 6 Recon2d: iteration number= 2 number of flips= 1 RECON: done. finish nodes in point set = 12 SGD: 0 total relaxations on point set of 12 indept. nodes. rmpoint/compress 42 points removed and 0 elements removed. The new point count is 31. The new element count is 29. finish nodes in point set = 12 recon/0/0.10000000E+00/ RECON2D reconnection loop2d: 1 RECON: done. finish nodes in point set = 12 SGD: 0 total relaxations on point set of 12 indept. nodes. rmpoint/compress 0 points removed and 0 elements removed. The new point count is 31. The new element count is 29. finish nodes in point set = 12 recon/0/0.10000000E+00/ RECON2D reconnection loop2d: 1 RECON: done. finish nodes in point set = 12 SGD: 0 total relaxations on point set of 12 indept. nodes. nodes in point set = 12 Face iteration, number of faces: 1 9 Face iteration, number of faces: 1 8 Face iteration, number of faces: 1 9 Face iteration, number of faces: 1 7 Face iteration, number of faces: 1 7 Face iteration, number of faces: 1 1 recon/0/0.10000000E+00/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 7 Recon2d: iteration number= 2 number of flips= 1 Recon2d: iteration number= 3 number of flips= 1 RECON: done. finish smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 72 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.2929688E+02 Max. area BEFORE smoothing= 0.2343750E+03 Iteration= 1; root mean square error= 0.1284957E+01 Iteration= 2; root mean square error= 0.4530739E+00 Iteration= 3; root mean square error= 0.1843337E+00 Iteration= 4; root mean square error= 0.1023544E+00 Iteration= 5; root mean square error= 0.6231207E-01 Iteration= 6; root mean square error= 0.3950751E-01 Iteration= 7; root mean square error= 0.2552321E-01 Iteration= 8; root mean square error= 0.1662822E-01 Iteration= 9; root mean square error= 0.1087580E-01 Min. area AFTER smoothing= 0.2929688E+02 Max. area AFTER smoothing= 0.1999616E+03 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 2 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 72 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.2929688E+02 Max. area BEFORE smoothing= 0.1999616E+03 Iteration= 1; root mean square error= 0.5605925E+00 Iteration= 2; root mean square error= 0.1623878E+00 Iteration= 3; root mean square error= 0.7837441E-01 Iteration= 4; root mean square error= 0.4543409E-01 Iteration= 5; root mean square error= 0.2877185E-01 Iteration= 6; root mean square error= 0.1861199E-01 Iteration= 7; root mean square error= 0.1210648E-01 Iteration= 8; root mean square error= 0.7888605E-02 Iteration= 9; root mean square error= 0.5143570E-02 Min. area AFTER smoothing= 0.2929688E+02 Max. area AFTER smoothing= 0.2006933E+03 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 1 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 72 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.2929688E+02 Max. area BEFORE smoothing= 0.2006933E+03 Iteration= 1; root mean square error= 0.4845180E+00 Iteration= 2; root mean square error= 0.1773996E+00 Iteration= 3; root mean square error= 0.6408765E-01 Iteration= 4; root mean square error= 0.3348785E-01 Iteration= 5; root mean square error= 0.2065250E-01 Iteration= 6; root mean square error= 0.1358303E-01 Iteration= 7; root mean square error= 0.9075100E-02 Iteration= 8; root mean square error= 0.6078992E-02 Iteration= 9; root mean square error= 0.4070452E-02 Min. area AFTER smoothing= 0.2929688E+02 Max. area AFTER smoothing= 0.2079616E+03 recon 0 RECON2D reconnection loop2d: 1 RECON: done. dump gmv tmp_tri05.gmv motri cmo/modatt/-def-/-def-/ioflag/x finish # Note the narrow tip with 25 edge length has neg ccoefs # that can not be resolved without adding nodes. # Fix this after joining with octree base mesh dump gmv top_tri_smooth.gmv motri cmo/modatt/-def-/-def-/ioflag/x finish dump avs top_tri_smooth.inp motri cmo/modatt/-def-/-def-/ioflag/x finish # this adds points where we do not want them recon 1 RECON2D reconnection loop2d: 1 Matbld0tri: total neg 28 mincoef= -2.3828125E+01 maxcoef= 4.0675799E+0 1 Negative coeff2d: sort edge vor_coeff row_node col_node Negative coeff2d: 1 99 -2.3828125E+01 18 19 Negative coeff2d: 2 102 -2.3828125E+01 19 18 Negative coeff2d: 3 97 -1.4322917E+01 18 17 Negative coeff2d: 4 93 -1.4322917E+01 17 18 Negative coeff2d: 5 163 -1.0825786E+01 30 36 Negative coeff2d: 6 197 -1.0825786E+01 36 30 Negative coeff2d: 7 75 -7.2206848E+00 14 15 Negative coeff2d: 8 79 -7.2206848E+00 15 14 Negative coeff2d: 9 91 -6.0937500E+00 17 16 Negative coeff2d: 10 87 -6.0937500E+00 16 17 Negative coeff2d: 11 73 -6.0836719E+00 14 13 Negative coeff2d: 12 69 -6.0836719E+00 13 14 Negative coeff2d: 13 29 -4.5674738E+00 6 7 Negative coeff2d: 14 32 -4.5674738E+00 7 6 Negative coeff2d: 15 67 -3.6780967E+00 13 12 Negative coeff2d: 16 63 -3.6780967E+00 12 13 Negative coeff2d: 17 191 -2.3413522E+00 35 24 Negative coeff2d: 18 129 -2.3413522E+00 24 35 Negative coeff2d: 19 61 -1.7705831E+00 12 11 Negative coeff2d: 20 58 -1.7705831E+00 11 12 Negative coeff2d: 21 85 -1.6741071E+00 16 15 Negative coeff2d: 22 81 -1.6741071E+00 15 16 Negative coeff2d: 23 22 -1.2249909E+00 5 6 Negative coeff2d: 24 27 -1.2249909E+00 6 5 Negative coeff2d: 25 15 -3.8225535E-01 4 3 Negative coeff2d: 26 13 -3.8225535E-01 3 4 Negative coeff2d: 27 52 -3.0984665E-01 10 11 Negative coeff2d: 28 56 -3.0984665E-01 11 10 Edge iteration refine_edge_2d: 1 14 Edge-refined tris: before= 96 after= 110 Matbld0tri: total neg 6 mincoef= -1.0381002E+00 maxcoef= 1.9165433E+0 1 Negative coeff2d: sort edge vor_coeff row_node col_node Negative coeff2d: 1 157 -1.0381002E+00 29 16 Negative coeff2d: 2 86 -1.0381002E+00 16 29 Negative coeff2d: 3 93 -9.2128722E-01 17 31 Negative coeff2d: 4 169 -9.2128722E-01 31 17 Negative coeff2d: 5 199 -1.1368684E-13 36 18 Negative coeff2d: 6 99 -1.1368684E-13 18 36 RECON2D reconnection loop2d: 2 Recon2d: iteration number= 1 number of flips= 2 Matbld0tri: total neg 2 mincoef= -1.1368684E-13 maxcoef= 1.9165433E+0 1 Negative coeff2d: sort edge vor_coeff row_node col_node Negative coeff2d: 1 195 -1.1368684E-13 36 18 Negative coeff2d: 2 97 -1.1368684E-13 18 36 RECON: done. cmo setatt motri itetclr 1 110 values reset for attribute itetclr cmo setatt motri imt 1 86 values reset for attribute imt1 dump gmv top_tri_recon.gmv motri cmo/modatt/-def-/-def-/ioflag/x finish dump avs top_tri_recon.inp motri cmo/modatt/-def-/-def-/ioflag/x finish finish LaGriT successfully completed