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
|
#! /usr/bin/env python
import openturns as ot
ot.TESTPREAMBLE()
# DEFAULT CONSTRUCTOR AND STRING CONVERTER
print("test : default constructor and string converter")
# Default constructor
triangularMatrix0 = ot.TriangularMatrix()
# String converter
print("triangularMatrix0 = ", repr(triangularMatrix0))
# CONSTRUCTOR WITH SIZE, OPERATOR() AND STRING CONVERTER
print("test : constructor with size, operator() and string converter")
# Constructor with size
triangularMatrix1 = ot.TriangularMatrix(2)
# Check operator() methods
triangularMatrix1[0, 0] = 1.0
triangularMatrix1[1, 0] = 3.0
triangularMatrix1[1, 1] = 4.0
# String converter
print("triangularMatrix1 = ", repr(triangularMatrix1))
# COPY CONSTRUCTOR AND STRING CONVERTER
print("test : copy constructor and string converter")
# Copy constructor
triangularMatrix2 = ot.TriangularMatrix(triangularMatrix1)
# String converter
print("triangularMatrix2 = ", repr(triangularMatrix2))
# GET DIMENSIONS METHODS
print("test : get dimensions methods")
# Get dimension methods
print("triangularMatrix1's nbRows = ", triangularMatrix1.getNbRows())
print("triangularMatrix1's nbColumns = ", triangularMatrix1.getNbColumns())
print("isLower? ", triangularMatrix1.isLowerTriangular())
# TRANSPOSE AND CONJUGATE METHOD
print("test : transpose and conjugate methods")
print("triangularMatrix1 transpose = ", repr(triangularMatrix1.transpose()))
print("is transpose Lower? ", triangularMatrix1.transpose().isLowerTriangular())
triangularMatrix4 = triangularMatrix1.transpose()
print("triangularMatrix1 conjugate and transpose = ", repr(triangularMatrix4))
print("is conjugate Lower? ", triangularMatrix4.isLowerTriangular())
# ADDITION METHOD
print("test : addition method")
# Check addition method : we check the operator and the symmetry of the
# operator, thus testing the comparison operator
sum1 = triangularMatrix1 + triangularMatrix4
sum2 = triangularMatrix4 + triangularMatrix1
print("sum1 = ", repr(sum1))
print("sum2 = ", repr(sum2))
print("sum1 equals sum2 = ", sum1 == sum2)
# SUBTRACTION METHOD
print("test : subtraction method")
# Check subtraction method
diff = triangularMatrix1 - triangularMatrix4
print("diff = ", repr(diff))
# MATRIX MULTIPLICATION METHOD
print("test : matrix multiplication method")
# Check multiplication method
prod = triangularMatrix1 * triangularMatrix4
print("prod = ", repr(prod))
# MULTIPLICATION WITH A NUMERICAL POINT METHOD
print("test : multiplication with a numerical point method")
# Create the numerical point
pt = ot.Point()
pt.add(1.0)
pt.add(2.0)
print("pt = ", repr(pt))
# Check the product method
ptResult = triangularMatrix1 * pt
print("ptResult = ", repr(ptResult))
# MULTIPLICATION AND DIVISION BY A NUMERICAL SCALAR METHODS
print("test : multiplication and division by a numerical scalar methods")
# Check the multiplication method
s = 3.0
scalprod = triangularMatrix1 * s
print("scalprod = ", repr(scalprod))
# Check the division method
scaldiv1 = triangularMatrix1 / s
scaldiv2 = triangularMatrix1 / s
print("scaldiv1 = ", repr(scaldiv1))
print("scaldiv2 = ", repr(scaldiv2))
print("scaldiv1 equals scaldiv2 = ", (scaldiv1 == scaldiv2))
# ISEMPTY METHOD
print("test : isEmpty method")
# Check method isEmpty
triangularMatrix5 = ot.TriangularMatrix()
triangularMatrix6 = ot.TriangularMatrix()
print("triangularMatrix0 is empty = ", triangularMatrix0.isEmpty())
print("triangularMatrix1 is empty = ", triangularMatrix1.isEmpty())
print("triangularMatrix5 is empty = ", triangularMatrix5.isEmpty())
|
