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
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
|
// File: GeomPlate_PlateG0Criterion.cxx
// Created: Wed Mar 5 10:04:41 1997
// Author: Joelle CHAUVET
// <jct@sgi38>
//Modified by Jerome LEMONIER
// Tue Jun 17 10:20:35 1997
// Modification de la comparaison dans la methode
// Value (Appel a des distances remplace par le calcul
// de la distance au carre)
#include <GeomPlate_PlateG0Criterion.ixx>
#include <AdvApp2Var_Criterion.hxx>
#include <TColgp_SequenceOfXY.hxx>
#include <TColgp_SequenceOfXYZ.hxx>
#include <gp_Vec.hxx>
#include <gp_XY.hxx>
#include <gp_Pnt.hxx>
#include <TColStd_HArray1OfReal.hxx>
#include <PLib.hxx>
//============================================================================
//function : GeomPlate_PlateG0Criterion
//purpose :
//============================================================================
GeomPlate_PlateG0Criterion::
GeomPlate_PlateG0Criterion(const TColgp_SequenceOfXY& Data,
const TColgp_SequenceOfXYZ& G0Data,
const Standard_Real Maximum,
const AdvApp2Var_CriterionType Type,
const AdvApp2Var_CriterionRepartition Repart)
{
myData=Data;
myXYZ=G0Data;
myMaxValue = Maximum;
myType = Type;
myRepartition = Repart;
}
//============================================================================
//function : Value
//purpose :
//============================================================================
void GeomPlate_PlateG0Criterion::Value(AdvApp2Var_Patch& P,
const AdvApp2Var_Context& C) const
{
Standard_Real UInt[2],VInt[2];
Standard_Integer MaxNbCoeff[2], NbCoeff[2];
Standard_Real * adrCoeff = NULL ;
adrCoeff = (Standard_Real *) &P.Coefficients(1,C) ->ChangeArray1()(P.Coefficients(1,C)->Lower());
MaxNbCoeff[0] = C.ULimit();
MaxNbCoeff[1] = C.VLimit();
NbCoeff[0] = P.NbCoeffInU();
NbCoeff[1] = P.NbCoeffInV();
UInt[0] = P.U0();
UInt[1] = P.U1();
VInt[0] = P.V0();
VInt[1] = P.V1();
Standard_Real up,vp, dist = 0.;
Standard_Integer dimension = 3 * NbCoeff[1];
TColStd_Array1OfReal Patch(1, NbCoeff[0] * dimension);
TColStd_Array1OfReal Curve(1, dimension);
TColStd_Array1OfReal Point(1, 3);
Standard_Real * Coeffs = (Standard_Real *) &Patch.ChangeValue(1);
Standard_Real * Digit = (Standard_Real *) &Point.ChangeValue(1);
Standard_Integer k1, k2, pos, ll=1;
for (k1 = 1; k1 <= NbCoeff[0]; k1++) {
// JAG 99.04.29 pos = 3*(MaxNbCoeff[0])*(k1-1);
pos = 3*(MaxNbCoeff[1])*(k1-1);
for (k2 = 1; k2 <= NbCoeff[1]; k2++, pos+=3 ) {
Patch(ll) = adrCoeff[pos];
Patch(ll+1) = adrCoeff[pos+1];
Patch(ll+2) = adrCoeff[pos+2];
ll += 3;
}
}
Standard_Integer i, NbCtr = myData.Length();
for(i=1; i<=NbCtr; i++) {
gp_XY P2d = myData.Value(i);
// gp_Pnt PP = myXYZ.Value(i);
gp_Pnt P3d;
if ( UInt[0]<P2d.X() && P2d.X()<UInt[1]
&& VInt[0]<P2d.Y() && P2d.Y()<VInt[1] ) {
// u,v recadres sur (-1,1)
up = (2*P2d.X()-UInt[0]-UInt[1])
/ (UInt[1]-UInt[0]) ;
vp = (2*P2d.Y()-VInt[0]-VInt[1])
/ (VInt[1]-VInt[0]) ;
PLib::EvalPoly2Var(up,vp,
0,0,
NbCoeff[0]-1,NbCoeff[1]-1,3,
Coeffs[0],
Digit[0]);
P3d.SetCoord(1,Digit[0]);
P3d.SetCoord(2,Digit[1]);
P3d.SetCoord(3,Digit[2]);
Standard_Real x=(P3d.Coord(1)-myXYZ.Value(i).Coord(1)),
y=(P3d.Coord(2)-myXYZ.Value(i).Coord(2)),
z=(P3d.Coord(3)-myXYZ.Value(i).Coord(3)),
DistTmp= x*x+y*y+z*z;
if( DistTmp>dist ) {
dist = DistTmp;
}
}
}
P.SetCritValue(Sqrt(dist));
}
//============================================================================
//function : IsSatisfied
//purpose :
//============================================================================
Standard_Boolean GeomPlate_PlateG0Criterion::IsSatisfied(const AdvApp2Var_Patch& P) const
{
return (P.CritValue() < myMaxValue);
}
|
