test10.di

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

135 139 170 
167 135 170 
32 34 171 
34 33 171 
168 167 171 
170 32 171 
167 170 171 
31 2 172 
33 31 172 
2 169 172 
169 168 172 
171 33 172 
168 171 172 
) # 2:0
) # end mesh info for topological entities of dimension 2

Thu Sep 02 18:13:12 Eastern Daylight Time 1999
-----> if {$show} {
    gmviz $m1
}


Thu Sep 02 18:13:12 Eastern Daylight Time 1999
-----> 
set p2 [gmcoarsetopo $p .3]

brep_v2.0 2 2
()
(
1 0 # cp 0
0.97814760073380569 0.20791169081775912 # cp 1
0.91354545764260109 0.40673664307579982 # cp 2
0.80901699437494778 0.58778525229247269 # cp 3
0.66913060635885879 0.74314482547739369 # cp 4
0.50000000000000089 0.86602540378443815 # cp 5
0.30901699437494851 0.9510565162951532 # cp 6
0.104528463267655 0.99452189536827318 # cp 7
-0.10452846326765179 0.99452189536827351 # cp 8
-0.30901699437494545 0.9510565162951542 # cp 9
-0.49999999999999828 0.8660254037844396 # cp 10
-0.66913060635885657 0.74314482547739569 # cp 11
-0.80901699437494601 0.58778525229247502 # cp 12
-0.91354545764259976 0.40673664307580287 # cp 13
-0.97814760073380502 0.20791169081776234 # cp 14
-1 3.231085104332676e-015 # cp 15
-0.97814760073380635 -0.20791169081775601 # cp 16
-0.91354545764260231 -0.40673664307579699 # cp 17
-0.80901699437494989 -0.58778525229246981 # cp 18
-0.66913060635886146 -0.74314482547739136 # cp 19
-0.50000000000000344 -0.8660254037844366 # cp 20
-0.30901699437495184 -0.9510565162951522 # cp 21
-0.10452846326765777 -0.99452189536827285 # cp 22
0.10452846326764857 -0.99452189536827384 # cp 23
0.30901699437494301 -0.95105651629515497 # cp 24
0.49999999999999545 -0.86602540378444126 # cp 25
0.66913060635885391 -0.74314482547739813 # cp 26
0.80901699437494368 -0.58778525229247836 # cp 27
0.91354545764259842 -0.40673664307580581 # cp 28
0.97814760073380425 -0.20791169081776595 # cp 29
) 
(   # begin topological vertices
 ) # end topological vertices
(   # begin topological edges

    Merge_of_$_e1_$_e2_$_e3_$_e4_$_e5_$_e6_$_e7_$_e8_$_e9_$_e10_$_e11_$_e12_$_e13_$_e14_$_e15_$_e16_$_e17_$_e18_$_e19_$_e20_$_e21_$_e22_$_e23_$_e24_$_e25_$_e26_$_e27_$_e28_$_e29_$_e30
      ()  () ()  ( (bezier_curve 1 0 1 )
    (bezier_curve 1 1 2 )
    (bezier_curve 1 2 3 )
    (bezier_curve 1 3 4 )
    (bezier_curve 1 4 5 )
    (bezier_curve 1 5 6 )
    (bezier_curve 1 6 7 )
    (bezier_curve 1 7 8 )
    (bezier_curve 1 8 9 )
    (bezier_curve 1 9 10 )
    (bezier_curve 1 10 11 )
    (bezier_curve 1 11 12 )
    (bezier_curve 1 12 13 )
    (bezier_curve 1 13 14 )
    (bezier_curve 1 14 15 )
    (bezier_curve 1 15 16 )
    (bezier_curve 1 16 17 )
    (bezier_curve 1 17 18 )
    (bezier_curve 1 18 19 )
    (bezier_curve 1 19 20 )
    (bezier_curve 1 20 21 )
    (bezier_curve 1 21 22 )
    (bezier_curve 1 22 23 )
    (bezier_curve 1 23 24 )
    (bezier_curve 1 24 25 )
    (bezier_curve 1 25 26 )
    (bezier_curve 1 26 27 )
    (bezier_curve 1 27 28 )
    (bezier_curve 1 28 29 )
    (bezier_curve 1 29 0 ) )
 ) # end topological edges
(   # begin topological regions
mregion  ()  (Merge_of_$_e1_$_e2_$_e3_$_e4_$_e5_$_e6_$_e7_$_e8_$_e9_$_e10_$_e11_$_e12_$_e13_$_e14_$_e15_$_e16_$_e17_$_e18_$_e19_$_e20_$_e21_$_e22_$_e23_$_e24_$_e25_$_e26_$_e27_$_e28_$_e29_$_e30
     ) ()  ( )
 ) # end topological regions

