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
|
#! /usr/bin/env python
from __future__ import print_function
import openturns as ot
dim = 2
# We create an empty series
ts1 = ot.TimeSeries(0, dim)
ts1.setName('Ts1')
# We populate the empty ts
for p in range(3):
pt = ot.NumericalPoint(dim)
for i in range(dim):
pt[i] = 10. * (p + 1) + i
ts1.add(pt)
print('ts1=', ts1)
# We get the second element of the ts
secondElement = ts1[1]
print('second element=', secondElement)
# We set the third element to a valid new element
newPoint = ot.NumericalPoint(dim + 1)
for i in range(dim):
newPoint[i + 1] = 1000. * (i + 1)
ts1[2] = newPoint
print('ts1=', ts1)
# We get the values of the second element of the ts
values1 = ts1.getValueAtIndex(1)
print('values of the second element=', values1)
# We set new values for the second element of the ts
newValues = values1 * 2.
ts1.setValueAtIndex(1, newValues)
values2 = ts1.getValueAtIndex(1)
print('ts1=', ts1)
print('modified values of the second element=',
values2, ' should be ', newValues)
# We get the values of the second element of the ts
values3 = ts1.getValueAtNearestTime(-1.0)
print('values at t=-1.0 =', values3)
values4 = ts1.getValueAtNearestTime(1.45)
print('values at t=1.45 =', values4)
values5 = ts1.getValueAtNearestTime(1.54)
print('values at t=1.54 =', values5)
values6 = ts1.getValueAtNearestTime(14.5)
print('values at t=14.5 =', values6)
# We set new values for the third element of the ts
ts1.setValueAtNearestTime(1.54, values6 * -1.0)
print('ts1=', ts1)
ts2 = ot.TimeSeries(10, dim)
print('ts2=', ts2)
# ts2[5] = point2
# print 'ts2=', ts2
try:
# We get the tenth element of the ts
# THIS SHOULD NORMALLY FAIL
tenthElement = ts1.at(9)
except:
print('Expected failure')
tg1 = ot.RegularGrid(0.0, 0.1, 11)
ts3 = ot.TimeSeries(tg1, dim)
print('ts3=', ts3)
tg2 = ot.RegularGrid(0.0, 0.2, 6)
ts4 = ot.TimeSeries(tg2, dim)
print('ts4=', ts4)
# We append a sample to a time series
ts5 = ot.TimeSeries(3, dim)
ns1 = ot.NumericalSample(3, [99.9] * dim)
print('ts5=', ts5)
ts5.add(ns1)
print('ts5=', ts5)
# We retrieve the values of the time series as a sample
ns2 = ts5.getSample()
print('ns2=', ns2)
# We build a time series from a time grid and a sample
ns3 = ts4.getSample()
ts6 = ot.TimeSeries(tg2, ns3)
print('ts6=', ts6)
|
