test8.di

算断裂的

Thu Aug 26 15:57:32 EDT 1999
-----> # QMG test 8: a cube in standard orientation versus


Thu Aug 26 15:57:32 EDT 1999
-----> # a tilted cube.


Thu Aug 26 15:57:32 EDT 1999
-----> 
global qmg_data


Thu Aug 26 15:57:32 EDT 1999
-----> gmset cube [gmoffread $qmg_data/cube.off]


Thu Aug 26 15:57:32 EDT 1999
-----> set show 0

0
Thu Aug 26 15:57:32 EDT 1999
-----> if {[llength [info globals interactive]]} {
    set show 1
}


Thu Aug 26 15:57:32 EDT 1999
-----> 
gmset mesh [gmmeshgen $cube  show $show]


Thu Aug 26 15:57:33 EDT 1999
-----> 
set asp [gmchecktri $cube $mesh]

Maximum aspect ratio =        2.70243 achieved in
simplex #12 of topological entity region (3:0) which has vertices 13 4 21 26 
Maximum global side length =  1.9105
Minimum global altitude =     0.597927
Number of nodes = 27 number of elements = 48
2.7024294625392167
Thu Aug 26 15:57:33 EDT 1999
-----> gmset {numvtx numelt} [gmmeshsize $mesh]


Thu Aug 26 15:57:33 EDT 1999
-----> 
if {[llength [info globals aspprod]]} {
    global aspprod
    set aspprod [expr $aspprod * $asp]
    global meshsizesum
    set meshsizesum [expr $meshsizesum + $numvtx]
}

3214
Thu Aug 26 15:57:33 EDT 1999
-----> 
gmset q {
    {  -0.84910409704140  -0.40704628877252   0.33665345859509 0}
    {   0.23026379507690  -0.85881459325243  -0.45762001605420 0}
    {   0.47539543231325  -0.31104792751790   0.82295101295431 0}
}


Thu Aug 26 15:57:33 EDT 1999
-----> 
gmset tiltcube [gmapply $q $cube]


Thu Aug 26 15:57:34 EDT 1999
-----> 
gmset mesh2 [gmmeshgen $tiltcube  show $show]


Thu Aug 26 15:57:37 EDT 1999
-----> 
## Display the properties of the second mesh


Thu Aug 26 15:57:37 EDT 1999
-----> 
set asp2 [gmchecktri $tiltcube $mesh2]

Maximum aspect ratio =        78.5844 achieved in
simplex #606 of topological entity region (3:0) which has vertices 108 6 110 195 
Maximum global side length =  1.55283
Minimum global altitude =     0.00634797
Number of nodes = 228 number of elements = 867
78.584388551085439
Thu Aug 26 15:57:38 EDT 1999
-----> gmset {numvtx2 numelt2} [gmmeshsize $mesh2]


Thu Aug 26 15:57:38 EDT 1999
-----> 
if {[llength [info globals aspprod]]} {
    global aspprod
    set aspprod [expr $aspprod * $asp2]
    global meshsizesum
    set meshsizesum [expr $meshsizesum + $numvtx2]
}

3442
Thu Aug 26 15:57:38 EDT 1999
-----> 
# ------------------------------------------------------------------


Thu Aug 26 15:57:38 EDT 1999
-----> # Copyright (c) 1999 by Cornell University.  All rights reserved


Thu Aug 26 15:57:38 EDT 1999
-----> # See the accompanying file 'Copyright' for authorship information,


Thu Aug 26 15:57:38 EDT 1999
-----> # the terms of the license governing this software, and disclaimers


Thu Aug 26 15:57:38 EDT 1999
-----> # concerning this software.


Thu Aug 26 15:57:38 EDT 1999
-----> # ------------------------------------------------------------------


Thu Aug 26 15:57:38 EDT 1999
-----> # This file is part of the QMG software.  


Thu Aug 26 15:57:38 EDT 1999
-----> # Version 2.0 of QMG, release date RELDATE.


Thu Aug 26 15:57:38 EDT 1999
-----> # ------------------------------------------------------------------


Thu Aug 26 15:57:38 EDT 1999