* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Jun 19 09:11:22 * * * 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 read avs hex_lvl3.inp mo geniee finish cmo/status/brief The current-mesh-object(CMO) is: mo 1 Mesh Object name: mo number of nodes = 631401 number of elements = 549456 dimensions geometry = 3 element type = hex dimensions topology = 3 8 nodes 6 faces 12 edges boundary flag = 22000000 status = active finish define CMO mo # infile subset.mlgi read avs tri_top_extrude.inp mo2 geniee finish cmo/status/brief The current-mesh-object(CMO) is: mo2 1 Mesh Object name: mo number of nodes = 631401 number of elements = 549456 dimensions geometry = 3 element type = hex dimensions topology = 3 8 nodes 6 faces 12 edges boundary flag = 22000000 status = inactive 2 Mesh Object name: mo2 number of nodes = 172 number of elements = 110 dimensions geometry = 3 element type = pri dimensions topology = 3 6 nodes 5 faces 9 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 631401 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 631401 zic 0.000000000E+00 6.000000000E+02 6.000000000E+02 631401 cmo printatt mo2 -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 3.000000000E+02 7.000000000E+02 4.000000000E+02 172 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 172 zic 5.500000000E+02 6.000000000E+02 5.000000000E+01 172 # Create volume mesh that can be used to find # delaunay tets that are created outside mesh boundaries addmesh merge/mo_merge/mo mo2 cmo setatt mo_merge ipolydat no 1 values reset for attribute ipolydat dump gmv mesh_merge.gmv mo_merge cmo/modatt/-def-/-def-/ioflag/x finish dump avs mesh_merge.inp mo_merge cmo/modatt/-def-/-def-/ioflag/x finish # Copy all points to new mo and connect delaunay # NOTE we are reading points from top, not elements cmo create mo_tet///tet read avs tri_top_prism_pts.inp mo_tet cmo/status/brief The current-mesh-object(CMO) is: mo_tet 1 Mesh Object name: mo number of nodes = 631401 number of elements = 549456 dimensions geometry = 3 element type = hex dimensions topology = 3 8 nodes 6 faces 12 edges boundary flag = 22000000 status = inactive 2 Mesh Object name: mo2 number of nodes = 172 number of elements = 110 dimensions geometry = 3 element type = pri dimensions topology = 3 6 nodes 5 faces 9 edges boundary flag = 16000000 status = inactive 3 Mesh Object name: mo_merge number of nodes = 631573 number of elements = 549566 dimensions geometry = 3 element type = hyb dimensions topology = 3 10 nodes 10 faces 12 edges boundary flag = 22000000 status = inactive 4 Mesh Object name: mo_tet number of nodes = 3526 number of elements = 0 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 16000000 status = active finish define CMO mo_tet # infile subset.mlgi copypts mo_tet mo/0,0/ copypts added nodes 3527 to 634927 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 634927 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 634927 zic 0.000000000E+00 6.000000000E+02 6.000000000E+02 634927 cmo/setatt/mo_tet/itetclr 1 Warning: nelements = 0 0 length attribute: itetclr cmo/setatt/mo_tet/imt 1 Warning: nelements = 0 634927 values reset for attribute imt1 # remove duplicate points filter 1,0,0 FILTER:Use internal epsilonl value = 0.341111394591E-09 FILTER:Dudding 779 duplicate points. rmpoint compress 779 points removed and 0 elements removed. The new point count is 634148. 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= 31321 node index= 5000 nlstfail= 0 %done= 0.82 ntets= 64204 node index= 10000 nlstfail= 0 %done= 1.69 ntets= 96624 node index= 15000 nlstfail= 0 %done= 2.54 ntets= 129528 node index= 20000 nlstfail= 0 %done= 3.40 ntets= 161969 node index= 25000 nlstfail= 0 %done= 4.26 ntets= 194519 node index= 30000 nlstfail= 0 %done= 5.11 ntets= 226978 node index= 35000 nlstfail= 0 %done= 5.97 ntets= 259065 node index= 40000 nlstfail= 0 %done= 6.81 ntets= 293129 node index= 45000 nlstfail= 0 %done= 7.70 ntets= 326144 node index= 50000 nlstfail= 0 %done= 8.57 ntets= 359162 node index= 55000 nlstfail= 0 %done= 9.44 ntets= 392407 node index= 60000 nlstfail= 0 %done= 10.31 ntets= 425186 node index= 65000 nlstfail= 0 %done= 11.17 ntets= 457737 node index= 70000 nlstfail= 0 %done= 12.03 ntets= 490503 node index= 75000 nlstfail= 0 %done= 12.89 ntets= 525620 node index= 80000 nlstfail= 0 %done= 13.81 ntets= 559261 node index= 85000 nlstfail= 0 %done= 14.70 ntets= 591879 node index= 90000 nlstfail= 0 %done= 15.56 ntets= 624408 node index= 95000 nlstfail= 0 %done= 16.41 ntets= 656925 node index= 100000 nlstfail= 0 %done= 17.27 ntets= 688801 node index= 105000 nlstfail= 0 %done= 18.10 ntets= 721399 node index= 110000 nlstfail= 0 %done= 18.96 ntets= 756095 node index= 115000 nlstfail= 0 %done= 19.87 ntets= 789274 node index= 120000 nlstfail= 0 %done= 20.74 ntets= 821821 node index= 125000 nlstfail= 0 %done= 21.60 ntets= 852961 node index= 130000 nlstfail= 0 %done= 22.42 ntets= 883973 node index= 135000 nlstfail= 0 %done= 23.23 ntets= 915004 node index= 140000 nlstfail= 0 %done= 24.05 ntets= 946356 node index= 145000 nlstfail= 0 %done= 24.87 ntets= 977191 node index= 150000 nlstfail= 0 %done= 25.68 ntets= 1009322 node index= 155000 nlstfail= 0 %done= 26.53 ntets= 1042416 node index= 160000 nlstfail= 0 %done= 27.40 ntets= 1073322 node index= 165000 nlstfail= 0 %done= 28.21 ntets= 1104124 node index= 170000 nlstfail= 0 %done= 29.02 ntets= 1134953 node index= 175000 nlstfail= 0 %done= 29.83 ntets= 1165554 node index= 180000 nlstfail= 0 %done= 30.63 ntets= 1196791 node index= 185000 nlstfail= 0 %done= 31.45 ntets= 1226823 node index= 190000 nlstfail= 0 %done= 32.24 ntets= 1258103 node index= 195000 nlstfail= 0 %done= 33.07 ntets= 1288738 node index= 200000 nlstfail= 0 %done= 33.87 ntets= 1318782 node index= 205000 nlstfail= 0 %done= 34.66 ntets= 1348880 node index= 210000 nlstfail= 0 %done= 35.45 ntets= 1378899 node index= 215000 nlstfail= 0 %done= 36.24 ntets= 1408973 node index= 220000 nlstfail= 0 %done= 37.03 ntets= 1438730 node index= 225000 nlstfail= 0 %done= 37.81 ntets= 1469810 node index= 230000 nlstfail= 0 %done= 38.63 ntets= 1500683 node index= 235000 nlstfail= 0 %done= 39.44 ntets= 1530853 node index= 240000 nlstfail= 0 %done= 40.23 ntets= 1561547 node index= 245000 nlstfail= 0 %done= 41.04 ntets= 1592206 node index= 250000 nlstfail= 0 %done= 41.85 ntets= 1622750 node index= 255000 nlstfail= 0 %done= 42.65 ntets= 1653048 node index= 260000 nlstfail= 0 %done= 43.45 ntets= 1682967 node index= 265000 nlstfail= 0 %done= 44.23 ntets= 1714036 node index= 270000 nlstfail= 0 %done= 45.05 ntets= 1744672 node index= 275000 nlstfail= 0 %done= 45.85 ntets= 1775143 node index= 280000 nlstfail= 0 %done= 46.65 ntets= 1805698 node index= 285000 nlstfail= 0 %done= 47.46 ntets= 1836156 node index= 290000 nlstfail= 0 %done= 48.26 ntets= 1866703 node index= 295000 nlstfail= 0 %done= 49.06 ntets= 1897920 node index= 300000 nlstfail= 0 %done= 49.88 ntets= 1929072 node index= 305000 nlstfail= 0 %done= 50.70 ntets= 1959915 node index= 310000 nlstfail= 0 %done= 51.51 ntets= 1990118 node index= 315000 nlstfail= 0 %done= 52.30 ntets= 2020831 node index= 320000 nlstfail= 0 %done= 53.11 ntets= 2051369 node index= 325000 nlstfail= 0 %done= 53.91 ntets= 2082175 node index= 330000 nlstfail= 0 %done= 54.72 ntets= 2112765 node index= 335000 nlstfail= 0 %done= 55.53 ntets= 2142760 node index= 340000 nlstfail= 0 %done= 56.32 ntets= 2173109 node index= 345000 nlstfail= 0 %done= 57.11 ntets= 2203240 node index= 350000 nlstfail= 0 %done= 57.91 ntets= 2233685 node index= 355000 nlstfail= 0 %done= 58.71 ntets= 2264091 node index= 360000 nlstfail= 0 %done= 59.50 ntets= 2294604 node index= 365000 nlstfail= 0 %done= 60.31 ntets= 2324883 node index= 370000 nlstfail= 0 %done= 61.10 ntets= 2355002 node index= 375000 nlstfail= 0 %done= 61.89 ntets= 2384759 node index= 380000 nlstfail= 0 %done= 62.68 ntets= 2414358 node index= 385000 nlstfail= 0 %done= 63.45 ntets= 2443635 node index= 390000 nlstfail= 0 %done= 64.22 ntets= 2473098 node index= 395000 nlstfail= 0 %done= 65.00 ntets= 2502646 node index= 400000 nlstfail= 0 %done= 65.77 ntets= 2532058 node index= 405000 nlstfail= 0 %done= 66.55 ntets= 2561340 node index= 410000 nlstfail= 0 %done= 67.32 ntets= 2590209 node index= 415000 nlstfail= 0 %done= 68.08 ntets= 2619521 node index= 420000 nlstfail= 0 %done= 68.85 ntets= 2648528 node index= 425000 nlstfail= 0 %done= 69.61 ntets= 2677719 node index= 430000 nlstfail= 0 %done= 70.38 ntets= 2706994 node index= 435000 nlstfail= 0 %done= 71.14 ntets= 2736351 node index= 440000 nlstfail= 0 %done= 71.92 ntets= 2765574 node index= 445000 nlstfail= 0 %done= 72.68 ntets= 2791298 node index= 450000 nlstfail= 0 %done= 73.36 ntets= 2819442 node index= 455000 nlstfail= 0 %done= 74.10 ntets= 2847346 node index= 460000 nlstfail= 0 %done= 74.83 ntets= 2874829 node index= 465000 nlstfail= 0 %done= 75.56 ntets= 2902225 node index= 470000 nlstfail= 0 %done= 76.28 ntets= 2929927 node index= 475000 nlstfail= 0 %done= 77.00 ntets= 2957363 node index= 480000 nlstfail= 0 %done= 77.72 ntets= 2984909 node index= 485000 nlstfail= 0 %done= 78.45 ntets= 3009398 node index= 490000 nlstfail= 0 %done= 79.09 ntets= 3036623 node index= 495000 nlstfail= 0 %done= 79.81 ntets= 3063081 node index= 500000 nlstfail= 0 %done= 80.50 ntets= 3089439 node index= 505000 nlstfail= 0 %done= 81.20 ntets= 3115825 node index= 510000 nlstfail= 0 %done= 81.89 ntets= 3142515 node index= 515000 nlstfail= 0 %done= 82.59 ntets= 3169065 node index= 520000 nlstfail= 0 %done= 83.29 ntets= 3192908 node index= 525000 nlstfail= 0 %done= 83.92 ntets= 3219268 node index= 530000 nlstfail= 0 %done= 84.61 ntets= 3245516 node index= 535000 nlstfail= 0 %done= 85.30 ntets= 3272026 node index= 540000 nlstfail= 0 %done= 85.99 ntets= 3298588 node index= 545000 nlstfail= 0 %done= 86.69 ntets= 3325083 node index= 550000 nlstfail= 0 %done= 87.39 ntets= 3351625 node index= 555000 nlstfail= 0 %done= 88.09 ntets= 3377673 node index= 560000 nlstfail= 0 %done= 88.77 ntets= 3401031 node index= 565000 nlstfail= 0 %done= 89.39 ntets= 3426256 node index= 570000 nlstfail= 0 %done= 90.05 ntets= 3451976 node index= 575000 nlstfail= 0 %done= 90.72 ntets= 3477896 node index= 580000 nlstfail= 0 %done= 91.41 ntets= 3503853 node index= 585000 nlstfail= 0 %done= 92.09 ntets= 3529662 node index= 590000 nlstfail= 0 %done= 92.77 ntets= 3555538 node index= 595000 nlstfail= 0 %done= 93.45 ntets= 3578217 node index= 600000 nlstfail= 0 %done= 94.04 ntets= 3600663 node index= 605000 nlstfail= 0 %done= 94.63 ntets= 3622737 node index= 610000 nlstfail= 0 %done= 95.21 ntets= 3643940 node index= 615000 nlstfail= 0 %done= 95.77 ntets= 3665135 node index= 620000 nlstfail= 0 %done= 96.33 ntets= 3686635 node index= 625000 nlstfail= 0 %done= 96.89 ntets= 3707937 node index= 630000 nlstfail= 0 %done= 97.45 The mesh is now complete! LaGriT FINISH: connect rmpoint compress 0 points removed and 0 elements removed. The new point count is 634148. The new element count is 3696596. quality epsilonl, epsilonaspect: 3.4111139E-10 3.9690693E-29 -------------------------------------------- elements with aspect ratio < .01: 10275 elements with aspect ratio b/w .01 and .02: 108 elements with aspect ratio b/w .02 and .05: 8 elements with aspect ratio b/w .05 and .1 : 14 elements with aspect ratio b/w .1 and .2 : 124 elements with aspect ratio b/w .2 and .5 : 10923 elements with aspect ratio b/w .5 and 1. : 3675144 min aspect ratio = 0.1381E-03 max aspect ratio = 0.1000E+01 epsilonvol: 1.3322676E-04 --------------------------------------- element volumes b/w 0.3756E+01 and 0.1596E+02: 2922656 element volumes b/w 0.1596E+02 and 0.6784E+02: 355864 element volumes b/w 0.6784E+02 and 0.2883E+03: 111213 element volumes b/w 0.2883E+03 and 0.1225E+04: 123083 element volumes b/w 0.1225E+04 and 0.5208E+04: 183780 min volume = 3.7560096E+00 max volume = 5.2083333E+03 ----------------------------------------------------------- 3696596 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 = 3696596 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 73933 73585 73585 2 % 147865 147507 147507 4 % 221797 221436 221436 6 % 295729 294760 294760 8 % 369661 368507 368507 10 % 443593 442435 442435 12 % 517525 515837 515837 14 % 591457 589557 589557 16 % 665389 663483 663483 18 % 739321 736962 736962 20 % 813253 810645 810645 22 % 887185 884577 884577 24 % 961117 958509 958509 26 % 1035049 1031902 1031902 28 % 1108981 1105716 1105716 30 % 1182913 1179648 1179648 32 % 1256845 1253040 1253040 34 % 1330777 1326859 1326859 36 % 1404709 1400791 1400791 38 % 1478641 1474274 1474274 40 % 1552573 1548106 1548106 42 % 1626505 1622038 1622038 44 % 1700437 1695469 1695469 46 % 1774369 1769291 1769291 48 % 1848301 1843223 1843223 50 % 1922233 1916686 1916686 52 % 1996165 1990518 1990518 54 % 2070097 2064450 2064450 56 % 2144029 2137980 2137980 58 % 2217961 2211745 2211745 60 % 2291893 2285677 2285677 62 % 2365825 2359202 2359202 64 % 2439757 2432944 2432944 66 % 2513689 2506876 2506876 68 % 2587621 2580451 2580451 70 % 2661553 2654224 2654224 72 % 2735485 2728156 2728156 74 % 2809417 2801620 2801620 76 % 2883349 2875446 2875446 78 % 2957281 2949378 2949378 80 % 3031213 3022893 3022893 82 % 3105145 3096719 3096719 84 % 3179077 3170488 3170488 86 % 3253009 3244003 3244003 88 % 3326941 3317935 3317935 90 % 3400873 3391350 3391350 92 % 3474805 3465144 3465144 94 % 3548737 3539076 3539076 96 % 3622669 3612296 3612296 98 % 3696596 3686223 3686223 Total cmo select mo_tet finish WARNING: Sink points not inside source grid: 10373 Total sink points marked: 10373 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 10373 ELEMENTS rmpoint element eltset,get,edel Dudded 0 points plus parents, and 10373 elements removed. rmpoint compress 0 points removed and 0 elements removed. The new point count is 634148. The new element count is 3686223. 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 3686223 values reset for attribute itetclr cmo setatt mo_tet imt 1 634148 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 = 634148 ncoefs = 5033908 AMatbld3d_stor: Number of unique coefs = 610 AMatbld3d_stor: Maximum num. connections to a node = 20 AMatbld3d_stor: Volume min = 7.6293943E+00 AMatbld3d_stor: Volume max = 1.5625000E+04 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