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
|
"""
Aggregate processes
===================
"""
# %%
# In this example we are going to concatenate several processes that share the same mesh.
# %%
import openturns as ot
import openturns.viewer as otv
# %%
# Create processes to aggregate
myMesher = ot.IntervalMesher([100, 10])
lowerbound = [0.0, 0.0]
upperBound = [2.0, 4.0]
myInterval = ot.Interval(lowerbound, upperBound)
myMesh = myMesher.build(myInterval)
myProcess1 = ot.WhiteNoise(ot.Normal(), myMesh)
myProcess2 = ot.WhiteNoise(ot.Triangular(), myMesh)
# %%
# Draw values of a realization of the 2nd process
marginal = ot.HistogramFactory().build(myProcess1.getRealization().getValues())
graph = marginal.drawPDF()
view = otv.View(graph)
# %%
# Create an aggregated process
myAggregatedProcess = ot.AggregatedProcess([myProcess1, myProcess2])
# %%
# Draw values of the realization on the 2nd marginal
marginal = ot.HistogramFactory().build(
myAggregatedProcess.getRealization().getValues().getMarginal(0)
)
graph = marginal.drawPDF()
otv.View(graph)
# %%
# Display all figures
otv.View.ShowAll()
|
