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
|
// File: Law_Linear.cxx
// Created: Thu Jan 12 10:37:57 1995
// Author: Laurent BOURESCHE
// <lbo@phylox>
// pmn -> 17/01/1996 added : Continuity, (Nb)Interals, D2, Trim
#include <Law_Linear.ixx>
Law_Linear::Law_Linear () {}
void Law_Linear::Set (const Standard_Real Pdeb,
const Standard_Real Valdeb,
const Standard_Real Pfin,
const Standard_Real Valfin)
{
pdeb = Pdeb;
pfin = Pfin;
valdeb = Valdeb;
valfin = Valfin;
}
//=======================================================================
//function : Continuity
//purpose :
//=======================================================================
GeomAbs_Shape Law_Linear::Continuity() const
{
return GeomAbs_CN;
}
//=======================================================================
//function : NbIntervals
//purpose :
//=======================================================================
//Standard_Integer Law_Linear::NbIntervals(const GeomAbs_Shape S) const
Standard_Integer Law_Linear::NbIntervals(const GeomAbs_Shape ) const
{
return 1;
}
//=======================================================================
//function : Intervals
//purpose :
//=======================================================================
void Law_Linear::Intervals(TColStd_Array1OfReal& T,
// const GeomAbs_Shape S) const
const GeomAbs_Shape ) const
{
T.SetValue(T.Lower(), pdeb);
T.SetValue(T.Upper(), pfin);
}
Standard_Real Law_Linear::Value(const Standard_Real X)
{
return ((X-pdeb)*valfin + (pfin-X)*valdeb)/(pfin-pdeb);
}
void Law_Linear::D1(const Standard_Real X,
Standard_Real& F,
Standard_Real& D)
{
F = ((X-pdeb)*valfin + (pfin-X)*valdeb)/(pfin-pdeb);
D = (valfin-valdeb)/(pfin-pdeb);
}
void Law_Linear::D2(const Standard_Real X,
Standard_Real& F,
Standard_Real& D,
Standard_Real& D2)
{
F = ((X-pdeb)*valfin + (pfin-X)*valdeb)/(pfin-pdeb);
D = (valfin-valdeb)/(pfin-pdeb);
D2 = 0;
}
//=======================================================================
//function : Trim
//purpose :
//=======================================================================
Handle(Law_Function) Law_Linear::Trim(const Standard_Real PFirst,
const Standard_Real PLast,
// const Standard_Real Tol) const
const Standard_Real ) const
{
Handle(Law_Linear) l = new (Law_Linear)();
Standard_Real Vdeb, Vfin;
Vdeb = (( PFirst-pdeb)*valfin + (pfin-PFirst)*valdeb)/(pfin-pdeb);
Vfin = (( PLast-pdeb)*valfin + (pfin-PLast)*valdeb)/(pfin-pdeb);
l->Set(PFirst, Vdeb, PLast, Vfin);
return l;
}
void Law_Linear::Bounds(Standard_Real& PFirst,
Standard_Real& PLast)
{
PFirst = pdeb;
PLast = pfin;
}
|
