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
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
|
#! /usr/bin/env python
import openturns as ot
ot.TESTPREAMBLE()
# Default constructor */
point1 = ot.Point()
ref_point1 = point1
# Check method add()
point1.add(0.0)
point1.add(1.0)
size = ref_point1.getDimension()
print("size of point1 = ", size)
val1 = ref_point1[0]
val2 = ref_point1[1]
print("point1 = ", repr(ref_point1))
print("point1[0] = %.6f" % val1)
print("point1[1] = %.6f" % val2)
# Constructor with size
point2 = ot.Point(2)
ref_point2 = point2
# Check operator[] methods
point2[0] = 10.0
point2[1] = 11.0
val1 = ref_point2[0]
val2 = ref_point2[1]
print("point2[0] = %.6f" % val1)
print("point2[1] = %.6f" % val2)
# Copy constructor
point3 = ot.Point(ref_point1)
ref_point3 = point3
val1 = ref_point3[0]
val2 = ref_point3[1]
print("point3[0] = %.6f" % val1)
print("point3[1] = %.6f" % val2)
# Constructor from indices
indices = ot.Indices(5)
indices.fill()
point = ot.Point(indices)
print("point=", point)
# Constructor from python sequence
point4 = ot.Point((1000.0, 1001.0))
i = 0
for coord in point4:
print("point4[%d] = %.6f" % (i, coord))
i += 1
# Comparison operator
if ref_point2 != ref_point2:
print(
"OT::Base::Type::Point.operator == does NOT return the correct value. Says that point2 and point2 are DIFFERENT though they are equal."
)
raise RuntimeError()
if ref_point2 == ref_point3:
print(
"OT::Base::Type::Point.operator == does NOT return the correct value. Says that point2 and point3 are EQUAL though they are different."
)
raise RuntimeError()
# Addition/Subtraction/Product operators
point5 = ref_point1 + ref_point2
ref_point5 = point5
print("point5 = ", repr(ref_point5))
point6 = ref_point1 - ref_point2
ref_point6 = point6
print("point6 = ", repr(ref_point6))
point7 = ref_point5
point7 += ref_point1 + ref_point2
ref_point7 = point7
print("point7 = ", repr(ref_point7))
point8 = ref_point6
point8 -= ref_point1 - ref_point2
ref_point8 = point8
print("point8 = ", repr(ref_point8))
point9 = ot.Point(2, 5.0) / 2.0
print("point9 = ", point9)
point10 = ot.Point(2, 5.0)
point10 /= 2.0
print("point10 = ", point10)
# We create an orthonormal base (O,i,j)
i = ot.Point(2)
i[0] = 1
i[1] = 0
print("i = ", repr(i))
j = ot.Point(2)
j[0] = 0
j[1] = 1
print("j = ", repr(j))
# BUG PYTHON for I = Point(2. * i)
II = ot.Point(i * 2.0)
print("I = ", repr(II))
J = ot.Point(j * 3.0)
print("J = ", repr(J))
dotProduct = II.dot(J)
print("dotProduct = %.6f" % dotProduct)
# slicing
point8 = ot.Point((1.0, 2.0, 3.0, 4.0, 5.0))
print("point8[2] = ", point8[2])
print("point8[1:3] = ", point8[1:3])
print("point8[:] = ", point8[:])
point8[1:3] = (88, 99)
print("point8 = ", point8)
point8[0:5] = ot.Point(5, 9)
print("point8 = ", point8)
# equality
pt1 = ot.Point([1.1, 2.2, 3.3])
pt2 = ot.Point((1.1, 2.2, 3.3))
print("Equality ? ", (pt1 == pt2))
# non-regression for #754
pt = ot.Point([-3, -5, 6])
print("pt = ", pt)
# last
print("pt[-1]=", pt[-1])
pt[-1] = 8.0
print("pt[-1]=", pt[-1])
print(0 + pt1)
try:
1 + pt1
except Exception:
print("1+Point -> ok")
# unary minus
x = ot.Point([1, 2, 3])
print(-x)
|
