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
95
96
|
#!/usr/bin/env python
import os
import sys
import persalys
# import openturns as ot
import openturns.testing as ott
ssh_hostname = ""
with open("input_logistic.txt.in", "w") as f:
f.write("y0=@y0@\n")
f.write("a=@a@\n")
f.write("b=@b@\n")
with open("program_logistic.py", "w") as f:
f.write("import math as m\n")
f.write("def logisticSolution(t, y0, a, b):\n")
f.write(" t0 = 1790.0\n")
f.write(" y = a * y0 / (b * y0 + (a - b * y0) * m.exp(-a * (t - t0)))\n")
f.write(" return y\n\n")
f.write('exec(open("input_logistic.txt").read())\n')
f.write("Y = [logisticSolution(t, y0, a, b) for t in range(1790, 2011)]\n")
f.write('with open("Y.txt", "w") as f:\n')
f.write(" for y in Y:\n")
f.write(' f.write("%.17g\\n" % y)\n')
with open("post_logistic.py", "w") as f:
f.write("import math as m\n")
f.write('with open("Y.txt") as f:\n')
f.write(" Y = [float(line) for line in f.readlines()]\n")
f.write("Y_min = min(Y)\n")
f.write("Y_max = max(Y)\n")
f.write("Y_mean = sum(Y) / len(Y)\n")
f.write("Y_last = Y[-1]\n")
f.write("try:\n idx = next(i for i,y in enumerate(Y) if y > 100e6)\n")
f.write(
" Y_mean_t = sum(Y[:idx]) / idx\nexcept StopIteration:\n Y_mean_t = 0.0\n"
)
f.write('with open("output.txt", "w") as f:\n')
f.write(' f.write("Y_min=%.17g\\n" % Y_min)\n')
f.write(' f.write("Y_max=%.17g\\n" % Y_max)\n')
f.write(' f.write("Y_mean=%.17g\\n" % Y_mean)\n')
f.write(' f.write("Y_last=%.17g\\n" % Y_last)\n')
f.write(' f.write("Y_mean_t=%.17g\\n" % Y_mean_t)\n')
input_file1 = persalys.CouplingInputFile("input_logistic.txt.in")
input_file1.setConfiguredPath("input_logistic.txt")
input_file1.setVariables(
["y0", "a", "b"], ["@y0@", "@a@", "@b@"], ["{:10.6f}", "{:10.6f}", "{:10.6f}"]
)
resource_file1 = persalys.CouplingResourceFile("program_logistic.py")
step1 = persalys.CouplingStep(
sys.executable + " program_logistic.py", [input_file1], [resource_file1], []
)
resource_file2 = persalys.CouplingResourceFile("post_logistic.py")
output_file2 = persalys.CouplingOutputFile("output.txt")
output_file2.setVariables(
["Y_min", "Y_max", "Y_mean", "Y_last", "Y_mean_t"],
["Y_min=", "Y_max=", "Y_mean=", "Y_last=", "Y_mean_t="],
[0] * 5,
[0] * 5,
[0] * 5,
)
step2 = persalys.CouplingStep(
sys.executable + " post_logistic.py", [], [resource_file2], [output_file2]
)
model = persalys.CouplingPhysicalModel("logistic", [step1, step2])
model.setSSHHostname(ssh_hostname)
print(model.getInputNames())
print(model.getOutputNames())
assert model.getInputNames() == ["y0", "a", "b"]
assert model.getOutputNames() == ["Y_min", "Y_max", "Y_mean", "Y_last", "Y_mean_t"]
# single evaluation
x = [3.9e6, 0.03134, 1.5887e-10]
f = model.getFunction()
y = f(x)
print(y)
ott.assert_almost_equal(y, [3.9e06, 3.8498e09, 5.64037e08, 3.8498e09, 2.94863e07])
# cleanup
os.remove("input_logistic.txt.in")
os.remove("program_logistic.py")
os.remove("post_logistic.py")
# script
myStudy = persalys.Study("myStudy")
myStudy.add(model)
script = myStudy.getPythonScript()
# print('script=', script)
exec(script)
|
