test7.di

来自「算断裂的」· DI 代码 · 共 714 行 · 第 1/3 页

Thu Sep 02 18:10:55 Eastern Daylight Time 1999
-----> # QMG test 7 - a circle with two elliptic holes and a crack.


Thu Sep 02 18:10:55 Eastern Daylight Time 1999
-----> # Make an approximation to a circle using 6 Bezier cubic arcs.


Thu Sep 02 18:10:55 Eastern Daylight Time 1999
-----> 
set numseg 6

6
Thu Sep 02 18:10:55 Eastern Daylight Time 1999
-----> global PI


Thu Sep 02 18:10:55 Eastern Daylight Time 1999
-----> global interactive


Thu Sep 02 18:10:55 Eastern Daylight Time 1999
-----> 
gmset {verts codes} [gm_circ_approx 0 [expr 2*$PI] $numseg]


Thu Sep 02 18:10:55 Eastern Daylight Time 1999
-----> gmset verts [lrange $verts 0 [expr 3*$numseg-1]]


Thu Sep 02 18:10:55 Eastern Daylight Time 1999
-----> set codes {}


Thu Sep 02 18:10:55 Eastern Daylight Time 1999
-----> for {set j 0} {$j < $numseg} {incr j} {
    if {$j == 0} {
	lappend codes 0
    } else {
	lappend codes 2
    } 
    lappend codes 3 3
}


Thu Sep 02 18:10:55 Eastern Daylight Time 1999
-----> 
set c0 [gm_cpoly $verts $codes]

brep_v2.0 2 2
()
(
1 0 # cp 0
1.0000000000000002 0.3568220897730896 # cp 1
0.80901699437494801 0.68761435889789313 # cp 2
0.50000000000000089 0.86602540378443815 # cp 3
0.19098300562505391 1.0444364486709834 # cp 4
-0.19098300562505077 1.0444364486709841 # cp 5
-0.49999999999999828 0.8660254037844396 # cp 6
-0.80901699437494579 0.68761435889789557 # cp 7
-0.99999999999999911 0.35682208977309243 # cp 8
-1 2.7869958944826134e-015 # cp 9
-1.0000000000000011 -0.35682208977308688 # cp 10
-0.80901699437495012 -0.68761435889789069 # cp 11
-0.50000000000000344 -0.8660254037844366 # cp 12
-0.19098300562505705 -1.0444364486709827 # cp 13
0.19098300562504739 -1.0444364486709845 # cp 14
0.49999999999999545 -0.86602540378444126 # cp 15
0.80901699437494368 -0.68761435889789813 # cp 16
0.99999999999999778 -0.35682208977309615 # cp 17
) 
(   # begin topological vertices
v0  ()  () ()  ( (vertex 0 ) )
 ) # end topological vertices
(   # begin topological edges
e18  ()  (v0 v0 ) ()  ( (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 ) )
 ) # end topological edges
(   # begin topological regions
mregion  ()  (e18 ) ()  ( )
 ) # end topological regions

Thu Sep 02 18:10:55 Eastern Daylight Time 1999
-----> 
## Make two elliptic holes.


Thu Sep 02 18:10:55 Eastern Daylight Time 1999
-----> 
set hole1 [gmapply {{.2 0 0} {0 .1 -.5}} $c0]

brep_v2.0 2 2
()
(
0.20000000000000001 0.5 # cp 0
0.20000000000000007 0.535682208977309 # cp 1
0.16180339887498962 0.56876143588978934 # cp 2
0.10000000000000019 0.58660254037844384 # cp 3
0.038196601125010783 0.60444364486709834 # cp 4
-0.038196601125010159 0.60444364486709845 # cp 5
-0.099999999999999659 0.58660254037844395 # cp 6
-0.16180339887498918 0.56876143588978956 # cp 7
-0.19999999999999984 0.53568220897730923 # cp 8
-0.20000000000000001 0.50000000000000033 # cp 9
-0.20000000000000023 0.46431779102269133 # cp 10
-0.16180339887499004 0.43123856411021094 # cp 11
-0.1000000000000007 0.41339745962155633 # cp 12
-0.038196601125011415 0.39555635513290172 # cp 13
0.038196601125009479 0.39555635513290155 # cp 14
0.09999999999999909 0.41339745962155588 # cp 15
0.16180339887498874 0.43123856411021017 # cp 16
0.19999999999999957 0.46431779102269038 # cp 17
) 
(   # begin topological vertices
v0  ()  () ()  ( (vertex 0 ) )
 ) # end topological vertices
(   # begin topological edges
e18  ()  (v0 v0 ) ()  ( (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 ) )
 ) # end topological edges
