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
|
#! /usr/bin/env python
from __future__ import print_function
from openturns import *
TESTPREAMBLE()
try:
# Left hand side of the composition
left = NumericalMathFunction(
['x1', 'x2'], ['y1', 'y2', 'y3'], ['x1*sin(x2)', 'cos(x1+x2)', '(x2+1)*exp(x1-2*x2)'])
# Right hand side of the composition
right = NumericalMathFunction(['x1', 'x2', 'x3', 'x4'], ['y1', 'y2'], [
'(x1*x1+x2^3*x1)/(2*x3*x3+x4^4+1)', 'cos(x2*x2+x4)/(x1*x1+1+x3^4)'])
# Compositon of left and right
composed = NumericalMathFunction(
ComposedNumericalMathFunction(left, right))
print("right=", repr(right))
print("left=", repr(left))
print("composed=", repr(composed))
# Does it worked?
x = NumericalPoint(right.getInputDimension(), 1.0)
y = NumericalPoint(right(x))
z = NumericalPoint(left(y))
Dy = Matrix(right.gradient(x))
Dz = Matrix(left.gradient(y))
print("x=", repr(x), " y=right(x)=", repr(y), " z=left(y)=", repr(z))
print("left(right(x))=", repr(composed(x)))
print("D(right)(x)=", repr(Dy), " D(left)(y)=", repr(Dz))
print(" prod=", repr(Dy * Dz))
print("D(left(right(x)))=", repr(composed.gradient(x)))
result = composed.hessian(x)
print("DD(left(right(x)))=")
for k in range(result.getNbSheets()):
for j in range(result.getNbColumns()):
for i in range(result.getNbRows()):
print("%.6f" % result[i, j, k])
print("")
print("")
except:
import sys
print("t_NumericalMathFunction_composition.py",
sys.exc_info()[0], sys.exc_info()[1])
|
