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
|
"""
Compute the L2 error between two functions
==========================================
"""
# %%
# %%
# In this example we compute the L2 error between two functions.
# %%
import openturns as ot
from openturns.usecases import ishigami_function
# %%
# This is easy using the arithmetic operator `-` on two
# instances of the :class:`~openturns.Function` class.
# %%
im = ishigami_function.IshigamiModel()
partOfIshigamiFunction = ot.SymbolicFunction(
["x1", "x2", "x3"], ["7 * (sin(x2))^2 + 0.1 * x3^4 * sin(x1)"]
)
sampleSize = 1000
experiment = ot.MonteCarloExperiment(im.distribution, sampleSize)
integration = ot.ExperimentIntegration(experiment)
functionError = integration.computeL2Norm(im.model - partOfIshigamiFunction)
print(functionError)
|
