test9.di

算断裂的

Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> 
# Topological vertices bounding the crack


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> 
lappend vertlist vcrk0 {} {} {} {{vertex 72}}

vbottom0 {} {} {} {{vertex 36}} vbottom1 {} {} {} {{vertex 39}} vbottom2 {} {} {} {{vertex 42}} vbottom3 {} {} {} {{vertex 45}} vbottom4 {} {} {} {{vertex 48}} vbottom5 {} {} {} {{vertex 51}} vtop0 {} {} {} {{vertex 54}} vtop1 {} {} {} {{vertex 57}} vtop2 {} {} {} {{vertex 60}} vtop3 {} {} {} {{vertex 63}} vtop4 {} {} {} {{vertex 66}} vtop5 {} {} {} {{vertex 69}} vcrk0 {} {} {} {{vertex 72}}
Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> lappend vertlist vcrk1 {} {} {} {{vertex 73}}

vbottom0 {} {} {} {{vertex 36}} vbottom1 {} {} {} {{vertex 39}} vbottom2 {} {} {} {{vertex 42}} vbottom3 {} {} {} {{vertex 45}} vbottom4 {} {} {} {{vertex 48}} vbottom5 {} {} {} {{vertex 51}} vtop0 {} {} {} {{vertex 54}} vtop1 {} {} {} {{vertex 57}} vtop2 {} {} {} {{vertex 60}} vtop3 {} {} {} {{vertex 63}} vtop4 {} {} {} {{vertex 66}} vtop5 {} {} {} {{vertex 69}} vcrk0 {} {} {} {{vertex 72}} vcrk1 {} {} {} {{vertex 73}}
Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> lappend vertlist vcrk2 {} {} {} {{vertex 74}}

vbottom0 {} {} {} {{vertex 36}} vbottom1 {} {} {} {{vertex 39}} vbottom2 {} {} {} {{vertex 42}} vbottom3 {} {} {} {{vertex 45}} vbottom4 {} {} {} {{vertex 48}} vbottom5 {} {} {} {{vertex 51}} vtop0 {} {} {} {{vertex 54}} vtop1 {} {} {} {{vertex 57}} vtop2 {} {} {} {{vertex 60}} vtop3 {} {} {} {{vertex 63}} vtop4 {} {} {} {{vertex 66}} vtop5 {} {} {} {{vertex 69}} vcrk0 {} {} {} {{vertex 72}} vcrk1 {} {} {} {{vertex 73}} vcrk2 {} {} {} {{vertex 74}}
Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> 

# Create the list of topological edges -- 6 on top, 6 on the bottom,


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> # and 6 on the side to bound the two hexagons and side faces,


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> # plus two circular curves and three edges bounding the crack.


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> 
set edgelist {}


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> 
# edges for bottom hex


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> for {set j 0} {$j < 6} {incr j} {
    set nextj [expr ($j+1)%6]
    lappend edgelist ebottom$j {} [list vbottom$j vbottom$nextj] {} \
	    [list [list bezier_curve 1 [expr 36+3*$j] [expr 36+3*$nextj]]]
}


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> 
#edges for top hex


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> for {set j 0} {$j < 6} {incr j} {
    set nextj [expr ($j+1)%6]
    lappend edgelist etop$j {} [list vtop$j vtop$nextj] {} \
	    [list [list bezier_curve 1 [expr 54+3*$j] [expr 54+3*$nextj]]]
}


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> 
#edges connecting top and bottom hex


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> 
for {set j 0} {$j < 6} {incr j} {
    lappend edgelist eside$j {} [list vbottom$j vtop$j] {} \
	    [list [list bezier_curve 1 [expr 36+3*$j] [expr 54+3*$j]]]
}


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


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


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> set curves {}


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> for {set j 0} {$j < 6} {incr j} {
    set nextj [expr ($j+1)%6]
    lappend curves \
	    [list bezier_curve 3 [expr 3*$j] [expr 3*$j+1] \
	    [expr 3*$j+2] [expr 3*$nextj]]
}


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> 
lappend edgelist ebottomcirc {} {} {} $curves

