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
|
#! /usr/bin/env python
from __future__ import print_function
from openturns import *
from math import *
TESTPREAMBLE()
try:
# Create an intance
inputVars = Description(["t"])
formula = Description(["sin(t)", "cos(t)"])
myFunc = NumericalMathFunction(inputVars, formula)
myInverseTrendFunc = InverseTrendTransform(myFunc)
print("myInverseTrendFunc=", myInverseTrendFunc)
# Get the input description and dimension
print("myInverseTrendFunc input description=",
myInverseTrendFunc.getInputDescription())
print("myInverseTrendFunc input dimension=",
myInverseTrendFunc.getInputDimension())
# Get the output description and dimension
print("myInverseTrendFunc output description=",
myInverseTrendFunc.getOutputDescription())
print("myInverseTrendFunc output dimension=",
myInverseTrendFunc.getOutputDimension())
# Create a TimeSeries
tg = RegularGrid(0.0, 0.1, 11)
data = NumericalSample(tg.getN(), formula.getSize())
for i in range(data.getSize()):
t = tg.getStart() + i * tg.getStep()
data[i, 0] = i + 1 + sin(t)
data[i, 1] = i + cos(t)
ts = TimeSeries(tg, data)
print("input time series = ")
print(ts)
print("output time series = ")
print(myInverseTrendFunc(ts))
# Get the number of calls
print("called ", myInverseTrendFunc.getCallsNumber(), " times")
except:
import sys
print("t_InverseTrendTransform_std.py",
sys.exc_info()[0], sys.exc_info()[1])
|
