18e25fc9 |
1 | if { [array get Draw_Groups "TOPOLOGY Check commands"] == "" } { |
2 | pload TOPTEST |
3 | pload AISV |
4 | } |
5 | # To prevent loops limit to 10 minutes |
6 | cpulimit 600 |
7 | |
8 | if { [info exists imagedir] == 0 } { |
9 | set imagedir . |
10 | } |
11 | |
12 | if { [info exists test_image ] == 0 } { |
13 | set test_image photo |
14 | } |
15 | |
16 | |
17 | # |
18 | #evaluate n points beetween u1 and u2 on the curve 2d c |
19 | # |
20 | |
21 | proc val2d { c u1 u2 n } { |
22 | |
23 | dset du ($u2-$u1)/$n |
18e25fc9 |
24 | |
25 | set i 1 |
26 | |
27 | for {dset u $u1} { [dval u] <= $u2} {dset u ($u1+$i*[dval du])} { |
1951a27c |
28 | 2dcvalue $c u x y dx dy d2x d2y ; |
18e25fc9 |
29 | global p_$i d1_$i d2_$i |
30 | point p_$i x y; |
31 | puts "u = [dval u]" |
32 | puts "p_$i [dval x ] [dval y]" |
33 | puts "d1_$i [dval dx ] [dval dy]"; |
34 | puts "d2_$i [dval d2x] [dval d2y]"; |
35 | |
36 | copy p_$i . ; |
37 | point d1_$i dx dy |
38 | point d2_$i d2x d2y |
39 | incr i 1 |
40 | } |
41 | } |
42 | |
43 | # |
44 | #evaluate n points beetween u1 and u2 on the curve 3d c |
45 | # |
46 | proc val3d { c u1 u2 n } { |
47 | |
48 | dset du ($u2-$u1)/$n |
18e25fc9 |
49 | |
50 | set i 1 |
51 | |
52 | for {dset u $u1} { [dval u] <= $u2} {dset u (u+[dval du])} { |
1951a27c |
53 | cvalue $c u x y z dx dy dz d2x d2y d2z ; |
18e25fc9 |
54 | point p_$i x y z; |
55 | puts "u = [dval u]" |
56 | puts "p_$i [dval x ] [dval y ] [dval z]"; |
57 | puts "d1_$i [dval dx ] [dval dy ] [dval dz]"; |
58 | puts "d2_$i [dval d2x] [dval d2y] [dval d2z]"; |
59 | copy p_$i . ; |
60 | point d1_$i dx dy dz |
61 | point d2_$i d2x d2y d2z |
62 | incr i |
63 | } |
64 | } |
65 | |
66 | proc compare {r1 r2 tol} { |
67 | if {$r1 - $r2 >= $tol} { |
68 | puts "Error : evalution" |
69 | } |
70 | if {$r2 - $r1 >= $tol} { |
71 | puts "Error : evalution" |
72 | } |
73 | } |
74 | |
75 | proc comparepnt2d {p1 p2 tol} { |
1951a27c |
76 | coord $p1 x1 y1 |
77 | coord $p2 x2 y2 |
18e25fc9 |
78 | compare [dval x1] [dval x2] $tol |
79 | compare [dval y1] [dval y2] $tol |
80 | } |
81 | |
82 | proc comparepnt3d {p1 p2 tol} { |
83 | upvar $p1 pp1 |
84 | upvar $p2 pp2 |
85 | coord pp1 x1 y1 z1 |
86 | coord pp2 x2 y2 z2 |
87 | compare [dval x1] [dval x2] $tol |
88 | compare [dval y1] [dval y2] $tol |
89 | compare [dval z1] [dval z2] $tol |
90 | } |
91 | |