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
|
#! /usr/bin/env python
import openturns as ot
ref_values = [[1.0, 0.0], [0.0, 0.5]]
mats = [
ot.Matrix(ref_values),
ot.SquareMatrix(ref_values),
ot.TriangularMatrix(ref_values),
ot.SymmetricMatrix(ref_values),
ot.CovarianceMatrix(ref_values),
ot.CorrelationMatrix(ref_values),
]
mats.extend(
[
ot.ComplexMatrix(ref_values),
ot.HermitianMatrix(ref_values),
ot.TriangularComplexMatrix(ref_values),
ot.SquareComplexMatrix(ref_values),
]
)
for a in mats:
# try conversion
ref = ot.Matrix([[1.0, 0.0], [0.0, 0.5]])
iname = a.__class__.__name__
print("a=", a)
# try scalar mul
try:
s = 5.0
ats = a * s
print("a*s=", ats)
sta = s * a
print("s*a=", sta)
except Exception:
print("no scalar mul for", iname)
# try scalar div
try:
s = 5.0
ads = a / s
print("a/s=", ads)
except Exception:
print("no scalar div for", iname)
# try vec mul
try:
x = ot.Point(range(6, 6 + a.getNbColumns()))
ax = a * x
print("a*x=", ax)
except Exception:
print("no vec mul for", iname)
try:
a3 = a**3
print("a**3=", a3)
except Exception:
print("no pow for", iname)
for b in mats:
jname = b.__class__.__name__
try:
ab = a * b
print("a*b=", ab)
except Exception:
print("no mul for", iname, jname)
try:
apb = a + b
print("a+b=", apb)
except Exception:
print("no add for", iname, jname)
try:
amb = a - b
print("a-b=", amb)
except Exception:
print("no sub for", jname)
sm = ot.SymmetricMatrix(2)
sm[0, 0] = 2.0
sm[1, 1] = 3.0
sm[0, 1] = 0.3
w = ot.Point([1.0, 2.0])
sq = ot.SquareMatrix(sm)
print("sm * w=", sm * w)
print("sq * w=", sq * w)
sm = ot.HermitianMatrix(2)
sm[0, 0] = 2.0
sm[1, 1] = 3.0
sm[1, 0] = 0.3j
w = ot.Point([1.0, 2.0])
sq = ot.SquareComplexMatrix(sm)
print("sm * w=", sm * w)
print("sq * w=", sq * w)
|