ebottom0 {} {vbottom0 vbottom1} {} {{bezier_curve 1 36 39}} ebottom1 {} {vbottom1 vbottom2} {} {{bezier_curve 1 39 42}} ebottom2 {} {vbottom2 vbottom3} {} {{bezier_curve 1 42 45}} ebottom3 {} {vbottom3 vbottom4} {} {{bezier_curve 1 45 48}} ebottom4 {} {vbottom4 vbottom5} {} {{bezier_curve 1 48 51}} ebottom5 {} {vbottom5 vbottom0} {} {{bezier_curve 1 51 36}} etop0 {} {vtop0 vtop1} {} {{bezier_curve 1 54 57}} etop1 {} {vtop1 vtop2} {} {{bezier_curve 1 57 60}} etop2 {} {vtop2 vtop3} {} {{bezier_curve 1 60 63}} etop3 {} {vtop3 vtop4} {} {{bezier_curve 1 63 66}} etop4 {} {vtop4 vtop5} {} {{bezier_curve 1 66 69}} etop5 {} {vtop5 vtop0} {} {{bezier_curve 1 69 54}} eside0 {} {vbottom0 vtop0} {} {{bezier_curve 1 36 54}} eside1 {} {vbottom1 vtop1} {} {{bezier_curve 1 39 57}} eside2 {} {vbottom2 vtop2} {} {{bezier_curve 1 42 60}} eside3 {} {vbottom3 vtop3} {} {{bezier_curve 1 45 63}} eside4 {} {vbottom4 vtop4} {} {{bezier_curve 1 48 66}} eside5 {} {vbottom5 vtop5} {} {{bezier_curve 1 51 69}} ebottomcirc {} {} {} {{bezier_curve 3 0 1 2 3} {bezier_curve 3 3 4 5 6} {bezier_curve 3 6 7 8 9} {bezier_curve 3 9 10 11 12} {bezier_curve 3 12 13 14 15} {bezier_curve 3 15 16 17 0}}
Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> 
# top circle


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> set curves {}


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> for {set j 0} {$j < 6} {incr j} {
    set nextj [expr ($j+1)%6]
    lappend curves \
	    [list bezier_curve 3 [expr 3*$j+18] [expr 3*$j+19] \
	    [expr 3*$j+20] [expr 3*$nextj +18]]
}


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> lappend edgelist etopcirc {} {} {} $curves

ebottom0 {} {vbottom0 vbottom1} {} {{bezier_curve 1 36 39}} ebottom1 {} {vbottom1 vbottom2} {} {{bezier_curve 1 39 42}} ebottom2 {} {vbottom2 vbottom3} {} {{bezier_curve 1 42 45}} ebottom3 {} {vbottom3 vbottom4} {} {{bezier_curve 1 45 48}} ebottom4 {} {vbottom4 vbottom5} {} {{bezier_curve 1 48 51}} ebottom5 {} {vbottom5 vbottom0} {} {{bezier_curve 1 51 36}} etop0 {} {vtop0 vtop1} {} {{bezier_curve 1 54 57}} etop1 {} {vtop1 vtop2} {} {{bezier_curve 1 57 60}} etop2 {} {vtop2 vtop3} {} {{bezier_curve 1 60 63}} etop3 {} {vtop3 vtop4} {} {{bezier_curve 1 63 66}} etop4 {} {vtop4 vtop5} {} {{bezier_curve 1 66 69}} etop5 {} {vtop5 vtop0} {} {{bezier_curve 1 69 54}} eside0 {} {vbottom0 vtop0} {} {{bezier_curve 1 36 54}} eside1 {} {vbottom1 vtop1} {} {{bezier_curve 1 39 57}} eside2 {} {vbottom2 vtop2} {} {{bezier_curve 1 42 60}} eside3 {} {vbottom3 vtop3} {} {{bezier_curve 1 45 63}} eside4 {} {vbottom4 vtop4} {} {{bezier_curve 1 48 66}} eside5 {} {vbottom5 vtop5} {} {{bezier_curve 1 51 69}} ebottomcirc {} {} {} {{bezier_curve 3 0 1 2 3} {bezier_curve 3 3 4 5 6} {bezier_curve 3 6 7 8 9} {bezier_curve 3 9 10 11 12} {bezier_curve 3 12 13 14 15} {bezier_curve 3 15 16 17 0}} etopcirc {} {} {} {{bezier_curve 3 18 19 20 21} {bezier_curve 3 21 22 23 24} {bezier_curve 3 24 25 26 27} {bezier_curve 3 27 28 29 30} {bezier_curve 3 30 31 32 33} {bezier_curve 3 33 34 35 18}}
Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> 
# crack


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> lappend edgelist ecrk0 {} {vcrk0 vcrk1} {} {{bezier_curve 1 72 73}}

