blob: 50ea724b3510ef9a56cab4552480cd352537a907 (
plain)
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
|
// File: Plate_Plate.lxx
// Created: Thu Oct 21 9:00:00 1995
// Author: Andre LIEUTIER
// <ds@sgi63>
#include <gp_XY.hxx>
#include <gp_XYZ.hxx>
inline Standard_Real Plate_Plate::Polm(const gp_XY& point2d, const Standard_Integer iu, const Standard_Integer iv,
const Standard_Integer idu, const Standard_Integer idv) const
{
if(idu>iu) return 0;
if(idv>iv) return 0;
Standard_Real u = point2d.X();
Standard_Real v = point2d.Y();
Standard_Real value = 1;
Standard_Integer degu = iu - idu;
Standard_Integer i ;
for( i=0;i<degu;i++) value *= u;
for(i=iu;i>degu;i--) value *= i;
Standard_Integer degv = iv - idv;
for(i=0;i<degv;i++) value *= v;
for(i=iv;i>degv;i--) value *= i;
// le produit par ddu[iu]*ddv[iv] n'est pas indispensable !! (il change
// les valeurs claculles pour la partie coef polynomiaux de Sol
// de telle facon que les methodes Evaluate et EvaluateDerivative donnent
//en theorie les memes valeurs. Toutefois, il nous semble que ce produit
// ameliore le conditionnmenet de la matrice
return value*ddu[iu]*ddv[iv];
//return value;
}
inline Standard_Integer& Plate_Plate::Deru(const Standard_Integer index) const
{
return ((Standard_Integer *)deru)[index];
}
inline Standard_Integer& Plate_Plate::Derv(const Standard_Integer index) const
{
return ((Standard_Integer *)derv)[index];
}
inline gp_XYZ& Plate_Plate::Solution(const Standard_Integer index) const
{
return ((gp_XYZ *)solution)[index];
}
inline gp_XY& Plate_Plate::Points(const Standard_Integer index) const
{
return ((gp_XY *)points)[index];
}
|
