test4.di

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

Thu Aug 26 15:34:19 EDT 1999
-----> ## Test 4: mesh a cylinder with two small spheres removed (like bubbles)


Thu Aug 26 15:34:19 EDT 1999
-----> 

global qmg_data


Thu Aug 26 15:34:19 EDT 1999
-----> 
gmset cyl1 [gmmake_cyl 4]


Thu Aug 26 15:34:19 EDT 1999
-----> gmset s [gm_read $qmg_data/sphere3.geo]


Thu Aug 26 15:34:19 EDT 1999
-----> gmset bubble1 [gmapply \
	[gmcompose [gmtranslate -.02 0 0] [gmdilate .01 .01 .01]] \
	$s]


Thu Aug 26 15:34:19 EDT 1999
-----> gmset bubble2 [gmapply \
	[gmcompose [gmtranslate .02 0 0] [gmdilate .01 .01 .01]] \
	$s]


Thu Aug 26 15:34:19 EDT 1999
-----> gmset cyl2 [gmcavity $cyl1 $bubble1]


Thu Aug 26 15:34:19 EDT 1999
-----> gmset cyl3 [gmcavity $cyl2 $bubble2]


Thu Aug 26 15:34:20 EDT 1999
-----> 
set show 0

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


Thu Aug 26 15:34:20 EDT 1999
-----> 
gmset mesh [gmmeshgen $cyl3 curve "(const 1.2)" show $show]


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

Maximum aspect ratio =        49.0823 achieved in
simplex #2922 of topological entity cylinder (3:0) which has vertices 107 80 137 739 
Maximum global side length =  1.03923
Minimum global altitude =     0.000131217
Number of nodes = 1179 number of elements = 6054
49.082293552010853
Thu Aug 26 15:57:07 EDT 1999
-----> 
if {[llength [info globals aspprod]]} {
    global aspprod
    set aspprod [expr $aspprod * $asp]
    gmset {numvtx numelt} [gmmeshsize $mesh]
    global meshsizesum
    set meshsizesum [expr $meshsizesum + $numvtx]
}

2394
Thu Aug 26 15:57:07 EDT 1999
-----> 
# ------------------------------------------------------------------


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


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


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


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


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


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


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


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


Thu Aug 26 15:57:07 EDT 1999