ebottom0 {} {vbottom0 vbottom1} {} {{bezier_curve 1 36 39}} ebottom1 {} {vbottom1 vbottom2} {} {{bezier_curve 1 39 42}} ebottom2 {} {vbottom2 vbottom3} {} {{bezier_curve 1 42 45}} ebottom3 {} {vbottom3 vbottom4} {} {{bezier_curve 1 45 48}} ebottom4 {} {vbottom4 vbottom5} {} {{bezier_curve 1 48 51}} ebottom5 {} {vbottom5 vbottom0} {} {{bezier_curve 1 51 36}} etop0 {} {vtop0 vtop1} {} {{bezier_curve 1 54 57}} etop1 {} {vtop1 vtop2} {} {{bezier_curve 1 57 60}} etop2 {} {vtop2 vtop3} {} {{bezier_curve 1 60 63}} etop3 {} {vtop3 vtop4} {} {{bezier_curve 1 63 66}} etop4 {} {vtop4 vtop5} {} {{bezier_curve 1 66 69}} etop5 {} {vtop5 vtop0} {} {{bezier_curve 1 69 54}} eside0 {} {vbottom0 vtop0} {} {{bezier_curve 1 36 54}} eside1 {} {vbottom1 vtop1} {} {{bezier_curve 1 39 57}} eside2 {} {vbottom2 vtop2} {} {{bezier_curve 1 42 60}} eside3 {} {vbottom3 vtop3} {} {{bezier_curve 1 45 63}} eside4 {} {vbottom4 vtop4} {} {{bezier_curve 1 48 66}} eside5 {} {vbottom5 vtop5} {} {{bezier_curve 1 51 69}} ebottomcirc {} {} {} {{bezier_curve 3 0 1 2 3} {bezier_curve 3 3 4 5 6} {bezier_curve 3 6 7 8 9} {bezier_curve 3 9 10 11 12} {bezier_curve 3 12 13 14 15} {bezier_curve 3 15 16 17 0}} etopcirc {} {} {} {{bezier_curve 3 18 19 20 21} {bezier_curve 3 21 22 23 24} {bezier_curve 3 24 25 26 27} {bezier_curve 3 27 28 29 30} {bezier_curve 3 30 31 32 33} {bezier_curve 3 33 34 35 18}} ecrk0 {} {vcrk0 vcrk1} {} {{bezier_curve 1 72 73}}
Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> lappend edgelist ecrk1 {} {vcrk1 vcrk2} {} {{bezier_curve 1 73 74}}

ebottom0 {} {vbottom0 vbottom1} {} {{bezier_curve 1 36 39}} ebottom1 {} {vbottom1 vbottom2} {} {{bezier_curve 1 39 42}} ebottom2 {} {vbottom2 vbottom3} {} {{bezier_curve 1 42 45}} ebottom3 {} {vbottom3 vbottom4} {} {{bezier_curve 1 45 48}} ebottom4 {} {vbottom4 vbottom5} {} {{bezier_curve 1 48 51}} ebottom5 {} {vbottom5 vbottom0} {} {{bezier_curve 1 51 36}} etop0 {} {vtop0 vtop1} {} {{bezier_curve 1 54 57}} etop1 {} {vtop1 vtop2} {} {{bezier_curve 1 57 60}} etop2 {} {vtop2 vtop3} {} {{bezier_curve 1 60 63}} etop3 {} {vtop3 vtop4} {} {{bezier_curve 1 63 66}} etop4 {} {vtop4 vtop5} {} {{bezier_curve 1 66 69}} etop5 {} {vtop5 vtop0} {} {{bezier_curve 1 69 54}} eside0 {} {vbottom0 vtop0} {} {{bezier_curve 1 36 54}} eside1 {} {vbottom1 vtop1} {} {{bezier_curve 1 39 57}} eside2 {} {vbottom2 vtop2} {} {{bezier_curve 1 42 60}} eside3 {} {vbottom3 vtop3} {} {{bezier_curve 1 45 63}} eside4 {} {vbottom4 vtop4} {} {{bezier_curve 1 48 66}} eside5 {} {vbottom5 vtop5} {} {{bezier_curve 1 51 69}} ebottomcirc {} {} {} {{bezier_curve 3 0 1 2 3} {bezier_curve 3 3 4 5 6} {bezier_curve 3 6 7 8 9} {bezier_curve 3 9 10 11 12} {bezier_curve 3 12 13 14 15} {bezier_curve 3 15 16 17 0}} etopcirc {} {} {} {{bezier_curve 3 18 19 20 21} {bezier_curve 3 21 22 23 24} {bezier_curve 3 24 25 26 27} {bezier_curve 3 27 28 29 30} {bezier_curve 3 30 31 32 33} {bezier_curve 3 33 34 35 18}} ecrk0 {} {vcrk0 vcrk1} {} {{bezier_curve 1 72 73}} ecrk1 {} {vcrk1 vcrk2} {} {{bezier_curve 1 73 74}}
Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> lappend edgelist ecrk2 {} {vcrk2 vcrk0} {} {{bezier_curve 1 74 72}}

