Note
Go to the end to download the full example code.
Simulate an EventΒΆ
In this example we are going to define an Event from a scalar variable in the form:
with a scalar threshold
import openturns as ot
ot.Log.Show(ot.Log.NONE)
Create model f(x) = x1 + 2*x2
model = ot.SymbolicFunction(["x1", "x2"], ["x1+2*x2"])
# Create the input distribution and random vector X
inputDist = ot.Normal(2)
inputDist.setDescription(["X1", "X2"])
inputVector = ot.RandomVector(inputDist)
# Create the output random vector Y=f(X)
outputVector = ot.CompositeRandomVector(model, inputVector)
Create the event Y > 3
threshold = 3.0
event = ot.ThresholdEvent(outputVector, ot.Greater(), threshold)
Realization as a Bernoulli
print("realization=", event.getRealization())
realization= [0]
Sample of 10 realizations as a Bernoulli
print("sample=", event.getSample(10))
sample= [ y0 ]
0 : [ 0 ]
1 : [ 0 ]
2 : [ 0 ]
3 : [ 0 ]
4 : [ 0 ]
5 : [ 0 ]
6 : [ 0 ]
7 : [ 0 ]
8 : [ 0 ]
9 : [ 0 ]
Build a standard event based on an event
standardEvent = ot.StandardEvent(event)