File: SmolyakExperiment.py

package info (click to toggle)
openturns 1.26-4
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 67,708 kB
  • sloc: cpp: 261,605; python: 67,030; ansic: 4,378; javascript: 406; sh: 185; xml: 164; makefile: 101
file content (26 lines) | stat: -rw-r--r-- 680 bytes parent folder | download
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
import openturns as ot
from matplotlib import pyplot as plt
import openturns.viewer as otv


experiment1 = ot.GaussProductExperiment(ot.Uniform(0.0, 1.0))
experiment2 = ot.GaussProductExperiment(ot.Uniform(0.0, 1.0))
collection = [experiment1, experiment2]
level = 3
smolyak = ot.SmolyakExperiment(collection, level)

sample = smolyak.generate()

# Create an empty graph
graph = ot.Graph("Smolyak experiment", "x1", "x2", True, "")

# Create the cloud
cloud = ot.Cloud(sample, "blue", "fsquare", "")

# Then, draw it
graph.add(cloud)

fig = plt.figure(figsize=(4, 4))
axis = fig.add_subplot(111)
axis.set_xlim(auto=True)
otv.View(graph, figure=fig, axes=[axis], add_legend=False)