File: bug25292_35

package info (click to toggle)
oce 0.18.3-2
  • links: PTS, VCS
  • area: main
  • in suites: bookworm
  • size: 301,548 kB
  • sloc: cpp: 1,190,609; ansic: 67,225; sh: 11,630; tcl: 7,954; cs: 5,221; python: 2,867; java: 1,522; makefile: 342; xml: 292; perl: 37
file content (93 lines) | stat: -rw-r--r-- 2,293 bytes parent folder | download | duplicates (3)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
puts "================"
puts "OCC25292"
puts "================"
puts ""
#######################################################################
# Face/Face intersection algorithm gives different results for different order of the arguments
#######################################################################

puts "##############################"
puts "#!!!Search \"Attention\" keyword on this web-page for additional checking!!!"
puts "##############################"
puts ""
puts ""

# intersect command for trimmed surfaces

#This value must be equal to the analogical value in bug25292_35 and bug25292_36 of "bugs modalg_5" testgrid. 
set GoodNbCurv 2

restore [locate_data_file bug25292_f1.brep] f1
restore [locate_data_file bug25292_f2.brep] f2

mksurface s1 f1
mksurface s2 f2

# trim s1 in accordance with the UV-bounds of f1:
trimv s1t s1 19.1600000005 19.6600000005

# trim s2 in accordance with the UV-bounds of f2:
trim s2t s2 0. 1.570796326795 -275 275

###################
intersect res s1t s2t
###################

set che [whatis res]
set ind [string first "3d curve" $che]
if {${ind} >= 0} {
  #Only variable "res" exists
  
  copy res res_1
}

set ic 1
set AllowRepeate 1
while { $AllowRepeate != 0 } {
  set che [whatis res_$ic]
  set ind [string first "3d curve" $che]
  if {${ind} < 0} {
    set AllowRepeate 0
  } else {
    set le [length res_$ic]
    regexp "The length res_$ic is +(\[-0-9.+eE\]+)" ${le} full ll
    
    if { $ll < 1.0e-7 } {
      puts "Error: Curve is too small!"
    }
    
    bounds res_$ic U1 U2
    
    if {[dval U2-U1] < 1.0e-9} {
      puts "Error: Wrong curve's range!"
    }
    
    dlog reset
    dlog on
    xdistcs res_$ic s1 U1 U2 10
    set Log1 [dlog get]
    set List1 [split ${Log1} {TD= \t\n}]
    set Tolerance 1.0e-7
    set Limit_Tol 1.0e-7
    set D_good 0.
    checkList ${List1} ${Tolerance} ${D_good} ${Limit_Tol}
    
    dlog reset
    dlog on
    xdistcs res_$ic s2 U1 U2 10
    set Log1 [dlog get]
    set List1 [split ${Log1} {TD= \t\n}]
    set Tolerance 1.0e-7
    set Limit_Tol 1.0e-7
    set D_good 0.
    checkList ${List1} ${Tolerance} ${D_good} ${Limit_Tol}
    
    incr ic
  }
}
  
if {[expr {$ic - 1}] == $GoodNbCurv} {
  puts "OK: Curve Number is good!"
} else {
  puts "Error: Curve Number is bad!"
}