Thu Sep 02 18:13:12 Eastern Daylight Time 1999
-----> set m2 [gmmeshgen $p2 show $show]

mesh_v2.01 2 2
()
( # begin list of mesh nodes
0 -0.99942422800197739 0.005478104631726266 
1 0.95456934360066015 0.28047810463172618 
2 0.95812923178640563 -0.26952189536827376 
3 0.10000000000000009 0.99452189536827307 
4 0.10000000000000009 -0.99452189536827384 
5 -0.44999999999999996 -0.8882868380498653 
6 -0.72499999999999998 -0.68109557770362383 
7 -0.17499999999999993 -0.97954270782484321 
8 0.65000000000000013 -0.75704402458363507 
9 0.375 -0.92167898944800553 
10 -0.44999999999999996 0.88828683804986552 
11 -0.72499999999999998 0.68109557770362361 
12 -0.17499999999999993 0.97954270782484287 
13 0.65000000000000013 0.75704402458363429 
14 0.375 0.92167898944800575 
15 0.99942422800197706 0.005478104631726266 
16 -0.83399510514118114 -0.54452189536827378 
17 -0.5640065882224885 -0.8195218953682738 
18 -0.95812923178640519 -0.26952189536827376 
19 0.83399510514118047 -0.54452189536827378 
20 0.56400658822248817 -0.8195218953682738 
21 -0.82766953477362848 0.55547810463172631 
22 -0.95456934360066026 0.28047810463172618 
23 -0.54892665987537315 0.83047810463172644 
24 0.8276695347736287 0.55547810463172631 
25 0.54892665987537348 0.83047810463172644 
26 0.10000000000000009 0.005478104631726266 
27 -0.72499999999999998 -0.54452189536827378 
28 -0.72499999999999998 -0.26952189536827376 
29 -0.72499999999999998 0.005478104631726266 
30 -0.72499999999999998 0.28047810463172618 
31 -0.72499999999999998 0.55547810463172631 
32 -0.44999999999999996 -0.8195218953682738 
33 -0.44999999999999996 -0.54452189536827378 
34 -0.44999999999999996 -0.26952189536827376 
35 -0.44999999999999996 0.005478104631726266 
36 -0.44999999999999996 0.28047810463172618 
37 -0.44999999999999996 0.55547810463172631 
38 -0.44999999999999996 0.83047810463172644 
39 -0.17499999999999993 -0.8195218953682738 
40 -0.17499999999999993 -0.54452189536827378 
41 -0.17499999999999993 0.55547810463172631 
42 -0.17499999999999993 0.83047810463172644 
43 0.10000000000000009 -0.8195218953682738 
44 0.10000000000000009 -0.54452189536827378 
45 0.10000000000000009 0.55547810463172631 
46 0.10000000000000009 0.83047810463172644 
47 0.375 -0.8195218953682738 
48 0.375 -0.54452189536827378 
49 0.375 0.55547810463172631 
50 0.375 0.83047810463172644 
51 0.65000000000000013 -0.54452189536827378 
52 0.65000000000000013 -0.26952189536827376 
53 0.65000000000000013 0.005478104631726266 
54 0.65000000000000013 0.28047810463172618 
55 0.65000000000000013 0.55547810463172631 
56 0.92500000000000027 0.005478104631726266 
) # end list of mesh nodes

( # begin mesh info for topological entities of dimension 0
) # end mesh info for topological entities of dimension 0

( # begin mesh info for topological entities of dimension 1
( # mesh vertices lying on top. entity 1:0
0 14 0.97365177200869979 
1 1 0.36497639249920849 
2 28 0.69012840767315475 
3 7 0.021661388324725583 
4 22 0.97833861167530523 
5 20 0.26180339887500909 
6 18 0.60060879093745034 
7 21 0.65537658102025231 
8 25 0.88688856044030984 
9 24 0.34549150281254964 
10 9 0.73819660112501961 
11 11 0.39939120906258174 
12 8 0.34462341897977761 
13 4 0.11311143955971946 
14 5 0.65450849718747905 
15 0 0.026348227991315747 
16 17 0.761040103476189 
17 19 0.6215552607511079 
18 16 0.30987159232689354 
19 27 0.23895989652385857 
20 25 0.37844473924893651 
21 12 0.17844460557045616 
22 13 0.63502360750080733 
23 10 0.28928329962681715 
24 2 0.82155539442955738 
25 4 0.7107167003731959 
) ( # mesh entities lying on top. entity 1:0
16 6 
0 18 
18 16 
22 0 
21 22 
11 21 
17 5 
6 17 
10 23 
23 11 
7 4 
5 7 
3 12 
12 10 
4 9 
20 8 
9 20 
14 3 
13 25 
25 14 
19 2 
8 19 
1 24 
24 13 
15 1 
2 15 
) # 1:0
) # end mesh info for topological entities of dimension 1

