* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.3.3 Linux * * * Date Compile: 2023/12/20 * * * Run Time: 2023/Dec 20 13:23:46 * * * Manual: https://lagrit.lanl.gov * * * * * * * * * * * * * * * * * * * * * * * * * * * * * -----oOo----- LaGriT V3 LA-CC-15-069 https://github.com/lanl/LaGriT Copyright 2016. Triad National Security, LLC. All rights reserved. This program was produced under U.S. Government contract 89233218CNA000001 for Los Alamos National Laboratory (LANL), which is operated by Triad National Security, LLC for the U.S. Department of Energy/National Nuclear Security Administration. All rights in the program are reserved by Triad National Security, LLC, and the U.S. Department of Energy/National Nuclear Security Administration. The Government is granted for itself and others acting on its behalf a nonexclusive, paid-up, irrevocable worldwide license in this material to reproduce, prepare derivative works, distribute copies to the public, perform publicly and display publicly, and to permit others to do so. This software is open source and available under the BSD-3 License. -----oOo----- Output log file: outx3dgen Command log file: logx3dgen # LaGriT Test 2D smooth and massage # read the line polygon and triangulate # mesh object must be planar # uncomment dump/avs for further evaluations read avs input_polygon.inp moin geniee finish cmo/status/brief The current-mesh-object(CMO) is: moin 1 Mesh Object name: moin number of nodes = 7 number of elements = 7 dimensions geometry = 3 element type = lin dimensions topology = 1 2 nodes 2 faces 1 edges boundary flag = 16000000 status = active finish cmo/printatt/moin/-xyz-/minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 0.000000000E+00 2.000000000E+04 2.000000000E+04 7 yic 0.000000000E+00 1.000000000E+04 1.000000000E+04 7 zic 0.000000000E+00 0.000000000E+00 0.000000000E+00 7 cmo create cmotri///triplane copypts cmotri moin copypts added nodes 1 to 7 cmo/delete/moin Released Mesh Object: moin cmo select cmotri # RUN 00 triangulate # triangulate inside the polygon triangulate/clockwise TRIANGULATE: Number of triangle cells = 5 filter/1,0,0 FILTER:Use internal epsilonl value = 0.496506830649E-08 FILTER:Dudding duplicate points: 0 finish rmpoint/compress 0 points removed and 0 elements removed. RMPOINT: new point count is 7 RMPOINT: new element count is 5 finish resetpts/itp geniee finish finish cmo setatt cmotri imt 1 7 values reset for attribute imt1 cmo setatt cmotri itetclr 1 5 values reset for attribute itetclr resetpts/itp geniee finish # CHECK 00 coarse triangulation quality epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 0 elements with aspect ratio b/w .02 and .05: 1 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 0 elements with aspect ratio b/w .2 and .5 : 2 elements with aspect ratio b/w .5 and 1. : 2 min aspect ratio = 0.3703E-01 max aspect ratio = 0.8284E+00 epsilonvol: 1.1102230E-04 --------------------------------------- element volumes b/w 0.1000E+07 and 0.2187E+07: 0 element volumes b/w 0.2187E+07 and 0.4782E+07: 0 element volumes b/w 0.4782E+07 and 0.1046E+08: 0 element volumes b/w 0.1046E+08 and 0.2287E+08: 2 element volumes b/w 0.2287E+08 and 0.5000E+08: 3 min volume = 1.0000000E+06 max volume = 5.0000000E+07 ----------------------------------------------------------- 5 total elements evaluated. quality edge_max y cmo/addatt//edgemax/VDOUBLE/scalar/nelements finish epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- Edge max normalized by dividing with max value. element norm max edge length < .01: 0 element norm edge length b/w .01 and .02: 0 element norm edge length b/w .02 and .05: 0 element norm edge length b/w .05 and .1 : 0 element norm edge length b/w .1 and .2 : 0 element norm edge length b/w .2 and .5 : 2 element norm edge length b/w .5 and 1. : 3 min max edge length = 0.1000E+03 max max edge length = 0.2062E+05 ----------------------------------------------------------- 5 total elements evaluated. # dump avs out_tri00.inp cmotri # RUN 01 refine edge 1000 # break into triangles of nicer sizes define EDGELEN 1000. refine rivara///edge/1,0,0/EDGELEN///inclusive 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 4 Face iteration, number of faces: 1 5 Face iteration, number of faces: 1 6 Face iteration, number of faces: 1 8 Face iteration, number of faces: 1 12 Face iteration, number of faces: 1 16 Face iteration, number of faces: 1 20 Face iteration, number of faces: 1 28 Face iteration, number of faces: 1 34 Face iteration, number of faces: 1 41 Face iteration, number of faces: 1 56 Face iteration, number of faces: 1 56 Face iteration, number of faces: 1 62 Face iteration, number of faces: 1 52 Face iteration, number of faces: 1 31 Face iteration, number of faces: 1 12 Face iteration, number of faces: 1 2 recon/0/0.49650683E-08/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 26 Recon2d: iteration number= 2 number of flips= 5 RECON: done. finish # CHECK 01 refine edge 1000 quality epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 0 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 0 elements with aspect ratio b/w .2 and .5 : 64 elements with aspect ratio b/w .5 and 1. : 766 min aspect ratio = 0.2983E+00 max aspect ratio = 0.9889E+00 epsilonvol: 1.1102230E-04 --------------------------------------- element volumes b/w 0.2441E+05 and 0.4251E+05: 8 element volumes b/w 0.4251E+05 and 0.7401E+05: 60 element volumes b/w 0.7401E+05 and 0.1289E+06: 319 element volumes b/w 0.1289E+06 and 0.2244E+06: 417 element volumes b/w 0.2244E+06 and 0.3906E+06: 26 min volume = 2.4414062E+04 max volume = 3.9062500E+05 ----------------------------------------------------------- 830 total elements evaluated. quality edge_max y cmo/addatt//edgemax/VDOUBLE/scalar/nelements CMO_ADDATT warning: attribute already exist: cmotri edgemax finish epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- Edge max normalized by dividing with max value. element norm max edge length < .01: 0 element norm edge length b/w .01 and .02: 0 element norm edge length b/w .02 and .05: 0 element norm edge length b/w .05 and .1 : 0 element norm edge length b/w .1 and .2 : 0 element norm edge length b/w .2 and .5 : 6 element norm edge length b/w .5 and 1. : 824 min max edge length = 0.1000E+03 max max edge length = 0.9882E+03 ----------------------------------------------------------- 830 total elements evaluated. cmo printatt cmotri edgemax minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH edgemax 3.493856215E+02 9.882117688E+02 6.388261473E+02 830 # dump avs out_tri01a.inp cmotri # RUN 01b smooth # smooth to make nice distribution cmo select cmotri smooth/position/esug/1 0 0 SMOOTH: position option nodes in point set = 456 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.2441406E+05 Max. area BEFORE smoothing= 0.3906250E+06 Iteration= 1; root mean square error= 0.7503653E+02 Iteration= 2; root mean square error= 0.3945801E+02 Iteration= 3; root mean square error= 0.2431157E+02 Iteration= 4; root mean square error= 0.1850852E+02 Iteration= 5; root mean square error= 0.1567687E+02 Iteration= 6; root mean square error= 0.1383700E+02 Iteration= 7; root mean square error= 0.1246247E+02 Iteration= 8; root mean square error= 0.1136265E+02 Iteration= 9; root mean square error= 0.1044543E+02 Iteration= 10; root mean square error= 0.9659485E+01 Iteration= 11; root mean square error= 0.8973466E+01 Iteration= 12; root mean square error= 0.8366795E+01 Iteration= 13; root mean square error= 0.7825089E+01 Iteration= 14; root mean square error= 0.7337750E+01 Iteration= 15; root mean square error= 0.6896632E+01 Iteration= 16; root mean square error= 0.6495270E+01 Iteration= 17; root mean square error= 0.6128410E+01 Iteration= 18; root mean square error= 0.5791707E+01 Iteration= 19; root mean square error= 0.5481523E+01 Iteration= 20; root mean square error= 0.5194785E+01 Iteration= 21; root mean square error= 0.4928877E+01 Iteration= 22; root mean square error= 0.4681561E+01 Iteration= 23; root mean square error= 0.4450915E+01 Iteration= 24; root mean square error= 0.4235277E+01 Iteration= 25; root mean square error= 0.4033210E+01 Min. area AFTER smoothing= 0.8222425E+05 Max. area AFTER smoothing= 0.2584581E+06 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 74 Recon2d: iteration number= 2 number of flips= 2 RECON: done. smooth/position/esug/1 0 0 SMOOTH: position option nodes in point set = 456 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.3690220E+05 Max. area BEFORE smoothing= 0.2584581E+06 Iteration= 1; root mean square error= 0.5033363E+02 Iteration= 2; root mean square error= 0.2419419E+02 Iteration= 3; root mean square error= 0.1435078E+02 Iteration= 4; root mean square error= 0.1142554E+02 Iteration= 5; root mean square error= 0.9845981E+01 Iteration= 6; root mean square error= 0.8756779E+01 Iteration= 7; root mean square error= 0.7930349E+01 Iteration= 8; root mean square error= 0.7270073E+01 Iteration= 9; root mean square error= 0.6723729E+01 Iteration= 10; root mean square error= 0.6259833E+01 Iteration= 11; root mean square error= 0.5858123E+01 Iteration= 12; root mean square error= 0.5504926E+01 Iteration= 13; root mean square error= 0.5190643E+01 Iteration= 14; root mean square error= 0.4908298E+01 Iteration= 15; root mean square error= 0.4652664E+01 Iteration= 16; root mean square error= 0.4419707E+01 Iteration= 17; root mean square error= 0.4206247E+01 Iteration= 18; root mean square error= 0.4009718E+01 Iteration= 19; root mean square error= 0.3828015E+01 Iteration= 20; root mean square error= 0.3659389E+01 Iteration= 21; root mean square error= 0.3502366E+01 Iteration= 22; root mean square error= 0.3355694E+01 Iteration= 23; root mean square error= 0.3218302E+01 Iteration= 24; root mean square error= 0.3089266E+01 Iteration= 25; root mean square error= 0.2967783E+01 Min. area AFTER smoothing= 0.8224881E+05 Max. area AFTER smoothing= 0.2617622E+06 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 17 RECON: done. smooth/position/esug/1 0 0 SMOOTH: position option nodes in point set = 456 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.4372314E+05 Max. area BEFORE smoothing= 0.3373561E+06 Iteration= 1; root mean square error= 0.3401122E+02 Iteration= 2; root mean square error= 0.1482737E+02 Iteration= 3; root mean square error= 0.8052817E+01 Iteration= 4; root mean square error= 0.6112795E+01 Iteration= 5; root mean square error= 0.5230872E+01 Iteration= 6; root mean square error= 0.4713861E+01 Iteration= 7; root mean square error= 0.4355679E+01 Iteration= 8; root mean square error= 0.4078773E+01 Iteration= 9; root mean square error= 0.3849724E+01 Iteration= 10; root mean square error= 0.3652172E+01 Iteration= 11; root mean square error= 0.3477179E+01 Iteration= 12; root mean square error= 0.3319385E+01 Iteration= 13; root mean square error= 0.3175321E+01 Iteration= 14; root mean square error= 0.3042595E+01 Iteration= 15; root mean square error= 0.2919467E+01 Iteration= 16; root mean square error= 0.2804620E+01 Iteration= 17; root mean square error= 0.2697021E+01 Iteration= 18; root mean square error= 0.2595839E+01 Iteration= 19; root mean square error= 0.2500389E+01 Iteration= 20; root mean square error= 0.2410101E+01 Iteration= 21; root mean square error= 0.2324490E+01 Iteration= 22; root mean square error= 0.2243141E+01 Iteration= 23; root mean square error= 0.2165696E+01 Iteration= 24; root mean square error= 0.2091843E+01 Iteration= 25; root mean square error= 0.2021308E+01 Min. area AFTER smoothing= 0.7704365E+05 Max. area AFTER smoothing= 0.2599663E+06 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 11 RECON: done. filter/1,0,0 FILTER:Use internal epsilonl value = 0.496506830649E-08 FILTER:Dudding duplicate points: 0 rmpoint compress 0 points removed and 0 elements removed. RMPOINT: new point count is 456 RMPOINT: new element count is 830 resetpts/itp geniee finish # dump avs out_tri01b.inp cmotri # CHECK 01b smooth esug quality epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 0 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 0 elements with aspect ratio b/w .2 and .5 : 0 elements with aspect ratio b/w .5 and 1. : 830 min aspect ratio = 0.5793E+00 max aspect ratio = 0.1000E+01 epsilonvol: 1.1102230E-04 --------------------------------------- element volumes b/w 0.7704E+05 and 0.9826E+05: 72 element volumes b/w 0.9826E+05 and 0.1253E+06: 131 element volumes b/w 0.1253E+06 and 0.1598E+06: 278 element volumes b/w 0.1598E+06 and 0.2038E+06: 306 element volumes b/w 0.2038E+06 and 0.2600E+06: 43 min volume = 7.7043648E+04 max volume = 2.5996627E+05 ----------------------------------------------------------- 830 total elements evaluated. quality edge_max y cmo/addatt//edgemax/VDOUBLE/scalar/nelements CMO_ADDATT warning: attribute already exist: cmotri edgemax finish epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- Edge max normalized by dividing with max value. element norm max edge length < .01: 0 element norm edge length b/w .01 and .02: 0 element norm edge length b/w .02 and .05: 0 element norm edge length b/w .05 and .1 : 0 element norm edge length b/w .1 and .2 : 0 element norm edge length b/w .2 and .5 : 0 element norm edge length b/w .5 and 1. : 830 min max edge length = 0.1000E+03 max max edge length = 0.9199E+03 ----------------------------------------------------------- 830 total elements evaluated. cmo printatt cmotri edgemax minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH edgemax 4.621758514E+02 9.198557544E+02 4.576799030E+02 830 # RUN 02 massage 200 define MAXEDGE 300. define MINEDGE 100. define DAMAGE .01 massage/MAXEDGE MINEDGE DAMAGE/1,0,0/ nodes in point set = 456 recon/0/0.10000000E-01/ RECON2D reconnection loop2d: 1 RECON: done. finish nodes in point set = 456 SGD: 8 total relaxations on point set of 456 indept. nodes. rmpoint/compress 0 points removed and 0 elements removed. RMPOINT: new point count is 456 RMPOINT: new element count is 830 finish nodes in point set = 456 recon/0/0.10000000E-01/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 2 RECON: done. finish nodes in point set = 456 SGD: 3 total relaxations on point set of 456 indept. nodes. rmpoint/compress 0 points removed and 0 elements removed. RMPOINT: new point count is 456 RMPOINT: new element count is 830 finish nodes in point set = 456 recon/0/0.10000000E-01/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 1 RECON: done. finish nodes in point set = 456 SGD: 0 total relaxations on point set of 456 indept. nodes. nodes in point set = 456 Face iteration, number of faces: 1 165 Face iteration, number of faces: 1 207 Face iteration, number of faces: 1 276 Face iteration, number of faces: 1 352 Face iteration, number of faces: 1 452 Face iteration, number of faces: 1 582 Face iteration, number of faces: 1 664 Face iteration, number of faces: 1 652 Face iteration, number of faces: 1 466 Face iteration, number of faces: 1 215 Face iteration, number of faces: 1 49 Face iteration, number of faces: 1 7 recon/0/0.10000000E-01/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 898 Recon2d: iteration number= 2 number of flips= 4 RECON: done. finish rmpoint compress 0 points removed and 0 elements removed. RMPOINT: new point count is 4543 RMPOINT: new element count is 8764 # CHECK 02 massage 200 quality epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 0 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 0 elements with aspect ratio b/w .2 and .5 : 10 elements with aspect ratio b/w .5 and 1. : 8754 min aspect ratio = 0.4659E+00 max aspect ratio = 0.1000E+01 epsilonvol: 1.1102230E-04 --------------------------------------- element volumes b/w 0.5401E+04 and 0.7936E+04: 263 element volumes b/w 0.7936E+04 and 0.1166E+05: 3137 element volumes b/w 0.1166E+05 and 0.1713E+05: 2920 element volumes b/w 0.1713E+05 and 0.2518E+05: 2282 element volumes b/w 0.2518E+05 and 0.3699E+05: 162 min volume = 5.4006058E+03 max volume = 3.6993526E+04 ----------------------------------------------------------- 8764 total elements evaluated. quality edge_max y cmo/addatt//edgemax/VDOUBLE/scalar/nelements CMO_ADDATT warning: attribute already exist: cmotri edgemax finish epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- Edge max normalized by dividing with max value. element norm max edge length < .01: 0 element norm edge length b/w .01 and .02: 0 element norm edge length b/w .02 and .05: 0 element norm edge length b/w .05 and .1 : 0 element norm edge length b/w .1 and .2 : 0 element norm edge length b/w .2 and .5 : 12 element norm edge length b/w .5 and 1. : 8752 min max edge length = 0.1000E+03 max max edge length = 0.3065E+03 ----------------------------------------------------------- 8764 total elements evaluated. cmo printatt cmotri edgemax minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH edgemax 1.508860489E+02 3.065087170E+02 1.556226682E+02 8764 # dump avs out_tri02a.inp cmotri # RUN 02b smooth 200 smooth/position/esug/1 0 0 SMOOTH: position option nodes in point set = 4543 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.5400606E+04 Max. area BEFORE smoothing= 0.3699353E+05 Iteration= 1; root mean square error= 0.2356595E+02 Iteration= 2; root mean square error= 0.9872066E+01 Iteration= 3; root mean square error= 0.5655145E+01 Iteration= 4; root mean square error= 0.4274960E+01 Iteration= 5; root mean square error= 0.3522895E+01 Iteration= 6; root mean square error= 0.3027468E+01 Iteration= 7; root mean square error= 0.2670318E+01 Iteration= 8; root mean square error= 0.2398249E+01 Iteration= 9; root mean square error= 0.2182845E+01 Iteration= 10; root mean square error= 0.2007242E+01 Iteration= 11; root mean square error= 0.1860717E+01 Iteration= 12; root mean square error= 0.1736109E+01 Iteration= 13; root mean square error= 0.1628456E+01 Iteration= 14; root mean square error= 0.1534209E+01 Iteration= 15; root mean square error= 0.1450774E+01 Iteration= 16; root mean square error= 0.1376203E+01 Iteration= 17; root mean square error= 0.1309014E+01 Iteration= 18; root mean square error= 0.1248054E+01 Iteration= 19; root mean square error= 0.1192413E+01 Iteration= 20; root mean square error= 0.1141364E+01 Iteration= 21; root mean square error= 0.1094317E+01 Iteration= 22; root mean square error= 0.1050788E+01 Iteration= 23; root mean square error= 0.1010374E+01 Iteration= 24; root mean square error= 0.9727394E+00 Iteration= 25; root mean square error= 0.9375972E+00 Min. area AFTER smoothing= 0.6804840E+04 Max. area AFTER smoothing= 0.2787275E+05 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 209 RECON: done. smooth/position/esug/1 0 0 SMOOTH: position option nodes in point set = 4543 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.6804840E+04 Max. area BEFORE smoothing= 0.2787275E+05 Iteration= 1; root mean square error= 0.1172705E+02 Iteration= 2; root mean square error= 0.5022222E+01 Iteration= 3; root mean square error= 0.2356728E+01 Iteration= 4; root mean square error= 0.1661420E+01 Iteration= 5; root mean square error= 0.1342065E+01 Iteration= 6; root mean square error= 0.1152057E+01 Iteration= 7; root mean square error= 0.1023827E+01 Iteration= 8; root mean square error= 0.9303592E+00 Iteration= 9; root mean square error= 0.8584720E+00 Iteration= 10; root mean square error= 0.8009427E+00 Iteration= 11; root mean square error= 0.7534795E+00 Iteration= 12; root mean square error= 0.7133760E+00 Iteration= 13; root mean square error= 0.6788418E+00 Iteration= 14; root mean square error= 0.6486450E+00 Iteration= 15; root mean square error= 0.6219087E+00 Iteration= 16; root mean square error= 0.5979900E+00 Iteration= 17; root mean square error= 0.5764054E+00 Iteration= 18; root mean square error= 0.5567833E+00 Iteration= 19; root mean square error= 0.5388325E+00 Iteration= 20; root mean square error= 0.5223207E+00 Iteration= 21; root mean square error= 0.5070601E+00 Iteration= 22; root mean square error= 0.4928968E+00 Iteration= 23; root mean square error= 0.4797030E+00 Iteration= 24; root mean square error= 0.4673718E+00 Iteration= 25; root mean square error= 0.4558124E+00 Min. area AFTER smoothing= 0.6889443E+04 Max. area AFTER smoothing= 0.2595500E+05 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 75 RECON: done. smooth/position/esug/1 0 0 SMOOTH: position option nodes in point set = 4543 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.6889443E+04 Max. area BEFORE smoothing= 0.2595500E+05 Iteration= 1; root mean square error= 0.7254046E+01 Iteration= 2; root mean square error= 0.3178381E+01 Iteration= 3; root mean square error= 0.1560877E+01 Iteration= 4; root mean square error= 0.1085901E+01 Iteration= 5; root mean square error= 0.8632930E+00 Iteration= 6; root mean square error= 0.7303082E+00 Iteration= 7; root mean square error= 0.6407725E+00 Iteration= 8; root mean square error= 0.5761905E+00 Iteration= 9; root mean square error= 0.5273974E+00 Iteration= 10; root mean square error= 0.4892348E+00 Iteration= 11; root mean square error= 0.4585573E+00 Iteration= 12; root mean square error= 0.4333326E+00 Iteration= 13; root mean square error= 0.4121892E+00 Iteration= 14; root mean square error= 0.3941696E+00 Iteration= 15; root mean square error= 0.3785868E+00 Iteration= 16; root mean square error= 0.3649376E+00 Iteration= 17; root mean square error= 0.3528460E+00 Iteration= 18; root mean square error= 0.3420269E+00 Iteration= 19; root mean square error= 0.3322609E+00 Iteration= 20; root mean square error= 0.3233773E+00 Iteration= 21; root mean square error= 0.3152410E+00 Iteration= 22; root mean square error= 0.3077444E+00 Iteration= 23; root mean square error= 0.3008007E+00 Iteration= 24; root mean square error= 0.2943390E+00 Iteration= 25; root mean square error= 0.2883011E+00 Min. area AFTER smoothing= 0.6962014E+04 Max. area AFTER smoothing= 0.2569392E+05 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 29 RECON: done. # CHECK 02b smooth after massage quality epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 0 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 0 elements with aspect ratio b/w .2 and .5 : 0 elements with aspect ratio b/w .5 and 1. : 8764 min aspect ratio = 0.5671E+00 max aspect ratio = 0.1000E+01 epsilonvol: 1.1102230E-04 --------------------------------------- element volumes b/w 0.6962E+04 and 0.9040E+04: 72 element volumes b/w 0.9040E+04 and 0.1174E+05: 1693 element volumes b/w 0.1174E+05 and 0.1524E+05: 4252 element volumes b/w 0.1524E+05 and 0.1979E+05: 2263 element volumes b/w 0.1979E+05 and 0.2569E+05: 484 min volume = 6.9620138E+03 max volume = 2.5693923E+04 ----------------------------------------------------------- 8764 total elements evaluated. quality edge_max y cmo/addatt//edgemax/VDOUBLE/scalar/nelements CMO_ADDATT warning: attribute already exist: cmotri edgemax finish epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- Edge max normalized by dividing with max value. element norm max edge length < .01: 0 element norm edge length b/w .01 and .02: 0 element norm edge length b/w .02 and .05: 0 element norm edge length b/w .05 and .1 : 0 element norm edge length b/w .1 and .2 : 0 element norm edge length b/w .2 and .5 : 8 element norm edge length b/w .5 and 1. : 8756 min max edge length = 0.1000E+03 max max edge length = 0.3035E+03 ----------------------------------------------------------- 8764 total elements evaluated. cmo printatt cmotri edgemax minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH edgemax 1.452597882E+02 3.034671813E+02 1.582073931E+02 8764 # dump avs out_tri02b.inp cmotri rmpoint compress 0 points removed and 0 elements removed. RMPOINT: new point count is 4543 RMPOINT: new element count is 8764 resetpts itp geniee finish pset/pedge/attribute itp/1,0,0/ge 10 THE PSET pedge HAS 320 POINTS cmo create cmotmp copypts/cmotmp/cmotri/0 0/pset,get,pedge copypts added nodes 1 to 320 compute/distance_field/cmotri/cmotmp/dfield cmo/addatt/cmotri/temp_tmp_var/VDOUBLE/scalar/nnodes/linear/temporary//-1.0 finish cmo/addatt/cmotri/xobj_tmp_var/VDOUBLE/scalar/nnodes/linear/temporary//-1.0 finish cmo/addatt/cmotri/yobj_tmp_var/VDOUBLE/scalar/nnodes/linear/temporary//-1.0 finish cmo/addatt/cmotri/zobj_tmp_var/VDOUBLE/scalar/nnodes/linear/temporary//-1.0 finish cmo/addatt/cmotri/dfield/VDOUBLE/scalar/nnodes/linear/permanent//-1.0 finish interpolate/voronoi/cmotri/xobj_tmp_var/1 0 0/cmotmp/xic/keepatt nodes in indexed point set = 4543 cmo/addatt/cmotri/pt_gtg/VINT/scalar/nnodes//permanent/agfx/0 finish INTRP METHOD: voronoi FUNCTION: linear TIEBREAKER: maxtie FLAG: plus1 0.200010000000E+05 pt_gtg attribute will be added and kept. Build kdtree0 done - assign nearest node for each. Sink point Points Searched Points Found Percent Done 1137 1137 1137 25 % 2273 2273 2273 50 % 3409 3409 3409 75 % 4543 4543 4543 100% Total Source Nodes: 320 Total Sink Nodes: 4543 interpolate/voronoi done. cmo select cmotri finish cmo kdtree release finish finish interpolate/voronoi/cmotri/yobj_tmp_var/1 0 0/cmotmp/yic/keepatt nodes in indexed point set = 4543 INTRP METHOD: voronoi FUNCTION: linear TIEBREAKER: maxtie FLAG: plus1 0.100010000000E+05 SKIPPING POINT SEARCH... using lookup attribute pt_gtg 1137 1137 1137 25 % 2273 2273 2273 50 % 3409 3409 3409 75 % 4543 4543 4543 100% Total Source Nodes: 320 Total Sink Nodes: 4543 interpolate/voronoi done. cmo select cmotri finish cmo kdtree release finish finish interpolate/voronoi/cmotri/zobj_tmp_var/1 0 0/cmotmp/zic nodes in indexed point set = 4543 INTRP METHOD: voronoi FUNCTION: linear TIEBREAKER: maxtie FLAG: plus1 0.100000000000E+01 SKIPPING POINT SEARCH... using lookup attribute pt_gtg 1137 1137 1137 25 % 2273 2273 2273 50 % 3409 3409 3409 75 % 4543 4543 4543 100% Total Source Nodes: 320 Total Sink Nodes: 4543 interpolate/voronoi done. cmo DELATT cmotri/pt_gtg finish cmo select cmotri finish cmo kdtree release finish finish math/sub/cmotri/temp_tmp_var/1,0,0/cmotri/xic/cmotri/xobj_tmp_var finish math/power/cmotri/dfield/1,0,0/cmotri/temp_tmp_var/2.0 finish math/sub/cmotri/temp_tmp_var/1,0,0/cmotri/yic/cmotri/yobj_tmp_var finish math/power/cmotri/temp_tmp_var/1,0,0/cmotri/temp_tmp_var/2.0 finish math/add/cmotri/dfield/1 0 0/cmotri/dfield/cmotri/temp_tmp_var finish math/sub/cmotri/temp_tmp_var/1,0,0/cmotri/zic/cmotri/zobj_tmp_var finish math/power/cmotri/temp_tmp_var/1,0,0/cmotri/temp_tmp_var/2.0 finish math/add/cmotri/dfield/1 0 0/cmotri/dfield/cmotri/temp_tmp_var finish math/power/cmotri/dfield/1,0,0/cmotri/dfield/0.5 finish cmo/delatt/cmotri/temp_tmp_var finish cmo/delatt/cmotri/xobj_tmp_var finish cmo/delatt/cmotri/yobj_tmp_var finish cmo/delatt/cmotri/zobj_tmp_var finish COMPUTE distance_field: Done. cmo/delete/cmotmp Released Mesh Object: cmotmp cmo select cmotri # RUN 03 massage de-refine 600 # select nodes inside mesh to massage define MAXEDGE 1200. define MINEDGE 600. define DAMAGE .01 pset/pmass/attribute dfield/1,0,0/gt 800. THE PSET pmass HAS 2956 POINTS massage/MAXEDGE MINEDGE DAMAGE/pset,get,pmass/nosmooth nodes in point set = 2956 rmpoint/element Dudded 0 points plus parents, and 1474 elements removed. finish geniee finish Merged 737 nodes. rmpoint/element Dudded 0 points plus parents, and 1154 elements removed. finish geniee finish Merged 577 nodes. rmpoint/element Dudded 0 points plus parents, and 886 elements removed. finish geniee finish Merged 443 nodes. rmpoint/element Dudded 0 points plus parents, and 640 elements removed. finish geniee finish Merged 320 nodes. rmpoint/element Dudded 0 points plus parents, and 478 elements removed. finish geniee finish Merged 239 nodes. rmpoint/element Dudded 0 points plus parents, and 356 elements removed. finish geniee finish Merged 178 nodes. rmpoint/element Dudded 0 points plus parents, and 252 elements removed. finish geniee finish Merged 126 nodes. rmpoint/element Dudded 0 points plus parents, and 196 elements removed. finish geniee finish Merged 98 nodes. rmpoint/element Dudded 0 points plus parents, and 136 elements removed. finish geniee finish Merged 68 nodes. rmpoint/element Dudded 0 points plus parents, and 94 elements removed. finish geniee finish Merged 47 nodes. rmpoint/element Dudded 0 points plus parents, and 62 elements removed. finish geniee finish Merged 31 nodes. rmpoint/element Dudded 0 points plus parents, and 32 elements removed. finish geniee finish Merged 16 nodes. rmpoint/element Dudded 0 points plus parents, and 10 elements removed. finish geniee finish Merged 5 nodes. recon/0/0.10000000E-01/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 92 Recon2d: iteration number= 2 number of flips= 34 Recon2d: iteration number= 3 number of flips= 22 Recon2d: iteration number= 4 number of flips= 10 Recon2d: iteration number= 5 number of flips= 1 Recon2d: iteration number= 6 number of flips= 2 Recon2d: iteration number= 7 number of flips= 3 RECON: done. finish rmpoint/compress 2885 points removed and 0 elements removed. RMPOINT: new point count is 1658 RMPOINT: new element count is 2994 finish nodes in point set = 71 rmpoint/element Dudded 0 points plus parents, and 14 elements removed. finish geniee finish Merged 7 nodes. recon/0/0.10000000E-01/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 17 Recon2d: iteration number= 2 number of flips= 9 Recon2d: iteration number= 3 number of flips= 10 Recon2d: iteration number= 4 number of flips= 7 Recon2d: iteration number= 5 number of flips= 2 Recon2d: iteration number= 6 number of flips= 3 Recon2d: iteration number= 7 number of flips= 1 Recon2d: iteration number= 8 number of flips= 1 Recon2d: iteration number= 9 number of flips= 1 Recon2d: iteration number= 10 number of flips= 1 RECON: done. finish rmpoint/compress 7 points removed and 0 elements removed. RMPOINT: new point count is 1651 RMPOINT: new element count is 2980 finish nodes in point set = 64 rmpoint/element Dudded 0 points plus parents, and 2 elements removed. finish geniee finish Merged 1 nodes. recon/0/0.10000000E-01/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 4 Recon2d: iteration number= 2 number of flips= 2 Recon2d: iteration number= 3 number of flips= 1 Recon2d: iteration number= 4 number of flips= 1 Recon2d: iteration number= 5 number of flips= 1 Recon2d: iteration number= 6 number of flips= 2 Recon2d: iteration number= 7 number of flips= 2 Recon2d: iteration number= 8 number of flips= 1 Recon2d: iteration number= 9 number of flips= 1 RECON: done. finish rmpoint/compress 1 points removed and 0 elements removed. RMPOINT: new point count is 1650 RMPOINT: new element count is 2978 finish nodes in point set = 63 recon/0/0.10000000E-01/ RECON2D reconnection loop2d: 1 RECON: done. finish rmpoint/compress 0 points removed and 0 elements removed. RMPOINT: new point count is 1650 RMPOINT: new element count is 2978 finish MASSAGE: hit max. iteration limit. nodes in point set = 63 Face iteration, number of faces: 1 28 Face iteration, number of faces: 1 32 Face iteration, number of faces: 1 32 Face iteration, number of faces: 1 23 Face iteration, number of faces: 1 10 Face iteration, number of faces: 1 2 Face iteration, number of faces: 1 2 Face iteration, number of faces: 1 1 recon/0/0.10000000E-01/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 77 Recon2d: iteration number= 2 number of flips= 34 Recon2d: iteration number= 3 number of flips= 9 Recon2d: iteration number= 4 number of flips= 5 Recon2d: iteration number= 5 number of flips= 5 Recon2d: iteration number= 6 number of flips= 3 RECON: done. finish rmpoint compress 0 points removed and 0 elements removed. RMPOINT: new point count is 1780 RMPOINT: new element count is 3238 # CHECK 03 massage de-refine quality epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 0 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 28 elements with aspect ratio b/w .2 and .5 : 205 elements with aspect ratio b/w .5 and 1. : 3005 min aspect ratio = 0.1037E+00 max aspect ratio = 0.9999E+00 epsilonvol: 1.1102230E-04 --------------------------------------- element volumes b/w 0.5630E+04 and 0.1356E+05: 1229 element volumes b/w 0.1356E+05 and 0.3264E+05: 1478 element volumes b/w 0.3264E+05 and 0.7858E+05: 131 element volumes b/w 0.7858E+05 and 0.1892E+06: 201 element volumes b/w 0.1892E+06 and 0.4555E+06: 199 min volume = 5.6299363E+03 max volume = 4.5551568E+05 ----------------------------------------------------------- 3238 total elements evaluated. quality edge_max y cmo/addatt//edgemax/VDOUBLE/scalar/nelements CMO_ADDATT warning: attribute already exist: cmotri edgemax finish epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- Edge max normalized by dividing with max value. element norm max edge length < .01: 0 element norm edge length b/w .01 and .02: 0 element norm edge length b/w .02 and .05: 0 element norm edge length b/w .05 and .1 : 0 element norm edge length b/w .1 and .2 : 2106 element norm edge length b/w .2 and .5 : 591 element norm edge length b/w .5 and 1. : 541 min max edge length = 0.1000E+03 max max edge length = 0.1181E+04 ----------------------------------------------------------- 3238 total elements evaluated. cmo printatt cmotri edgemax minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH edgemax 1.511033628E+02 1.181202998E+03 1.030099635E+03 3238 # dump avs out_tri03a.inp cmotri # RUN 03b smooth after de-refine cmo select cmotri smooth/position/esug/1 0 0 SMOOTH: position option nodes in point set = 1780 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.5629936E+04 Max. area BEFORE smoothing= 0.4555157E+06 Iteration= 1; root mean square error= 0.6516577E+02 Iteration= 2; root mean square error= 0.5582389E+02 Iteration= 3; root mean square error= 0.5340996E+02 Iteration= 4; root mean square error= 0.4901882E+02 Iteration= 5; root mean square error= 0.4371524E+02 Iteration= 6; root mean square error= 0.3885613E+02 Iteration= 7; root mean square error= 0.3443594E+02 Iteration= 8; root mean square error= 0.3053791E+02 Iteration= 9; root mean square error= 0.2707585E+02 Iteration= 10; root mean square error= 0.2401504E+02 Iteration= 11; root mean square error= 0.2132917E+02 Iteration= 12; root mean square error= 0.1897595E+02 Iteration= 13; root mean square error= 0.1691286E+02 Iteration= 14; root mean square error= 0.1510337E+02 Iteration= 15; root mean square error= 0.1351571E+02 Iteration= 16; root mean square error= 0.1212208E+02 Iteration= 17; root mean square error= 0.1089818E+02 Iteration= 18; root mean square error= 0.9822759E+01 Iteration= 19; root mean square error= 0.8877222E+01 Iteration= 20; root mean square error= 0.8045322E+01 Iteration= 21; root mean square error= 0.7312842E+01 Iteration= 22; root mean square error= 0.6667339E+01 Iteration= 23; root mean square error= 0.6097918E+01 Iteration= 24; root mean square error= 0.5595031E+01 Iteration= 25; root mean square error= 0.5150308E+01 Min. area AFTER smoothing= 0.1156955E+05 Max. area AFTER smoothing= 0.2076847E+06 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 658 Recon2d: iteration number= 2 number of flips= 154 Recon2d: iteration number= 3 number of flips= 50 Recon2d: iteration number= 4 number of flips= 5 Recon2d: iteration number= 5 number of flips= 3 Recon2d: iteration number= 6 number of flips= 3 RECON: done. smooth/position/esug/1 0 0 SMOOTH: position option nodes in point set = 1780 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.3726787E+04 Max. area BEFORE smoothing= 0.2739801E+06 Iteration= 1; root mean square error= 0.5048696E+02 Iteration= 2; root mean square error= 0.2819378E+02 Iteration= 3; root mean square error= 0.1993504E+02 Iteration= 4; root mean square error= 0.1623195E+02 Iteration= 5; root mean square error= 0.1403982E+02 Iteration= 6; root mean square error= 0.1250124E+02 Iteration= 7; root mean square error= 0.1130486E+02 Iteration= 8; root mean square error= 0.1031458E+02 Iteration= 9; root mean square error= 0.9464934E+01 Iteration= 10; root mean square error= 0.8720411E+01 Iteration= 11; root mean square error= 0.8059698E+01 Iteration= 12; root mean square error= 0.7468648E+01 Iteration= 13; root mean square error= 0.6937085E+01 Iteration= 14; root mean square error= 0.6457208E+01 Iteration= 15; root mean square error= 0.6022744E+01 Iteration= 16; root mean square error= 0.5628476E+01 Iteration= 17; root mean square error= 0.5269968E+01 Iteration= 18; root mean square error= 0.4943381E+01 Iteration= 19; root mean square error= 0.4645366E+01 Iteration= 20; root mean square error= 0.4372972E+01 Iteration= 21; root mean square error= 0.4123590E+01 Iteration= 22; root mean square error= 0.3894904E+01 Iteration= 23; root mean square error= 0.3684853E+01 Iteration= 24; root mean square error= 0.3491599E+01 Iteration= 25; root mean square error= 0.3313501E+01 Min. area AFTER smoothing= 0.1196117E+05 Max. area AFTER smoothing= 0.1116543E+06 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 234 Recon2d: iteration number= 2 number of flips= 10 Recon2d: iteration number= 3 number of flips= 1 RECON: done. # CHECK 03b smooth after de-refine quality epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 0 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 0 elements with aspect ratio b/w .2 and .5 : 16 elements with aspect ratio b/w .5 and 1. : 3222 min aspect ratio = 0.3107E+00 max aspect ratio = 0.9999E+00 epsilonvol: 1.1102230E-04 --------------------------------------- element volumes b/w 0.8688E+04 and 0.1465E+05: 41 element volumes b/w 0.1465E+05 and 0.2471E+05: 551 element volumes b/w 0.2471E+05 and 0.4168E+05: 1459 element volumes b/w 0.4168E+05 and 0.7029E+05: 1043 element volumes b/w 0.7029E+05 and 0.1186E+06: 144 min volume = 8.6876689E+03 max volume = 1.1855536E+05 ----------------------------------------------------------- 3238 total elements evaluated. quality edge_max y cmo/addatt//edgemax/VDOUBLE/scalar/nelements CMO_ADDATT warning: attribute already exist: cmotri edgemax finish epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- Edge max normalized by dividing with max value. element norm max edge length < .01: 0 element norm edge length b/w .01 and .02: 0 element norm edge length b/w .02 and .05: 0 element norm edge length b/w .05 and .1 : 0 element norm edge length b/w .1 and .2 : 0 element norm edge length b/w .2 and .5 : 851 element norm edge length b/w .5 and 1. : 2387 min max edge length = 0.1000E+03 max max edge length = 0.6327E+03 ----------------------------------------------------------- 3238 total elements evaluated. quality aspect y cmo/addatt//aratio/VDOUBLE/scalar/nelements finish epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 0 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 0 elements with aspect ratio b/w .2 and .5 : 16 elements with aspect ratio b/w .5 and 1. : 3222 min aspect ratio = 0.3107E+00 max aspect ratio = 0.9999E+00 ----------------------------------------------------------- 3238 total elements evaluated. cmo/addatt/cmotri/ang_mind/ang_mind ADDATT/ang_mind: creating new attribute: ang_mind cmo printatt cmotri edgemax minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH edgemax 1.768115362E+02 6.327184833E+02 4.559069471E+02 3238 cmo printatt cmotri aratio minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH aratio 3.107099577E-01 9.999339295E-01 6.892239718E-01 3238 cmo printatt cmotri ang_mind minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH ang_mind 1.852724051E+01 5.962343086E+01 4.109619035E+01 3238 # dump avs out_tri03b.inp cmotri # RUN 04 massage merge small edges define MAXEDGE 1.e+20 define MINEDGE .1 define DAMAGE .1 massage/MAXEDGE MINEDGE DAMAGE/1,0,0/ nodes in point set = 1780 recon/0/0.10000000E+00/ RECON2D reconnection loop2d: 1 RECON: done. finish nodes in point set = 1780 SGD: 137 total relaxations on point set of 1780 indept. nodes. rmpoint/compress 0 points removed and 0 elements removed. RMPOINT: new point count is 1780 RMPOINT: new element count is 3238 finish nodes in point set = 1780 recon/0/0.10000000E+00/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 39 Recon2d: iteration number= 2 number of flips= 1 RECON: done. finish nodes in point set = 1780 SGD: 23 total relaxations on point set of 1780 indept. nodes. rmpoint/compress 0 points removed and 0 elements removed. RMPOINT: new point count is 1780 RMPOINT: new element count is 3238 finish nodes in point set = 1780 recon/0/0.10000000E+00/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 6 RECON: done. finish nodes in point set = 1780 SGD: 0 total relaxations on point set of 1780 indept. nodes. nodes in point set = 1780 rmpoint compress 0 points removed and 0 elements removed. RMPOINT: new point count is 1780 RMPOINT: new element count is 3238 resetpts itp geniee finish # CHECK 04 massage merge small edges quality epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 0 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 0 elements with aspect ratio b/w .2 and .5 : 0 elements with aspect ratio b/w .5 and 1. : 3238 min aspect ratio = 0.5165E+00 max aspect ratio = 0.9999E+00 epsilonvol: 1.1102230E-04 --------------------------------------- element volumes b/w 0.1196E+05 and 0.1879E+05: 146 element volumes b/w 0.1879E+05 and 0.2952E+05: 843 element volumes b/w 0.2952E+05 and 0.4638E+05: 1451 element volumes b/w 0.4638E+05 and 0.7287E+05: 674 element volumes b/w 0.7287E+05 and 0.1145E+06: 124 min volume = 1.1961165E+04 max volume = 1.1448430E+05 ----------------------------------------------------------- 3238 total elements evaluated. quality edge_max y cmo/addatt//edgemax/VDOUBLE/scalar/nelements CMO_ADDATT warning: attribute already exist: cmotri edgemax finish epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- Edge max normalized by dividing with max value. element norm max edge length < .01: 0 element norm edge length b/w .01 and .02: 0 element norm edge length b/w .02 and .05: 0 element norm edge length b/w .05 and .1 : 0 element norm edge length b/w .1 and .2 : 0 element norm edge length b/w .2 and .5 : 902 element norm edge length b/w .5 and 1. : 2336 min max edge length = 0.1000E+03 max max edge length = 0.6327E+03 ----------------------------------------------------------- 3238 total elements evaluated. quality aspect y cmo/addatt//aratio/VDOUBLE/scalar/nelements CMO_ADDATT warning: attribute already exist: cmotri aratio finish epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 0 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 0 elements with aspect ratio b/w .2 and .5 : 0 elements with aspect ratio b/w .5 and 1. : 3238 min aspect ratio = 0.5165E+00 max aspect ratio = 0.9999E+00 ----------------------------------------------------------- 3238 total elements evaluated. cmo/addatt/cmotri/ang_mind/ang_mind CMO_ADDATT warning: attribute already exist: cmotri ang_mind ADDATT/ang_mind: writing to existing attribute: ang_mind cmo printatt cmotri edgemax minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH edgemax 2.063933970E+02 6.327184833E+02 4.263250862E+02 3238 cmo printatt cmotri aratio minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH aratio 5.164627823E-01 9.999216272E-01 4.834588449E-01 3238 cmo printatt cmotri ang_mind minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH ang_mind 1.852724051E+01 5.946247654E+01 4.093523603E+01 3238 # dump avs out_tri04a.inp cmotri # RUN 04b smooth after massage merge cmo select cmotri smooth/position/esug/1 0 0 SMOOTH: position option nodes in point set = 1780 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.1196117E+05 Max. area BEFORE smoothing= 0.1144843E+06 Iteration= 1; root mean square error= 0.2809953E+02 Iteration= 2; root mean square error= 0.1229607E+02 Iteration= 3; root mean square error= 0.7458821E+01 Iteration= 4; root mean square error= 0.5881224E+01 Iteration= 5; root mean square error= 0.5087326E+01 Iteration= 6; root mean square error= 0.4584816E+01 Iteration= 7; root mean square error= 0.4219348E+01 Iteration= 8; root mean square error= 0.3930712E+01 Iteration= 9; root mean square error= 0.3690999E+01 Iteration= 10; root mean square error= 0.3485416E+01 Iteration= 11; root mean square error= 0.3305231E+01 Iteration= 12; root mean square error= 0.3144866E+01 Iteration= 13; root mean square error= 0.3000523E+01 Iteration= 14; root mean square error= 0.2869484E+01 Iteration= 15; root mean square error= 0.2749720E+01 Iteration= 16; root mean square error= 0.2639661E+01 Iteration= 17; root mean square error= 0.2538057E+01 Iteration= 18; root mean square error= 0.2443890E+01 Iteration= 19; root mean square error= 0.2356313E+01 Iteration= 20; root mean square error= 0.2274611E+01 Iteration= 21; root mean square error= 0.2198174E+01 Iteration= 22; root mean square error= 0.2126474E+01 Iteration= 23; root mean square error= 0.2059054E+01 Iteration= 24; root mean square error= 0.1995512E+01 Iteration= 25; root mean square error= 0.1935497E+01 Min. area AFTER smoothing= 0.1220703E+05 Max. area AFTER smoothing= 0.8651605E+05 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 85 RECON: done. smooth/position/esug/1 0 0 SMOOTH: position option nodes in point set = 1780 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.1220703E+05 Max. area BEFORE smoothing= 0.8651605E+05 Iteration= 1; root mean square error= 0.1700973E+02 Iteration= 2; root mean square error= 0.7831059E+01 Iteration= 3; root mean square error= 0.4478101E+01 Iteration= 4; root mean square error= 0.3409824E+01 Iteration= 5; root mean square error= 0.2916082E+01 Iteration= 6; root mean square error= 0.2625400E+01 Iteration= 7; root mean square error= 0.2425743E+01 Iteration= 8; root mean square error= 0.2274483E+01 Iteration= 9; root mean square error= 0.2152415E+01 Iteration= 10; root mean square error= 0.2049697E+01 Iteration= 11; root mean square error= 0.1960763E+01 Iteration= 12; root mean square error= 0.1882190E+01 Iteration= 13; root mean square error= 0.1811719E+01 Iteration= 14; root mean square error= 0.1747775E+01 Iteration= 15; root mean square error= 0.1689209E+01 Iteration= 16; root mean square error= 0.1635158E+01 Iteration= 17; root mean square error= 0.1584953E+01 Iteration= 18; root mean square error= 0.1538067E+01 Iteration= 19; root mean square error= 0.1494076E+01 Iteration= 20; root mean square error= 0.1452634E+01 Iteration= 21; root mean square error= 0.1413455E+01 Iteration= 22; root mean square error= 0.1376301E+01 Iteration= 23; root mean square error= 0.1340970E+01 Iteration= 24; root mean square error= 0.1307291E+01 Iteration= 25; root mean square error= 0.1275117E+01 Min. area AFTER smoothing= 0.1220703E+05 Max. area AFTER smoothing= 0.8178562E+05 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 46 Recon2d: iteration number= 2 number of flips= 1 RECON: done. filter/1,0,0 FILTER:Use internal epsilonl value = 0.496506830649E-08 FILTER:Dudding duplicate points: 0 rmpoint compress 0 points removed and 0 elements removed. RMPOINT: new point count is 1780 RMPOINT: new element count is 3238 resetpts/itp geniee finish # CHECK 04b smooth after massage merge quality epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 0 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 0 elements with aspect ratio b/w .2 and .5 : 2 elements with aspect ratio b/w .5 and 1. : 3236 min aspect ratio = 0.4309E+00 max aspect ratio = 0.1000E+01 epsilonvol: 1.1102230E-04 --------------------------------------- element volumes b/w 0.1221E+05 and 0.1786E+05: 47 element volumes b/w 0.1786E+05 and 0.2612E+05: 448 element volumes b/w 0.2612E+05 and 0.3822E+05: 1115 element volumes b/w 0.3822E+05 and 0.5591E+05: 1366 element volumes b/w 0.5591E+05 and 0.8179E+05: 262 min volume = 1.2207031E+04 max volume = 8.1785622E+04 ----------------------------------------------------------- 3238 total elements evaluated. quality edge_max y cmo/addatt//edgemax/VDOUBLE/scalar/nelements CMO_ADDATT warning: attribute already exist: cmotri edgemax finish epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- Edge max normalized by dividing with max value. element norm max edge length < .01: 0 element norm edge length b/w .01 and .02: 0 element norm edge length b/w .02 and .05: 0 element norm edge length b/w .05 and .1 : 0 element norm edge length b/w .1 and .2 : 0 element norm edge length b/w .2 and .5 : 136 element norm edge length b/w .5 and 1. : 3102 min max edge length = 0.1000E+03 max max edge length = 0.5192E+03 ----------------------------------------------------------- 3238 total elements evaluated. quality aspect y cmo/addatt//aratio/VDOUBLE/scalar/nelements CMO_ADDATT warning: attribute already exist: cmotri aratio finish epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 0 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 0 elements with aspect ratio b/w .2 and .5 : 2 elements with aspect ratio b/w .5 and 1. : 3236 min aspect ratio = 0.4309E+00 max aspect ratio = 0.1000E+01 ----------------------------------------------------------- 3238 total elements evaluated. cmo/addatt/cmotri/ang_mind/ang_mind CMO_ADDATT warning: attribute already exist: cmotri ang_mind ADDATT/ang_mind: writing to existing attribute: ang_mind cmo printatt cmotri edgemax minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH edgemax 2.207399709E+02 5.192331048E+02 2.984931339E+02 3238 cmo printatt cmotri aratio minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH aratio 4.309068587E-01 9.999916634E-01 5.690848047E-01 3238 cmo printatt cmotri ang_mind minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH ang_mind 2.057960050E+01 5.984399222E+01 3.926439172E+01 3238 cmo printatt cmotri edgemax minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH edgemax 2.207399709E+02 5.192331048E+02 2.984931339E+02 3238 # dump avs out_tri04b.inp cmotri # RUN FINAL reconnect delaunay recon/1/ RECON2D reconnection loop2d: 1 RECON: done. # CHECK FINAL triangulation quality epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 0 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 0 elements with aspect ratio b/w .2 and .5 : 2 elements with aspect ratio b/w .5 and 1. : 3236 min aspect ratio = 0.4309E+00 max aspect ratio = 0.1000E+01 epsilonvol: 1.1102230E-04 --------------------------------------- element volumes b/w 0.1221E+05 and 0.1786E+05: 47 element volumes b/w 0.1786E+05 and 0.2612E+05: 448 element volumes b/w 0.2612E+05 and 0.3822E+05: 1115 element volumes b/w 0.3822E+05 and 0.5591E+05: 1366 element volumes b/w 0.5591E+05 and 0.8179E+05: 262 min volume = 1.2207031E+04 max volume = 8.1785622E+04 ----------------------------------------------------------- 3238 total elements evaluated. quality edge_max y cmo/addatt//edgemax/VDOUBLE/scalar/nelements CMO_ADDATT warning: attribute already exist: cmotri edgemax finish epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- Edge max normalized by dividing with max value. element norm max edge length < .01: 0 element norm edge length b/w .01 and .02: 0 element norm edge length b/w .02 and .05: 0 element norm edge length b/w .05 and .1 : 0 element norm edge length b/w .1 and .2 : 0 element norm edge length b/w .2 and .5 : 136 element norm edge length b/w .5 and 1. : 3102 min max edge length = 0.1000E+03 max max edge length = 0.5192E+03 ----------------------------------------------------------- 3238 total elements evaluated. quality aspect y cmo/addatt//aratio/VDOUBLE/scalar/nelements CMO_ADDATT warning: attribute already exist: cmotri aratio finish epsilonl, epsilonaspect: 4.9650683E-09 1.2239838E-25 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 0 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 0 elements with aspect ratio b/w .2 and .5 : 2 elements with aspect ratio b/w .5 and 1. : 3236 min aspect ratio = 0.4309E+00 max aspect ratio = 0.1000E+01 ----------------------------------------------------------- 3238 total elements evaluated. cmo/addatt/cmotri/ang_mind/ang_mind CMO_ADDATT warning: attribute already exist: cmotri ang_mind ADDATT/ang_mind: writing to existing attribute: ang_mind cmo printatt cmotri edgemax minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH edgemax 2.207399709E+02 5.192331048E+02 2.984931339E+02 3238 cmo printatt cmotri aratio minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH aratio 4.309068587E-01 9.999916634E-01 5.690848047E-01 3238 cmo printatt cmotri ang_mind minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH ang_mind 2.057960050E+01 5.984399222E+01 3.926439172E+01 3238 dump avs out_final.inp cmotri cmo/modatt/-def-/-def-/ioflag/x finish cmo printatt cmotri -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 1780 scalar 1 1 0 1 vector 3 3 0 1 nnodes 1780 1780 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 3238 3238 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 2 2 0 1 ndimensions_geom 2 2 0 1 nodes_per_element 3 3 0 1 edges_per_element 3 3 0 1 faces_per_element 3 3 0 1 isetwd 0 2 2 1780 ialias 0 4543 4543 1780 imt1 1 1 0 1780 itp1 0 10 10 1780 icr1 0 0 0 1780 isn1 0 0 0 1780 xic 0.000000000E+00 2.000000000E+04 2.000000000E+04 1780 yic 0.000000000E+00 1.000000000E+04 1.000000000E+04 1780 zic 0.000000000E+00 0.000000000E+00 0.000000000E+00 1780 xtetwd 0 0 0 3238 itetclr 1 1 0 3238 itettyp 3 3 0 3238 itetoff 0 9711 9711 3238 jtetoff 0 9711 9711 3238 itet 1 1780 1779 3238x3 jtet 1 16000000 15999999 3238x3 epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 4.965068306E-09 4.965068306E-09 0.000000000E+00 1 epsilona 1.110223025E-04 1.110223025E-04 0.000000000E+00 1 epsilonv 1.110223025E-04 1.110223025E-04 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 1780 1780 0 1 idebug 0 0 0 1 itypconv_sm 1 1 0 1 maxiter_sm 25 25 0 1 tolconv_sm 1.000000000E+00 1.000000000E+00 0.000000000E+00 1 nnfreq 1 1 0 1 ivoronoi 1 1 0 1 iopt2to2 2 2 0 1 xmin 0.000000000E+00 0.000000000E+00 0.000000000E+00 1 ymin 0.000000000E+00 0.000000000E+00 0.000000000E+00 1 zmin 0.000000000E+00 0.000000000E+00 0.000000000E+00 1 xmax 2.000000000E+04 2.000000000E+04 0.000000000E+00 1 ymax 1.000000000E+04 1.000000000E+04 0.000000000E+00 1 zmax 0.000000000E+00 0.000000000E+00 0.000000000E+00 1 kdtree_level 0 0 0 1 max_number_sets 64 64 0 1 number_of_psets 2 2 0 1 number_of_eltsets 0 0 0 1 number_of_fsets 0 0 0 1 edgemax 2.207399709E+02 5.192331048E+02 2.984931339E+02 3238 dfield 0.000000000E+00 4.007218110E+03 4.007218110E+03 1780 aratio 4.309068587E-01 9.999916634E-01 5.690848047E-01 3238 ang_mind 2.057960050E+01 5.984399222E+01 3.926439172E+01 3238 finish LaGriT successfully completed