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
|
#! /usr/bin/env python
import openturns as ot
ot.TESTPREAMBLE()
basis = ot.Basis()
print("basis =", basis)
basisSize = 3
phis = []
for j in range(basisSize):
phis.append(ot.SymbolicFunction(["x"], ["x^" + str(j + 1)]))
basis = ot.Basis(phis)
print("basis =", basis)
assert basis.getSize() == 3, "wrong size"
assert basis.getInputDimension() == 1, "wrong input dim"
assert basis.getOutputDimension() == 1, "wrong output dim"
assert basis.getSubBasis([1, 2]).getSize() == 2, "wrong subsize"
assert basis.isFinite(), "!finite"
assert not basis.isOrthogonal(), "orthogonal"
try:
degree = 2
poly = []
factory = ot.MonomialFunctionFactory()
for i in range(degree + 1):
poly.append(factory.build(i))
basis = ot.Basis(poly)
raise AssertionError("should not go there")
except RuntimeError:
print("ok")
# check basis dimensions
f1 = ot.SymbolicFunction(["x", "y"], ["x+y"])
f2 = ot.SymbolicFunction(["x", "y"], ["x+y", "x-y"])
f3 = ot.SymbolicFunction(["x"], ["x+1"])
try:
b = ot.Basis([f1, f2])
raise AssertionError("should not go there")
except TypeError:
print("ok")
try:
b = ot.Basis([f1, f3])
raise AssertionError("should not go there")
except TypeError:
print("ok")
|
