* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Jun 19 16:32:26 * * * 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 # join top and bottom pieces and write final files # Can not use prism top, it was formed with stack and # the axis switched to make turn y to Z upward # the points can be used for connect # This is the merged 25m mesh for interpolate # it will be used to find elements connected # across the slope that we need to remove read avs tet_25m.inp mo_merge cmo/addatt/mo_merge/top/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/mo_merge/bottom/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/mo_merge/left_w/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/mo_merge/right_e/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/mo_merge/back_n/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/mo_merge/front_s/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish geniee finish cmo/status/brief The current-mesh-object(CMO) is: mo_merge 1 Mesh Object name: mo_merge number of nodes = 634148 number of elements = 3686223 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 22000000 status = active finish # This is the base hex mesh refined with octree read avs hex_lvl3.inp mo geniee finish cmo/status/brief The current-mesh-object(CMO) is: mo 1 Mesh Object name: mo_merge number of nodes = 634148 number of elements = 3686223 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 22000000 status = inactive 2 Mesh Object name: mo number of nodes = 155546 number of elements = 132980 dimensions geometry = 3 element type = hex dimensions topology = 3 8 nodes 6 faces 12 edges boundary flag = 16000000 status = active finish # This is the tet top cut from 25m mesh read avs tet_top.inp mo2 geniee finish cmo/status/brief The current-mesh-object(CMO) is: mo2 1 Mesh Object name: mo_merge number of nodes = 634148 number of elements = 3686223 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 22000000 status = inactive 2 Mesh Object name: mo number of nodes = 155546 number of elements = 132980 dimensions geometry = 3 element type = hex dimensions topology = 3 8 nodes 6 faces 12 edges boundary flag = 16000000 status = inactive 3 Mesh Object name: mo2 number of nodes = 5330 number of elements = 23954 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 16000000 status = active finish cmo printatt mo -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 0.000000000E+00 1.000000000E+03 1.000000000E+03 155546 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 155546 zic 0.000000000E+00 6.000000000E+02 6.000000000E+02 155546 cmo printatt mo2 -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 2.500000000E+02 7.000000000E+02 4.500000000E+02 5330 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 5330 zic 5.000000000E+02 6.000000000E+02 1.000000000E+02 5330 cmo printatt mo_merge -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 0.000000000E+00 1.000000000E+03 1.000000000E+03 634148 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 634148 zic 0.000000000E+00 6.000000000E+02 6.000000000E+02 634148 # Copy all points to new mo and connect delaunay # NOTE we are reading points from top, not elements cmo create mo_tet///tet copypts/mo_tet/mo2 copypts added nodes 1 to 5330 copypts mo_tet mo/0,0/ copypts added nodes 5331 to 160876 dump gmv tmp_points_all.gmv mo_tet cmo/modatt/-def-/-def-/ioflag/x finish cmo select mo_tet cmo printatt mo_tet -xyz- minmax Warning: nelements = 0 ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 0.000000000E+00 1.000000000E+03 1.000000000E+03 160876 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 160876 zic 0.000000000E+00 6.000000000E+02 6.000000000E+02 160876 cmo/setatt/mo_tet/itetclr 1 Warning: nelements = 0 0 length attribute: itetclr cmo/setatt/mo_tet/imt 1 Warning: nelements = 0 160876 values reset for attribute imt1 # remove duplicate points filter 1,0,0 FILTER:Use internal epsilonl value = 0.341111394591E-09 FILTER:Dudding 293 duplicate points. rmpoint compress 293 points removed and 0 elements removed. The new point count is 160583. The new element count is 0. dump gmv points_for_tet.gmv mo_tet cmo/modatt/-def-/-def-/ioflag/x finish dump avs2 points_for_tet.inp mo_tet cmo/modatt/-def-/-def-/ioflag/x finish connect noadd Coordinates of enclosing tetrahedron are: -0.75000D+03 -0.75000D+03 -0.12000D+04 -0.75000D+03 0.42500D+04 0.30000D+03 -0.75000D+03 -0.75000D+03 0.18000D+04 0.42500D+04 0.50000D+03 0.30000D+03 ntets= 32167 node index= 5000 nlstfail= 0 %done= 3.34 ntets= 64728 node index= 10000 nlstfail= 0 %done= 6.72 ntets= 98064 node index= 15000 nlstfail= 0 %done= 10.18 ntets= 131621 node index= 20000 nlstfail= 0 %done= 13.66 ntets= 164288 node index= 25000 nlstfail= 0 %done= 17.05 ntets= 197078 node index= 30000 nlstfail= 0 %done= 20.45 ntets= 228482 node index= 35000 nlstfail= 0 %done= 23.71 ntets= 259589 node index= 40000 nlstfail= 0 %done= 26.94 ntets= 290326 node index= 45000 nlstfail= 0 %done= 30.13 ntets= 320259 node index= 50000 nlstfail= 0 %done= 33.24 ntets= 350094 node index= 55000 nlstfail= 0 %done= 36.33 ntets= 379758 node index= 60000 nlstfail= 0 %done= 39.41 ntets= 409405 node index= 65000 nlstfail= 0 %done= 42.49 ntets= 439197 node index= 70000 nlstfail= 0 %done= 45.58 ntets= 469146 node index= 75000 nlstfail= 0 %done= 48.69 ntets= 498563 node index= 80000 nlstfail= 0 %done= 51.74 ntets= 528513 node index= 85000 nlstfail= 0 %done= 54.85 ntets= 557367 node index= 90000 nlstfail= 0 %done= 57.85 ntets= 586457 node index= 95000 nlstfail= 0 %done= 60.87 ntets= 615988 node index= 100000 nlstfail= 0 %done= 63.93 ntets= 645815 node index= 105000 nlstfail= 0 %done= 67.03 ntets= 675824 node index= 110000 nlstfail= 0 %done= 70.14 ntets= 705058 node index= 115000 nlstfail= 0 %done= 73.18 ntets= 734218 node index= 120000 nlstfail= 0 %done= 76.20 ntets= 762626 node index= 125000 nlstfail= 0 %done= 79.15 ntets= 791242 node index= 130000 nlstfail= 0 %done= 82.12 ntets= 818427 node index= 135000 nlstfail= 0 %done= 84.94 ntets= 844788 node index= 140000 nlstfail= 0 %done= 87.68 ntets= 870014 node index= 145000 nlstfail= 0 %done= 90.30 ntets= 893625 node index= 150000 nlstfail= 0 %done= 92.75 ntets= 915820 node index= 155000 nlstfail= 0 %done= 95.05 ntets= 937122 node index= 160000 nlstfail= 0 %done= 97.26 The mesh is now complete! LaGriT FINISH: connect rmpoint compress 0 points removed and 0 elements removed. The new point count is 160583. The new element count is 926401. quality epsilonl, epsilonaspect: 3.4111139E-10 3.9690693E-29 -------------------------------------------- elements with aspect ratio < .01: 5316 elements with aspect ratio b/w .01 and .02: 2 elements with aspect ratio b/w .02 and .05: 6 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 165 elements with aspect ratio b/w .2 and .5 : 4485 elements with aspect ratio b/w .5 and 1. : 916427 min aspect ratio = 0.2180E-03 max aspect ratio = 0.1000E+01 epsilonvol: 1.3322676E-04 --------------------------------------- element volumes b/w 0.3756E+01 and 0.2106E+02: 88 element volumes b/w 0.2106E+02 and 0.1181E+03: 746446 element volumes b/w 0.1181E+03 and 0.6624E+03: 108782 element volumes b/w 0.6624E+03 and 0.3715E+04: 43368 element volumes b/w 0.3715E+04 and 0.2083E+05: 27717 min volume = 3.7560096E+00 max volume = 2.0833333E+04 ----------------------------------------------------------- 926401 total elements evaluated. resetpts itp geniee finish dump gmv tet_connect.gmv mo_tet cmo/modatt/-def-/-def-/ioflag/x finish dump avs tet_connect.inp mo_tet cmo/modatt/-def-/-def-/ioflag/x finish # find tets formed outside boundaries define IFLAG 10 interpolate/map/mo_tet itetclr 1,0,0/mo_merge itetclr/IFLAG elements in indexed set = 926401 INTRP METHOD: map FUNCTION: user TIEBREAKER: maxtie FLAG: user 0.100000000000E+02 cmo select mo_merge finish cmo kdtree build cmo/addatt//v2/INT/scalar/scalar/constant/temporary//2.0 finish cmo/addatt//linkt/VINT/v2/nelements//temporary/x/0.0 finish cmo/addatt//v12/INT/scalar/scalar/constant/temporary//12.0 finish cmo/addatt//sbox/VDOUBLE/v12/nelements/linear/temporary/x/0.0 finish finish Sink point Elems Searched Elements Found Percent Done 18529 54070 18414 2 % 37057 111884 36942 4 % 55585 169658 55361 6 % 74113 220120 73680 8 % 92641 278973 92200 10 % 111169 337370 110728 12 % 129697 389041 128964 14 % 148225 445250 147455 16 % 166753 505258 165983 18 % 185281 561492 184327 20 % 203809 612678 202698 22 % 222337 673537 221215 24 % 240865 733554 239634 26 % 259393 781514 257888 28 % 277921 842051 276408 30 % 296449 900502 294931 32 % 314977 949023 313171 34 % 333505 1008242 331659 36 % 352033 1067350 350187 38 % 370561 1117249 368476 40 % 389089 1177129 386956 42 % 407617 1236476 405484 44 % 426145 1286354 423760 46 % 444673 1346416 442253 48 % 463201 1407737 460781 50 % 481729 1457262 479059 52 % 500257 1517353 497543 54 % 518785 1577750 516071 56 % 537313 1627302 534358 58 % 555841 1686641 552848 60 % 574369 1746181 571376 62 % 592897 1795770 589620 64 % 611425 1855694 608107 66 % 629953 1915143 626635 68 % 648481 1966016 644909 70 % 667009 2024687 663390 72 % 685537 2085104 681918 74 % 704065 2135892 700195 76 % 722593 2194188 718666 78 % 741121 2254138 737194 80 % 759649 2303879 755417 82 % 778177 2362593 773903 84 % 796705 2423038 792431 86 % 815233 2471694 810606 88 % 833761 2529593 829134 90 % 852289 2584879 847540 92 % 870817 2633042 865849 94 % 889345 2686520 884377 96 % 907873 2732311 902561 98 % 926401 2781976 921085 Total cmo select mo_tet finish WARNING: Sink points not inside source grid: 5316 Total sink points marked: 5316 with flag: 10.00 intrp/map done. cmo kdtree release cmo/DELATT/mo_merge/sbox finish cmo/DELATT/mo_merge/v12 finish cmo/DELATT/mo_merge/linkt finish cmo/DELATT/mo_merge/v2 finish finish dump gmv tet_connect.gmv mo_tet cmo/modatt/-def-/-def-/ioflag/x finish dump avs tet_connect.inp mo_tet cmo/modatt/-def-/-def-/ioflag/x finish cmo select mo_tet eltset edel itetclr eq IFLAG THE ELTSET edel HAS 5316 ELEMENTS rmpoint element eltset,get,edel Dudded 0 points plus parents, and 5316 elements removed. rmpoint compress 0 points removed and 0 elements removed. The new point count is 160583. The new element count is 921085. resetpts itp geniee finish cmo setatt mo_tet ipolydat no 1 values reset for attribute ipolydat dump gmv tet_final.gmv mo_tet cmo/modatt/-def-/-def-/ioflag/x finish dump avs tet_final.inp mo_tet cmo/modatt/-def-/-def-/ioflag/x finish # CHECK ccoef cmo setatt mo_tet itetclr 1 921085 values reset for attribute itetclr cmo setatt mo_tet imt 1 160583 values reset for attribute imt1 resetpts itp geniee finish dump stor tmp mo_tet *** Construct and Compress Sparse Matrix:3D *** *** Compress Area Coefficient Values *** AMatbld3d_stor: *****Zero Negative Coefficients ****** AMatbld3d_stor: Number of 'zero' (< 10e-8 *max) coefs 0 AMatbld3d_stor: npoints = 160583 ncoefs = 1281451 AMatbld3d_stor: Number of unique coefs = 632 AMatbld3d_stor: Maximum num. connections to a node = 23 AMatbld3d_stor: Volume min = 6.1035156E+01 AMatbld3d_stor: Volume max = 1.2500000E+05 AMatbld3d_stor: Total Volume: 5.7500000E+08 AMatbld3d_stor: abs(Aij/xij) min = 0.0000000E+00 AMatbld3d_stor: abs(Aij/xij) max = 1.3333333E+02 AMatbld3d_stor: (Aij/xij) max = 0.0000000E+00 AMatbld3d_stor: (Aij/xij) min = -1.3333333E+02 AMatbld3d_stor Matrix coefficient values stored as scalar area/distance AMatbld3d_stor Matrix compression used for graph and coefficient values ascii STOR file written with name tmp.stor *** SPARSE COEFFICIENT MATRIX _astor SUCCESSFUL *** 3D Matrix Coefficient file written with name tmp.stor dump gmv tet_ccoef.gmv mo_tet cmo/modatt/-def-/-def-/ioflag/x finish finish LaGriT successfully completed