test11.di

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

Thu Aug 26 16:03:23 EDT 1999
-----> # QMG test 11: miscellaneous tests.


Thu Aug 26 16:03:23 EDT 1999
-----> 
set show 0

0
Thu Aug 26 16:03:23 EDT 1999
-----> if {[llength [info globals interactive]]} {
    set show 1
}


Thu Aug 26 16:03:23 EDT 1999
-----> 

# test gm_addpropval & size control property


Thu Aug 26 16:03:23 EDT 1999
-----> 
# make a square


Thu Aug 26 16:03:23 EDT 1999
-----> gmset gradedsquare [gm_cpoly {{0 0} {0 1} {1 1} {1 0}} {0 1 1 1}]


Thu Aug 26 16:03:23 EDT 1999
-----> gmset gradedsquare \
	[gm_addpropval $gradedsquare e1 {sizecontrol} {"(formula 1.005-%1)"}]


Thu Aug 26 16:03:23 EDT 1999
-----> 
set pr1 [gm_lookup_prop [gm_obj2list $gradedsquare] 1 0 sizecontrol]

(formula 1.005-%1)
Thu Aug 26 16:03:23 EDT 1999
-----> if {$show} {
    puts "prop = $pr1"
}


Thu Aug 26 16:03:23 EDT 1999
-----> 

gmset mesh [gmmeshgen $gradedsquare show $show]


Thu Aug 26 16:03:24 EDT 1999
-----> 
if {$show} {
    gmviz $mesh
}


Thu Aug 26 16:03:24 EDT 1999
-----> 
set asp [gmchecktri $gradedsquare $mesh]

Maximum aspect ratio =        11.0909 achieved in
simplex #72 of topological entity mregion (2:0) which has vertices 20 21 55 
Maximum global side length =  0.672681
Minimum global altitude =     0.00146843
Number of nodes = 62 number of elements = 95
11.090909090909712
Thu Aug 26 16:03:24 EDT 1999
-----> 
if {[llength [info globals meshsizesum]]} {
    gmset {numvtx numelt} [gmmeshsize $mesh]
    global meshsizesum
    set meshsizesum [expr $meshsizesum + $numvtx]
    global aspprod
    set aspprod [expr $aspprod * $asp]
}

72499010538872080.0
Thu Aug 26 16:03:24 EDT 1999
-----> 
# test gmboundary


Thu Aug 26 16:03:24 EDT 1999
-----> 
gmset m2 [gmboundary $mesh]


Thu Aug 26 16:03:24 EDT 1999
-----> 
if {[llength [info globals meshsizesum]]} {
    gmset {numvtx numelt} [gmmeshsize $m2]
    set meshsizesum [expr $meshsizesum + $numvtx]
}

4564
Thu Aug 26 16:03:24 EDT 1999
-----> 
# test gmq11read


Thu Aug 26 16:03:24 EDT 1999
-----> 
global qmg_data


Thu Aug 26 16:03:24 EDT 1999
-----> gmset tubeobj [gmq11read $qmg_data/tubeobj.brp]


Thu Aug 26 16:03:25 EDT 1999
-----> gmset m [gmmeshgen $tubeobj show $show]


Thu Aug 26 16:03:58 EDT 1999
-----> gmset {numvtx numelt} [gmmeshsize $m]


Thu Aug 26 16:03:58 EDT 1999
-----> 

# test xdr read/write


Thu Aug 26 16:03:58 EDT 1999
-----> 
gmxdr_write $tubeobj tubeobj.xdrgeo


Thu Aug 26 16:03:58 EDT 1999
-----> gmset tubeobj2 [gmxdr_read tubeobj.xdrgeo]


Thu Aug 26 16:03:58 EDT 1999
-----> 
gmxdr_write $m tubeobj.xdrmsh


Thu Aug 26 16:03:59 EDT 1999
-----> gmset m2 [gmxdr_read tubeobj.xdrmsh]


Thu Aug 26 16:03:59 EDT 1999
-----> set asp [gmchecktri $tubeobj2 $m2]

Maximum aspect ratio =        121.688 achieved in
simplex #16005 of topological entity 3:0 (3:0) which has vertices 1928 1816 1960 3859 
Maximum global side length =  0.346472
Minimum global altitude =     0.000619213
Number of nodes = 3909 number of elements = 16497
121.68843227020716
Thu Aug 26 16:04:18 EDT 1999
-----> 
if {[llength [info globals meshsizesum]]} {
    set meshsizesum [expr $meshsizesum + $numvtx]
    set aspprod [expr $aspprod * $asp]
}

8.8222909336165704e+18
Thu Aug 26 16:04:18 EDT 1999
-----> 


# ------------------------------------------------------------------


Thu Aug 26 16:04:18 EDT 1999
-----> # Copyright (c) 1999 by Cornell University.  All rights reserved


Thu Aug 26 16:04:18 EDT 1999
-----> # See the accompanying file 'Copyright' for authorship information,


Thu Aug 26 16:04:18 EDT 1999
-----> # the terms of the license governing this software, and disclaimers


Thu Aug 26 16:04:18 EDT 1999
-----> # concerning this software.


Thu Aug 26 16:04:18 EDT 1999
-----> # ------------------------------------------------------------------


Thu Aug 26 16:04:18 EDT 1999
-----> # This file is part of the QMG software.  


Thu Aug 26 16:04:18 EDT 1999
-----> # Version 2.0 of QMG, release date RELDATE.


Thu Aug 26 16:04:18 EDT 1999
-----> # ------------------------------------------------------------------


Thu Aug 26 16:04:18 EDT 1999