blob: 373cc4592a511923063c4891ce7a24c504611183 (
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
|
-- File: FEmTool_ProfileMatrix.cdl
-- Created: Wed Oct 29 16:54:05 1997
-- Author: Roman BORISOV
-- <rbv@velox.nnov.matra-dtv.fr>
---Copyright: Matra Datavision 1997
class ProfileMatrix from FEmTool inherits SparseMatrix from FEmTool
---Purpose: Symmetric Sparse ProfileMatrix useful for 1D Finite
-- Element methods
uses
HArray1OfInteger from TColStd,
Array1OfInteger from TColStd,
Array2OfInteger from TColStd,
HArray1OfReal from TColStd,
Vector from math
raises
NotDone from StdFail,
NotImplemented from Standard,
OutOfRange from Standard
is
Create(FirstIndexes : Array1OfInteger)
returns mutable ProfileMatrix from FEmTool;
Init(me: mutable; Value : Real);
ChangeValue(me: mutable; I, J : Integer)
---C++: return &
returns Real
raises OutOfRange;
Decompose(me : mutable)
---Purpose: To make a Factorization of <me>
returns Boolean;
Solve(me; B : Vector; X : in out Vector)
---Purpose: Direct Solve of AX = B
raises NotDone from StdFail; -- if <me> is not decomposed
Prepare(me : mutable)
---Purpose: Make Preparation to iterative solve
returns Boolean
raises NotImplemented from Standard;
Solve(me; B : Vector;
Init : Vector;
X : out Vector;
Residual : out Vector;
Tolerance : Real = 1.0e-8;
NbIterations: Integer = 50)
---Purpose: Iterative solve of AX = B
raises NotDone from StdFail; -- if <me> is not prepared;
Multiplied(me; X: Vector; MX : in out Vector);
---Purpose: returns the product of a SparseMatrix by a vector.
-- An exception is raised if the dimensions are different
RowNumber(me)
---Purpose: returns the row range of a matrix.
returns Integer;
ColNumber(me)
---Purpose: returns the column range of the matrix.
returns Integer;
IsInProfile(me; i, j : Integer)
returns Boolean;
-- for debug
OutM(me);
OutS(me);
fields
profile : Array2OfInteger; -- Like MPOSIT in Fortran
ProfileMatrix : HArray1OfReal; -- Like AMATRI in Fortran
SMatrix : HArray1OfReal; -- Like SMATRI in Fortran
NextCoeff : HArray1OfInteger; -- Like POSUIV in Fortran
IsDecomp : Boolean;
end ProfileMatrix;
|