ebottom0 {} {vbottom0 vbottom1} {} {{bezier_curve 1 36 39}} ebottom1 {} {vbottom1 vbottom2} {} {{bezier_curve 1 39 42}} ebottom2 {} {vbottom2 vbottom3} {} {{bezier_curve 1 42 45}} ebottom3 {} {vbottom3 vbottom4} {} {{bezier_curve 1 45 48}} ebottom4 {} {vbottom4 vbottom5} {} {{bezier_curve 1 48 51}} ebottom5 {} {vbottom5 vbottom0} {} {{bezier_curve 1 51 36}} etop0 {} {vtop0 vtop1} {} {{bezier_curve 1 54 57}} etop1 {} {vtop1 vtop2} {} {{bezier_curve 1 57 60}} etop2 {} {vtop2 vtop3} {} {{bezier_curve 1 60 63}} etop3 {} {vtop3 vtop4} {} {{bezier_curve 1 63 66}} etop4 {} {vtop4 vtop5} {} {{bezier_curve 1 66 69}} etop5 {} {vtop5 vtop0} {} {{bezier_curve 1 69 54}} eside0 {} {vbottom0 vtop0} {} {{bezier_curve 1 36 54}} eside1 {} {vbottom1 vtop1} {} {{bezier_curve 1 39 57}} eside2 {} {vbottom2 vtop2} {} {{bezier_curve 1 42 60}} eside3 {} {vbottom3 vtop3} {} {{bezier_curve 1 45 63}} eside4 {} {vbottom4 vtop4} {} {{bezier_curve 1 48 66}} eside5 {} {vbottom5 vtop5} {} {{bezier_curve 1 51 69}} ebottomcirc {} {} {} {{bezier_curve 3 0 1 2 3} {bezier_curve 3 3 4 5 6} {bezier_curve 3 6 7 8 9} {bezier_curve 3 9 10 11 12} {bezier_curve 3 12 13 14 15} {bezier_curve 3 15 16 17 0}} etopcirc {} {} {} {{bezier_curve 3 18 19 20 21} {bezier_curve 3 21 22 23 24} {bezier_curve 3 24 25 26 27} {bezier_curve 3 27 28 29 30} {bezier_curve 3 30 31 32 33} {bezier_curve 3 33 34 35 18}} ecrk0 {} {vcrk0 vcrk1} {} {{bezier_curve 1 72 73}} ecrk1 {} {vcrk1 vcrk2} {} {{bezier_curve 1 73 74}} ecrk2 {} {vcrk2 vcrk0} {} {{bezier_curve 1 74 72}}
Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> 
# Make the topological surfaces.  There is one on bottom, one on top,


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> # one for the inside, six for the outside, and one for the crack.


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> 
set surflist {}


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


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> set patches {}


Thu Sep 02 18:11:00 Eastern Daylight Time 1999
-----> for {set j 0} {$j < 6} {incr j} {
    set nextj [expr ($j+1)%6]
    lappend patches [list bezier_quad 3 1 [expr 36+3*$j] \
	    [expr 37+3*$j] [expr 38+3*$j] [expr 36+3*$nextj] \
	    [expr 3*$j] [expr 3*$j+1] [expr 3*$j+2] [expr 3*$nextj]]
}