(   # begin topological regions
mregion  ()  (e18 ) ()  ( )
 ) # end topological regions

Thu Sep 02 18:10:55 Eastern Daylight Time 1999
-----> set c [gmcavity $c0 $hole1]

brep_v2.0 2 2
()
(
1 0 # cp 0
1.0000000000000002 0.3568220897730896 # cp 1
0.80901699437494801 0.68761435889789313 # cp 2
0.50000000000000089 0.86602540378443815 # cp 3
0.19098300562505391 1.0444364486709834 # cp 4
-0.19098300562505077 1.0444364486709841 # cp 5
-0.49999999999999828 0.8660254037844396 # cp 6
-0.80901699437494579 0.68761435889789557 # cp 7
-0.99999999999999911 0.35682208977309243 # cp 8
-1 2.7869958944826134e-015 # cp 9
-1.0000000000000011 -0.35682208977308688 # cp 10
-0.80901699437495012 -0.68761435889789069 # cp 11
-0.50000000000000344 -0.8660254037844366 # cp 12
-0.19098300562505705 -1.0444364486709827 # cp 13
0.19098300562504739 -1.0444364486709845 # cp 14
0.49999999999999545 -0.86602540378444126 # cp 15
0.80901699437494368 -0.68761435889789813 # cp 16
0.99999999999999778 -0.35682208977309615 # cp 17
0.20000000000000001 0.5 # cp 18
0.20000000000000007 0.535682208977309 # cp 19
0.16180339887498962 0.56876143588978934 # cp 20
0.10000000000000019 0.58660254037844384 # cp 21
0.038196601125010783 0.60444364486709834 # cp 22
-0.038196601125010159 0.60444364486709845 # cp 23
-0.099999999999999659 0.58660254037844395 # cp 24
-0.16180339887498918 0.56876143588978956 # cp 25
-0.19999999999999984 0.53568220897730923 # cp 26
-0.20000000000000001 0.50000000000000033 # cp 27
-0.20000000000000023 0.46431779102269133 # cp 28
-0.16180339887499004 0.43123856411021094 # cp 29
-0.1000000000000007 0.41339745962155633 # cp 30
-0.038196601125011415 0.39555635513290172 # cp 31
0.038196601125009479 0.39555635513290155 # cp 32
0.09999999999999909 0.41339745962155588 # cp 33
0.16180339887498874 0.43123856411021017 # cp 34
0.19999999999999957 0.46431779102269038 # cp 35
) 
(   # begin topological vertices
v0  ()  () ()  ( (vertex 0 ) )
v0_cav1  ()  () ()  ( (vertex 18 ) )
 ) # end topological vertices
(   # begin topological edges
e18  ()  (v0 v0 ) ()  ( (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 ) )
e18_cav1  ()  (v0_cav1 v0_cav1 ) ()  ( (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 ) )
 ) # end topological edges
(   # begin topological regions
mregion  ()  (e18 e18_cav1 ) ()  ( )
 ) # end topological regions

Thu Sep 02 18:10:55 Eastern Daylight Time 1999
-----> set hole2 [gmapply {{.2 0 0} {0 .1 .5}} $c0]

brep_v2.0 2 2
()
(
0.20000000000000001 -0.5 # cp 0
0.20000000000000007 -0.46431779102269105 # cp 1
0.16180339887498962 -0.43123856411021066 # cp 2
0.10000000000000019 -0.41339745962155616 # cp 3
0.038196601125010783 -0.39555635513290166 # cp 4
-0.038196601125010159 -0.3955563551329016 # cp 5
-0.099999999999999659 -0.41339745962155605 # cp 6
-0.16180339887498918 -0.43123856411021044 # cp 7
-0.19999999999999984 -0.46431779102269077 # cp 8
-0.20000000000000001 -0.49999999999999972 # cp 9
-0.20000000000000023 -0.53568220897730867 # cp 10
-0.16180339887499004 -0.56876143588978911 # cp 11
-0.1000000000000007 -0.58660254037844362 # cp 12
-0.038196601125011415 -0.60444364486709823 # cp 13
0.038196601125009479 -0.60444364486709845 # cp 14
0.09999999999999909 -0.58660254037844417 # cp 15
0.16180339887498874 -0.56876143588978978 # cp 16
0.19999999999999957 -0.53568220897730967 # cp 17
) 
(   # begin topological vertices
v0  ()  () ()  ( (vertex 0 ) )
 ) # end topological vertices
(   # begin topological edges
e18  ()  (v0 v0 ) ()  ( (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 ) )
 ) # end topological edges
(   # begin topological regions
mregion  ()  (e18 ) ()  ( )
 ) # end topological regions

Thu Sep 02 18:10:55 Eastern Daylight Time 1999