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
122
123
124
125
126
127
128
129
130
|
#! /usr/bin/env python
from __future__ import print_function
from openturns import *
from math import *
TESTPREAMBLE()
try:
# DEFAULT CONSTRUCTOR AND STRING CONVERTER
print("test : default constructor and string converter")
# Default constructor
triangularMatrix0 = 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 = TriangularMatrix(2)
# Check operator() methods
triangularMatrix1[0, 0] = 1.
triangularMatrix1[1, 0] = 3.
triangularMatrix1[1, 1] = 4.
# String converter
print("triangularMatrix1 = ", repr(triangularMatrix1))
# COPY CONSTRUCTOR AND STRING CONVERTER
print("test : copy constructor and string converter")
# Copy constructor
triangularMatrix2 = 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)
# SUBSTRACTION METHOD
print("test : substraction method")
# Check substraction 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 = NumericalPoint()
pt.add(1.)
pt.add(2.)
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.
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 = TriangularMatrix()
triangularMatrix6 = TriangularMatrix()
print("triangularMatrix0 is empty = ", triangularMatrix0.isEmpty())
print("triangularMatrix1 is empty = ", triangularMatrix1.isEmpty())
print("triangularMatrix5 is empty = ", triangularMatrix5.isEmpty())
except:
import sys
print("t_TriangularMatrix_std.py", sys.exc_info()[0], sys.exc_info()[1])
|
