blob: e9e7a04f0aceea5555dc58ef4a157338fa42584a (
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
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
|
-- File: FunctionRoot.cdl
-- Created: Tue 14 11:02:43 1991
-- Author: Laurent PAINNOT
-- <lpa@topsn3>
---Copyright: Matra Datavision 1991, 1992
--
class FunctionRoot from math
---Purpose:
-- This class implements the computation of a root of a function of
-- a single variable which is near an initial guess using a minimization
-- algorithm.Knowledge of the derivative is required. The
-- algorithm used is the same as in
uses Vector from math, Matrix from math,
FunctionWithDerivative from math,
OStream from Standard
raises NotDone from StdFail
is
Create(F: in out FunctionWithDerivative;
Guess, Tolerance: Real;
NbIterations: Integer = 100)
---Purpose:
-- The Newton-Raphson method is done to find the root of the function F
-- from the initial guess Guess.The tolerance required on
-- the root is given by Tolerance. Iterations are stopped if
-- the expected solution does not stay in the range A..B.
-- The solution is found when abs(Xi - Xi-1) <= Tolerance;
-- The maximum number of iterations allowed is given by NbIterations.
returns FunctionRoot;
Create(F: in out FunctionWithDerivative;
Guess, Tolerance,A,B: Real;
NbIterations: Integer = 100)
---Purpose:
-- The Newton-Raphson method is done to find the root of the function F
-- from the initial guess Guess.
-- The tolerance required on the root is given by Tolerance.
-- Iterations are stopped if the expected solution does not stay in the
-- range A..B
-- The solution is found when abs(Xi - Xi-1) <= Tolerance;
-- The maximum number of iterations allowed is given by NbIterations.
returns FunctionRoot;
IsDone(me)
---Purpose: Returns true if the computations are successful, otherwise returns false.
---C++: inline
returns Boolean
is static;
Root(me)
---Purpose: returns the value of the root.
-- Exception NotDone is raised if the root was not found.
---C++: inline
returns Real
raises NotDone
is static;
Derivative(me)
---Purpose: returns the value of the derivative at the root.
-- Exception NotDone is raised if the root was not found.
---C++: inline
returns Real
raises NotDone
is static;
Value(me)
---Purpose: returns the value of the function at the root.
-- Exception NotDone is raised if the root was not found.
---C++: inline
returns Real
raises NotDone
is static;
NbIterations(me)
---Purpose: returns the number of iterations really done on the
-- computation of the Root.
-- Exception NotDone is raised if the root was not found.
---C++: inline
returns Integer
raises NotDone
is static;
Dump(me; o: in out OStream)
---Purpose: Prints on the stream o information on the current state
-- of the object.
-- Is used to redefine the operator <<.
is static;
fields
Done: Boolean;
TheRoot: Real ;
TheError: Real ;
TheDerivative: Real ;
NbIter: Integer;
end FunctionRoot;
|