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
|
#! /usr/bin/env python
import openturns as ot
import persalys
# Model with error
X0 = persalys.Input("X0", 1, ot.Normal())
X1 = persalys.Input("X1", 2)
Y0 = persalys.Output("Y0")
Y1 = persalys.Output("Y1")
model2 = persalys.PythonPhysicalModel(
"model2",
[X0, X1],
[Y0],
"from math import pi\n\ndef _exec(X0, X1):\n Y0 = 2*X1 - 1/X0\n return Y0\n",
)
model2.setParallel(True)
design5 = persalys.FixedDesignOfExperiment("design5", model2)
design5.setOriginalInputSample([[0.9, 1.1], [1.8, 2.2], [0, 2.2]])
design5.setBlockSize(3)
design5.run()
assert (
design5.getResult().getDesignOfExperiment().getOutputSample().getSize() == 2
), "whole block must no fail"
model3 = persalys.PythonPhysicalModel(
"model3",
[X0, X1],
[Y0, Y1],
"from math import pi\n\ndef _exec(X0, X1):\n Y0 = 2*X1 - 1/X0\n Y1=5\n return Y0, Y1\n",
)
model3.setParallel(True)
# same with multiple outputs
design6 = persalys.FixedDesignOfExperiment("design6", model3)
design6.setInterestVariables(["Y0"])
design6.setOriginalInputSample([[0.9, 1.1], [1.8, 2.2], [0, 2.2]])
design6.setBlockSize(3)
design6.run()
assert (
design6.getResult().getDesignOfExperiment().getOutputSample().getSize() == 2
), "whole block must no fail"
|
