test2.di

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

Thu Sep 02 17:52:26 Eastern Daylight Time 1999
-----> ## QMG test 2 This test generates a torus made up of 36


Thu Sep 02 17:52:26 Eastern Daylight Time 1999
-----> ## biquadratic patches and generates a mesh for


Thu Sep 02 17:52:26 Eastern Daylight Time 1999
-----> ## it.  The mesh size is restricted to 0.5


Thu Sep 02 17:52:26 Eastern Daylight Time 1999
-----> 
gmset t [gmmake_torus 1 2 6]


Thu Sep 02 17:52:26 Eastern Daylight Time 1999
-----> gmsummary $t

object type =                  brep
intrinsic dimension =          3
embedded dimension =           3
number of brep prop/vals =     0
global id =                    [none listed]
number of control points =     144
number of top vertices =       0
number of top edges =          0
number of bezier curves =      0
max degree among curves =      0
number of topo. surfaces =     1
number of triangle patches =   0
max degree among triangles =   0
number of quad patches =       36
max u-degree among quads =     2
max v-degree among quads =     2
max total degree among quads = 4
number of chambers =           1

Thu Sep 02 17:52:26 Eastern Daylight Time 1999
-----> set show 0

0
Thu Sep 02 17:52:26 Eastern Daylight Time 1999
-----> if {[llength [info globals interactive]]} {
    set show 1
}


Thu Sep 02 17:52:26 Eastern Daylight Time 1999
-----> gmset m [gmmeshgen $t size "(const .5)" curve "(const 1)" show $show]


Thu Sep 02 18:02:04 Eastern Daylight Time 1999
-----> gmset {numvtx numelt} [gmmeshsize $m]


Thu Sep 02 18:02:04 Eastern Daylight Time 1999
-----> set asp [gmchecktri $t $m 1 1e-10]

Maximum aspect ratio =        33.6416 achieved in
simplex #1876 of topological entity torus (3:0) which has vertices 96 241 561 584 
Maximum global side length =  0.537516
Minimum global altitude =     0.0127932
Number of nodes = 633 number of elements = 2421
33.641571813973883
Thu Sep 02 18:02:04 Eastern Daylight Time 1999
-----> 
if {[llength [info globals aspprod]]} {
    global aspprod
    set aspprod [expr $aspprod * $asp]
    global meshsizesum
    set meshsizesum [expr $meshsizesum + $numvtx]
}

892
Thu Sep 02 18:02:04 Eastern Daylight Time 1999
-----> 
# ------------------------------------------------------------------


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


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


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


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


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


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


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


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


Thu Sep 02 18:02:04 Eastern Daylight Time 1999