Note

Click here to download the full example code

Value function

A value function f_{value}: \mathcal{D} \times \mathbb{R}^d \rightarrow \mathcal{D} \times \mathbb{R}^q is a particular field function that lets invariant the mesh of a field and defined by a function g : \mathbb{R}^d \rightarrow \mathbb{R}^q such that:

\begin{aligned} f_{value}(\underline{t}, \underline{x})=(\underline{t}, g(\underline{x}))\end{aligned}

Let’s note that the input dimension of f_{value} still designs the dimension of \underline{x} : d. Its output dimension is equal to q.

The creation of the ValueFunction object requires the function g and the integer n: the dimension of the vertices of the mesh \mathcal{M}. This data is required for tests on the compatibility of dimension when a composite process is created using the spatial function.

The use case illustrates the creation of a spatial (field) function from the function g: \mathbb{R}^2 \rightarrow \mathbb{R}^2 such as :

\begin{aligned} g(\underline{x})=(x_1^2, x_1+x_2) \end{aligned}

from __future__ import print_function
import openturns as ot
import openturns.viewer as viewer
from matplotlib import pylab as plt
import math as m
ot.Log.Show(ot.Log.NONE)

Create a mesh

N = 100
mesh = ot.RegularGrid(0.0, 1.0, N)

Create the function that acts the values of the mesh

g = ot.SymbolicFunction(['x1', 'x2'],  ['x1^2', 'x1+x2'])

Create the field function

f = ot.ValueFunction(g, mesh)

Evaluate f

inF = ot.Normal(2).getSample(N)
outF = f(inF)

# print input/output at first mesh nodes
xy = inF
xy.stack(outF)
xy[:5]
X0X1y0y1
00.3255111.3342730.10595741.659784
1-1.7222110.19122892.96601-1.530982
20.6154579-0.063081960.37878840.5523759
31.20791.2796871.4590212.487586
4-2.5553950.88022156.530041-1.675173


Total running time of the script: ( 0 minutes 0.004 seconds)

Gallery generated by Sphinx-Gallery