* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Sep 5 13:11:48 * * * 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 # read subset truncated at top and bottom read avs tet_trunc.inp cmotet cmo/addatt/cmotet/xval/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/cmotet/nlayer/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/cmotet/imat/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish geniee finish cmo/status/brief The current-mesh-object(CMO) is: cmotet 1 Mesh Object name: cmotet number of nodes = 17200 number of elements = 94792 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 16000000 status = active finish resetpts itp geniee finish cmo printatt cmotet -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 17200 scalar 1 1 0 1 vector 3 3 0 1 nnodes 17200 17200 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 94792 94792 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 3 3 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 4 4 0 1 edges_per_element 6 6 0 1 faces_per_element 4 4 0 1 isetwd 0 0 0 17200 ialias 0 0 0 17200 imt1 1 1 0 17200 itp1 0 12 12 17200 icr1 0 0 0 17200 isn1 0 0 0 17200 xic 2.410500000E+05 2.885500000E+05 4.750000000E+04 17200 yic 4.044700000E+06 4.092200000E+06 4.750000000E+04 17200 zic -3.430279602E+03 1.523439774E+03 4.953719376E+03 17200 xtetwd 0 0 0 94792 itetclr 1 7 6 94792 itettyp 5 5 0 94792 itetoff 0 379164 379164 94792 jtetoff 0 379164 379164 94792 itet 1 17200 17199 94792x4 jtet 1 16376497 16376496 94792x4 epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 1.495638024E-08 1.495638024E-08 0.000000000E+00 1 epsilona 1.007425107E-03 1.007425107E-03 0.000000000E+00 1 epsilonv 2.481754655E+00 2.481754655E+00 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 17200 17200 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 2.410500000E+05 2.410500000E+05 0.000000000E+00 1 ymin 4.044700000E+06 4.044700000E+06 0.000000000E+00 1 zmin -3.430279602E+03 -3.430279602E+03 0.000000000E+00 1 xmax 2.885500000E+05 2.885500000E+05 0.000000000E+00 1 ymax 4.092200000E+06 4.092200000E+06 0.000000000E+00 1 zmax 1.523439774E+03 1.523439774E+03 0.000000000E+00 1 kdtree_level 0 0 0 1 max_number_sets 64 64 0 1 number_of_psets 0 0 0 1 number_of_eltsets 0 0 0 1 number_of_fsets 0 0 0 1 xval 1.000000000E+00 7.000000000E+00 6.000000000E+00 17200 nlayer 1 43 42 17200 imat 1 7 6 17200 ********* Write FEHM setup files ******** cmo/addatt/cmotet/matsav/vint/scalar/nelements/linear/permanent/afgx/0. cmo/addatt/cmotet/cc_neg/vint/scalar/nnodes/linear/permanent/afgx/0. cmo copyatt/cmotet cmotet/matsav itetclr 94792 copied from cmotet itetclr to -> cmotet matsav cmo setatt cmotet itetclr 1 94792 values reset for attribute itetclr cmo copyatt/cmotet cmotet/imt imat 17200 copied from cmotet imat to -> cmotet imt1 filter 1,0,0 FILTER:Use internal epsilonl value = 0.149563802384E-07 FILTER:Dudding 0 duplicate points. rmpoint compress 0 points removed and 0 elements removed. The new point count is 17200. The new element count is 94792. resetpts itp geniee finish quality epsilonl, epsilonaspect: 1.4956380E-08 3.3456422E-24 -------------------------------------------- elements with aspect ratio < .01: 3820 elements with aspect ratio b/w .01 and .02: 2 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 : 26842 elements with aspect ratio b/w .5 and 1. : 64128 min aspect ratio = 0.6001E-05 max aspect ratio = 0.1166E+00 epsilonvol: 2.4817547E+00 --------------------------------------- element volumes b/w 0.1473E+05 and 0.9895E+05: 67 element volumes b/w 0.9895E+05 and 0.6645E+06: 2270 element volumes b/w 0.6645E+06 and 0.4462E+07: 1485 element volumes b/w 0.4462E+07 and 0.2996E+08: 0 element volumes b/w 0.2996E+08 and 0.2012E+09: 90970 min volume = 1.4734271E+04 max volume = 2.0122611E+08 ----------------------------------------------------------- 94792 total elements evaluated. * Must be single defined mesh cmo select cmotet resetpts/parent rmpoint compress 0 points removed and 0 elements removed. The new point count is 17200. The new element count is 94792. cmo status cmotet The current-mesh-object(CMO) is: cmotet 1 Mesh Object name: cmotet number of nodes = 17200 number of elements = 94792 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 16000000 status = active NAME TYPE RANK LENGTH INTER PERSI IO VALUE 1 -def- VDOU scalar nnodes linea tempo agl 0.000E+00 2 scalar INT scalar scalar const perma l 1.000E+00 3 vector INT scalar scalar const perma l 3.000E+00 4 nnodes INT scalar scalar const perma l 1.720E+04 5 nedges INT scalar scalar const perma l 0.000E+00 6 nfaces INT scalar scalar const perma l 0.000E+00 7 nelements INT scalar scalar const perma l 9.479E+04 8 mbndry INT scalar scalar const perma l 1.600E+07 9 ndimensions_top INT scalar scalar const perma l 3.000E+00 10 ndimensions_geo INT scalar scalar const perma l 3.000E+00 11 nodes_per_eleme INT scalar scalar const perma l 4.000E+00 12 edges_per_eleme INT scalar scalar const perma l 6.000E+00 13 faces_per_eleme INT scalar scalar const perma l 4.000E+00 14 isetwd VDOU scalar nnodes or perma l 0.000E+00 15 ialias VINT scalar nnodes seque perma l 0.000E+00 16 imt1 VINT scalar nnodes max perma gal 0.000E+00 17 itp1 VINT scalar nnodes min perma gal 0.000E+00 18 icr1 VINT scalar nnodes min perma gal 0.000E+00 19 isn1 VINT scalar nnodes user perma gal 0.000E+00 20 xic VDOU scalar nnodes linea perma l 0.000E+00 21 yic VDOU scalar nnodes linea perma l 0.000E+00 22 zic VDOU scalar nnodes linea perma l 0.000E+00 23 xtetwd VDOU scalar nelements or perma l 0.000E+00 24 itetclr VINT scalar nelements user perma l 0.000E+00 25 itettyp VINT scalar nelements user perma l 0.000E+00 26 itetoff VINT scalar nelements user perma l 0.000E+00 27 jtetoff VINT scalar nelements user perma l 0.000E+00 28 itet VINT nodes_per_ nelements user perma l 0.000E+00 29 jtet VINT faces_per_ nelements user perma l 0.000E+00 30 ipolydat CHAR scalar scalar const perma l yes 31 vor2d CHAR scalar scalar const perma l yes 32 vor3d CHAR scalar scalar const perma l no 33 epsilon REAL scalar scalar const perma l 1.000E-15 34 epsilonl REAL scalar scalar const perma l 1.496E-08 35 epsilona REAL scalar scalar const perma l 1.007E-03 36 epsilonv REAL scalar scalar const perma l 2.482E+00 37 ipointi INT scalar scalar const perma l 1.000E+00 38 ipointj INT scalar scalar const perma l 1.720E+04 39 idebug INT scalar scalar const perma l 0.000E+00 40 itypconv_sm INT scalar scalar const perma l 1.000E+00 41 maxiter_sm INT scalar scalar const perma l 2.500E+01 42 tolconv_sm REAL scalar scalar const perma l 1.000E+00 43 nnfreq INT scalar scalar const perma l 1.000E+00 44 ivoronoi INT scalar scalar const perma l 1.000E+00 45 iopt2to2 INT scalar scalar const perma l 2.000E+00 46 dumptype CHAR scalar scalar const perma l binary 47 velname CHAR scalar scalar const perma l vels 48 densname CHAR scalar scalar const perma l ric 49 presname CHAR scalar scalar const perma l pic 50 enername CHAR scalar scalar const perma l eic 51 xmin REAL scalar scalar const perma l 2.411E+05 52 ymin REAL scalar scalar const perma l 4.045E+06 53 zmin REAL scalar scalar const perma l -3.430E+03 54 xmax REAL scalar scalar const perma l 2.886E+05 55 ymax REAL scalar scalar const perma l 4.092E+06 56 zmax REAL scalar scalar const perma l 1.523E+03 57 kdtree_level INT scalar scalar const perma l 0.000E+00 58 max_number_sets INT scalar scalar const perma l 6.400E+01 59 number_of_psets INT scalar scalar const perma l 0.000E+00 60 number_of_eltse INT scalar scalar const perma l 0.000E+00 61 psetnames VCHA scalar max_number_sets const perma l 62 eltsetnames VCHA scalar max_number_sets const perma l 63 geom_name CHAR scalar max_number_sets const perma l -defaultge 64 fsetnames VCHA scalar max_number_sets const perma l 65 number_of_fsets INT scalar scalar const perma l 0.000E+00 66 xval VDOU scalar nnodes linea perma gxa 0.000E+00 67 nlayer VINT scalar nnodes linea perma gxa 0.000E+00 68 imat VINT scalar nnodes linea perma gxa 0.000E+00 69 matsav VINT scalar nelements linea perma afg 0.000E+00 70 cc_neg VINT scalar nnodes linea perma afg 0.000E+00 cmo setatt cmotet ipolydat no 1 values reset for attribute ipolydat dump fehm tet cmotet ascii *** Write FEHMN GEOM AND ZONE FILES *** *********dump_material_lists******** Minimum material ID value = 1 Maximum material ID value = 7 Total possible materials = 7 Material 1 has 4000 nodes. #nodes/nnodes is 0.232558146119 Material 2 has 800 nodes. #nodes/nnodes is 0.465116277337E-01 Material 3 has 400 nodes. #nodes/nnodes is 0.232558138669E-01 Material 4 has 1200 nodes. #nodes/nnodes is 0.697674453259E-01 Material 5 has 6000 nodes. #nodes/nnodes is 0.348837196827 Material 6 has 800 nodes. #nodes/nnodes is 0.465116277337E-01 Material 7 has 4000 nodes. #nodes/nnodes is 0.232558146119 *********dump_interface_list******** *********dump_multi_material_list******** Material 1 has 1522 multi-material connections. Material 2 has 3046 multi-material connections. Material 3 has 3046 multi-material connections. Material 4 has 3046 multi-material connections. Material 5 has 3048 multi-material connections. Material 6 has 3054 multi-material connections. Material 7 has 1530 multi-material connections. log/tty/off finish *********dump_outside_list******** Voronoi Areas used for outside faces. assign_quadrant_id epsilon for zero: 1.0000000E-15 cmo/addatt/cmotet/top/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face top 1 has 400 nodes. top Sum Voronoi Area_x Area_y Area_z 0.9465189E+07 0.1127235E+08 0.2256250E+10 cmo/addatt/cmotet/bottom/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face bottom 2 has 400 nodes. bottom Sum Voronoi Area_x Area_y Area_z 0.1619625E+08 0.7915136E+07 0.2256250E+10 cmo/addatt/cmotet/left_w/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face left_w 3 has 860 nodes. left_w Sum Voronoi Area_x Area_y Area_z 0.2252806E+09 0.1234406E+08 0.1187500E+09 cmo/addatt/cmotet/right_e/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face right_e 5 has 860 nodes. right_e Sum Voronoi Area_x Area_y Area_z 0.2101958E+09 0.1125071E+08 0.1187500E+09 cmo/addatt/cmotet/back_n/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face back_n 6 has 860 nodes. back_n Sum Voronoi Area_x Area_y Area_z 0.1191237E+08 0.2182378E+09 0.1187500E+09 cmo/addatt/cmotet/front_s/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face front_s 4 has 860 nodes. front_s Sum Voronoi Area_x Area_y Area_z 0.1191992E+08 0.2157814E+09 0.1187500E+09 cmo/DELATT//top finish cmo/DELATT//bottom finish cmo/DELATT//left_w finish cmo/DELATT//right_e finish cmo/DELATT//back_n finish cmo/DELATT//front_s finish log/tty/on finish *********dump_parent_list******** *** Construct and Compress Sparse Matrix:3D *** *** Compress Area Coefficient Values *** AMatbld3d_stor: *****Negative Coefficients ****** AMatbld3d_stor: Total Number of Negative Coefficients 7 AMatbld3d_stor: Number of Significant Negative Coefs 7 AMatbld3d_stor: Number of 'zero' (< 10e-8 *max) coefs 0 Negative coef at row 4140 Column 4520 with value -8.1615288E-04 Negative coef at row 7598 Column 7978 with value -2.7384189E-03 Negative coef at row 7795 Column 8175 with value -8.1615288E-04 Negative coef at row 8636 Column 9056 with value -2.7384189E-03 Negative coef at row 11845 Column 12225 with value -2.7384189E-03 Negative coef at row 14355 Column 14735 with value -4.7908997E-03 Negative coef at row 16055 Column 16435 with value -2.7384189E-03 cmo/addatt/cmotet/ccoef/VDOUBLE/scalar/nnodes/linear/permanent/afgx/0.0/ finish cmo/addatt/cmotet/ij_ccoef/VINT/scalar/nnodes/linear/permanent/afgx/0.0/ finish AMatbld3d_stor: npoints = 17200 ncoefs = 204424 AMatbld3d_stor: Number of unique coefs = 54042 AMatbld3d_stor: Maximum num. connections to a node = 15 AMatbld3d_stor: Volume min = 5.2338774E+07 AMatbld3d_stor: Volume max = 8.8327354E+08 AMatbld3d_stor: Total Volume: 1.0313790E+13 AMatbld3d_stor: abs(Aij/xij) min = 0.0000000E+00 AMatbld3d_stor: abs(Aij/xij) max = 7.8377798E+04 AMatbld3d_stor: (Aij/xij) max = 4.7908997E-03 AMatbld3d_stor: (Aij/xij) min = -7.8377798E+04 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 tet.stor *** SPARSE COEFFICIENT MATRIX _astor SUCCESSFUL *** 3D Matrix Coefficient file written with name tet.stor cmo copyatt cmotet cmotet itetclr matsav 94792 copied from cmotet matsav to -> cmotet itetclr # check to see if neg ccoefs are on outside of mesh quality epsilonl, epsilonaspect: 1.4956380E-08 3.3456422E-24 -------------------------------------------- elements with aspect ratio < .01: 3820 elements with aspect ratio b/w .01 and .02: 2 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 : 26842 elements with aspect ratio b/w .5 and 1. : 64128 min aspect ratio = 0.6001E-05 max aspect ratio = 0.1166E+00 epsilonvol: 2.4817547E+00 --------------------------------------- element volumes b/w 0.1473E+05 and 0.9895E+05: 67 element volumes b/w 0.9895E+05 and 0.6645E+06: 2270 element volumes b/w 0.6645E+06 and 0.4462E+07: 1485 element volumes b/w 0.4462E+07 and 0.2996E+08: 0 element volumes b/w 0.2996E+08 and 0.2012E+09: 90970 min volume = 1.4734271E+04 max volume = 2.0122611E+08 ----------------------------------------------------------- 94792 total elements evaluated. pset ncoef attribute ccoef 1,0,0 lt -0.001 THE PSET ncoef HAS 10 POINTS pset pout attribute itp 1,0,0 eq 10. THE PSET pout HAS 3916 POINTS pset pin not pout THE PSET pin HAS 13284 POINTS pset pBAD inter pin ncoef THE PSET pBAD HAS 9 POINTS cmo setatt cmotet cc_neg pset,get,ncoef 1 10 values reset for attribute cc_neg dump gmv tet_fehm.gmv cmotet cmo/modatt/-def-/-def-/ioflag/x finish * translate to 0,0 for better voronoi volume precision trans 1,0,0 zero cmo/addatt/cmotet/transauto/VDOUBLE/scalar/vector/user/permanent/x/0.0 finish dump stor tet cmotet ascii *** Construct and Compress Sparse Matrix:3D *** *** Compress Area Coefficient Values *** AMatbld3d_stor: *****Negative Coefficients ****** AMatbld3d_stor: Total Number of Negative Coefficients 7 AMatbld3d_stor: Number of Significant Negative Coefs 7 AMatbld3d_stor: Number of 'zero' (< 10e-8 *max) coefs 0 Negative coef at row 4140 Column 4520 with value -8.1615390E-04 Negative coef at row 7598 Column 7978 with value -2.7384173E-03 Negative coef at row 7795 Column 8175 with value -8.1615390E-04 Negative coef at row 8636 Column 9056 with value -2.7384173E-03 Negative coef at row 11845 Column 12225 with value -2.7384173E-03 Negative coef at row 14355 Column 14735 with value -4.7909096E-03 Negative coef at row 16055 Column 16435 with value -2.7384173E-03 AMatbld3d_stor: npoints = 17200 ncoefs = 204424 AMatbld3d_stor: Number of unique coefs = 54042 AMatbld3d_stor: Maximum num. connections to a node = 15 AMatbld3d_stor: Volume min = 5.2338774E+07 AMatbld3d_stor: Volume max = 8.8327353E+08 AMatbld3d_stor: Total Volume: 1.0313790E+13 AMatbld3d_stor: abs(Aij/xij) min = 0.0000000E+00 AMatbld3d_stor: abs(Aij/xij) max = 7.8377798E+04 AMatbld3d_stor: (Aij/xij) max = 4.7909096E-03 AMatbld3d_stor: (Aij/xij) min = -7.8377798E+04 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 tet.stor *** SPARSE COEFFICIENT MATRIX _astor SUCCESSFUL *** 3D Matrix Coefficient file written with name tet.stor # Extract surface mesh for viewing voronoi polygons cmo select cmotet extract/surfmesh/1 0 0/cmos/cmotet/external cmo/addatt/cmos/itetclr0/vint/scalar/nelements////0 finish cmo/addatt/cmos/itetclr1/vint/scalar/nelements////0 finish cmo/addatt/cmos/idnode0/vint/scalar/nnodes////0 finish cmo/addatt/cmos/idelem0/vint/scalar/nelements////0 finish cmo/addatt/cmos/idelem1/vint/scalar/nelements////0 finish cmo/addatt/cmos/facecol/vint/scalar/nelements////0 finish cmo/addatt/cmos/idface0/vint/scalar/nelements////0 finish cmo/addatt/cmos/idface1/vint/scalar/nelements////0 finish cmo/addatt/cmos/xval/VDOUBLE/scalar/nnodes/linear/permanent/gxaf finish cmo/addatt/cmos/nlayer/VINT/scalar/nnodes/linear/permanent/gxaf finish cmo/addatt/cmos/imat/VINT/scalar/nnodes/linear/permanent/gxaf finish cmo/addatt/cmos/cc_neg/VINT/scalar/nnodes/linear/permanent/afgx finish cmo/addatt/cmos/ccoef/VDOUBLE/scalar/nnodes/linear/permanent/afgx finish cmo/addatt/cmos/ij_ccoef/VINT/scalar/nnodes/linear/permanent/afgx finish cmo/addatt/cmos/ncon50/int/scalar/scalar/constant/permanent/x/0 finish cmo/addatt/cmos/nconbnd/int/scalar/scalar/constant/permanent/x/0 finish cmo/addatt/cmos/icontab/vint/scalar/ncon50/constant/permanent/x/0 finish resetpts/itp geniee finish finish cmo copyatt/cmos cmos/itetclr itetclr1 7828 copied from cmos itetclr1 to -> cmos itetclr cmo setatt cmos ipolydat yes 1 values reset for attribute ipolydat dump gmv tet_surfmesh.gmv cmos cmo/modatt/-def-/-def-/ioflag/x finish finish LaGriT successfully completed