1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
|
import openturns as ot
from matplotlib import pyplot as plt
import openturns.viewer as otv
# Create a bivariate Gaussian process
myMesh = ot.IntervalMesher([39, 39]).build(ot.Interval([0.0] * 2, [1.0] * 2))
myCov = ot.GeneralizedExponential(2 * [0.1], 1.3)
myProcess = ot.GaussianProcess(myCov, myMesh)
myField = myProcess.getRealization()
graph = myField.drawMarginal(0, False)
graph.setTitle("A field")
fig = plt.figure(figsize=(8, 4))
axis = fig.add_subplot(111)
axis.set_xlim(auto=True)
otv.View(graph, figure=fig, axes=[axis], add_legend=True)
|
