test8.di

来自「算断裂的」· DI 代码 · 共 143 行

Thu Sep 02 18:10:57 Eastern Daylight Time 1999
-----> # QMG test 8: a cube in standard orientation versus


Thu Sep 02 18:10:57 Eastern Daylight Time 1999
-----> # a tilted cube.


Thu Sep 02 18:10:57 Eastern Daylight Time 1999
-----> 
global qmg_data


Thu Sep 02 18:10:57 Eastern Daylight Time 1999
-----> gmset cube [gmoffread $qmg_data/cube.off]


Thu Sep 02 18:10:57 Eastern Daylight Time 1999
-----> set show 0

0
Thu Sep 02 18:10:57 Eastern Daylight Time 1999
-----> if {[llength [info globals interactive]]} {
    set show 1
}


Thu Sep 02 18:10:57 Eastern Daylight Time 1999
-----> 
gmset mesh [gmmeshgen $cube  show $show]


Thu Sep 02 18:10:58 Eastern Daylight Time 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 Sep 02 18:10:58 Eastern Daylight Time 1999
-----> gmset {numvtx numelt} [gmmeshsize $mesh]


Thu Sep 02 18:10:58 Eastern Daylight Time 1999
-----> 
if {[llength [info globals aspprod]]} {
    global aspprod
    set aspprod [expr $aspprod * $asp]
    global meshsizesum
    set meshsizesum [expr $meshsizesum + $numvtx]
}

3214
Thu Sep 02 18:10:58 Eastern Daylight Time 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 Sep 02 18:10:58 Eastern Daylight Time 1999
-----> 
gmset tiltcube [gmapply $q $cube]


Thu Sep 02 18:10:58 Eastern Daylight Time 1999
-----> 
gmset mesh2 [gmmeshgen $tiltcube  show $show]


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> 
## Display the properties of the second mesh


Thu Sep 02 18:11:00 Eastern Daylight Time 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 Sep 02 18:11:00 Eastern Daylight Time 1999
-----> gmset {numvtx2 numelt2} [gmmeshsize $mesh2]


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> 
if {[llength [info globals aspprod]]} {
    global aspprod
    set aspprod [expr $aspprod * $asp2]
    global meshsizesum
    set meshsizesum [expr $meshsizesum + $numvtx2]
}

3442
Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> 
# ------------------------------------------------------------------


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> # Copyright (c) 1999 by Cornell University.  All rights reserved


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> # See the accompanying file 'Copyright' for authorship information,


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> # the terms of the license governing this software, and disclaimers


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> # concerning this software.


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> # ------------------------------------------------------------------


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> # This file is part of the QMG software.  


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> # Version 2.0 of QMG, release date September 3, 1999.


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> # ------------------------------------------------------------------


Thu Sep 02 18:11:00 Eastern Daylight Time 1999