( # begin mesh info for topological entities of dimension 2
() ( # mesh entities lying on top. entity 2:0
16 6 27 
0 18 28 
18 16 28 
16 27 28 
0 28 29 
22 0 30 
21 22 30 
0 29 30 
11 21 31 
21 30 31 
17 5 32 
6 17 32 
27 6 33 
6 32 33 
28 27 34 
33 26 34 
27 33 34 
29 28 35 
34 26 35 
28 34 35 
30 29 36 
35 26 36 
29 35 36 
11 31 37 
31 30 37 
36 26 37 
30 36 37 
10 23 38 
23 11 38 
11 37 38 
7 4 39 
5 7 39 
32 5 39 
26 33 40 
33 32 40 
32 39 40 
37 26 41 
3 12 42 
12 10 42 
10 38 42 
38 37 42 
37 41 42 
39 4 43 
26 40 44 
40 39 44 
39 43 44 
41 26 45 
3 42 46 
42 41 46 
41 45 46 
4 9 47 
20 8 47 
9 20 47 
43 4 47 
26 44 48 
44 43 48 
47 8 48 
43 47 48 
45 26 49 
14 3 50 
13 25 50 
25 14 50 
3 46 50 
46 45 50 
49 13 50 
45 49 50 
19 2 51 
8 19 51 
48 8 51 
26 48 51 
51 2 52 
26 51 52 
52 2 53 
26 52 53 
53 1 54 
26 53 54 
1 24 55 
24 13 55 
13 49 55 
49 26 55 
54 1 55 
26 54 55 
15 1 56 
2 15 56 
1 53 56 
53 2 56 
) # 2:0
) # end mesh info for topological entities of dimension 2

Thu Sep 02 18:13:12 Eastern Daylight Time 1999
-----> if {$show} {
    gmviz $m2
}


Thu Sep 02 18:13:12 Eastern Daylight Time 1999
-----> 
# How much does brep 2 deviate from G^1? 


Thu Sep 02 18:13:12 Eastern Daylight Time 1999
-----> 
set g1dev [gmchecknormals $p2]

0.20905692653531011 {Merge_of_$_e1_$_e2_$_e3_$_e4_$_e5_$_e6_$_e7_$_e8_$_e9_$_e10_$_e11_$_e12_$_e13_$_e14_$_e15_$_e16_$_e17_$_e18_$_e19_$_e20_$_e21_$_e22_$_e23_$_e24_$_e25_$_e26_$_e27_$_e28_$_e29_$_e30} 0
Thu Sep 02 18:13:12 Eastern Daylight Time 1999
-----> 
if {$show} {
    puts "deviation of brep2 from G^1 condition is $g1dev"
}


Thu Sep 02 18:13:13 Eastern Daylight Time 1999
-----> 

set asp [gmchecktri $p $m1]

Maximum aspect ratio =        13.4882 achieved in
simplex #272 of topological entity mregion (2:0) which has vertices 31 2 172 
Maximum global side length =  0.777817
Minimum global altitude =     0.00434737
Number of nodes = 173 number of elements = 278
13.48821015390233
Thu Sep 02 18:13:13 Eastern Daylight Time 1999
-----> 
gmset {numvtx numelt} [gmmeshsize $m1]


Thu Sep 02 18:13:13 Eastern Daylight Time 1999
-----> 
if {[llength [info globals meshsizesum]]} {
    global meshsizesum
    set meshsizesum [expr $meshsizesum + $numvtx]
    global aspprod
    set aspprod [expr $aspprod * $asp]
}

970611407963963.37
Thu Sep 02 18:13:13 Eastern Daylight Time 1999
-----> 
set asp [gmchecktri $p2 $m2]

Maximum aspect ratio =        6.73472 achieved in
simplex #51 of topological entity mregion (2:0) which has vertices 20 8 47 
Maximum global side length =  0.777817
Minimum global altitude =     0.0418737
Number of nodes = 57 number of elements = 86
6.7347199698642823
Thu Sep 02 18:13:13 Eastern Daylight Time 1999
-----> gmset {numvtx numelt} [gmmeshsize $m2]


Thu Sep 02 18:13:13 Eastern Daylight Time 1999
-----> 
if {[llength [info globals meshsizesum]]} {
    set meshsizesum [expr $meshsizesum + $numvtx]
    set aspprod [expr $aspprod * $asp]
}

6536796032192992.0
Thu Sep 02 18:13:13 Eastern Daylight Time 1999
-----> 


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


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


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


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


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


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


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


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


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


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