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
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
|
class=ImportanceSamplingAnalysis name=myIS limitState=class=LimitStateImplementation name=aLimitState physicalModel=aModelPhys outputName=Y0 operator=class=Greater name=Unnamed threshold=20 maximumCalls=1000 maximumCoefficientOfVariation=0.1 maximumElapsedTime=60 seed=2 blockSize=10 standardSpaceDesignPoint=class=Point name=Unnamed dimension=2 values=[0,0]
result= class=SimulationReliabilityResult name=Unnamed simulationResult=probabilityEstimate=6.900000e-02 varianceEstimate=7.073900e-05 standard deviation=8.41e-03 coefficient of variation=1.22e-01 confidenceLength(0.95)=3.30e-02 outerSampling=100 blockSize=10
class=ImportanceSamplingAnalysis name=myIS2 limitState=class=LimitStateImplementation name=aLimitState physicalModel=aModelPhys outputName=Y0 operator=class=Greater name=Unnamed threshold=20 maximumCalls=inf maximumCoefficientOfVariation=0.02 maximumElapsedTime=100000 seed=0 blockSize=100 standardSpaceDesignPoint=class=Point name=Unnamed dimension=2 values=[0,0]
result= class=SimulationReliabilityResult name=Unnamed simulationResult=probabilityEstimate=7.684211e-02 varianceEstimate=2.356834e-06 standard deviation=1.54e-03 coefficient of variation=2.00e-02 confidenceLength(0.95)=6.02e-03 outerSampling=304 blockSize=100
class=ImportanceSamplingAnalysis name=myIS3 limitState=class=LimitStateImplementation name=aLimitState physicalModel=aModelPhys outputName=Y0 operator=class=Greater name=Unnamed threshold=20 maximumCalls=1000 maximumCoefficientOfVariation=0.1 maximumElapsedTime=60 seed=0 blockSize=10 standardSpaceDesignPoint=class=Point name=Unnamed dimension=2 values=[0,0]
result= class=SimulationReliabilityResult name=Unnamed simulationResult=probabilityEstimate=1.272727e-01 varianceEstimate=1.586825e-04 standard deviation=1.26e-02 coefficient of variation=9.90e-02 confidenceLength(0.95)=4.94e-02 outerSampling=77 blockSize=10
class=ImportanceSamplingAnalysis name=myIS4 limitState=class=LimitStateImplementation name=aLimitState physicalModel=aModelPhys outputName=Y0 operator=class=Greater name=Unnamed threshold=20 maximumCalls=1000 maximumCoefficientOfVariation=-1 maximumElapsedTime=60 seed=0 blockSize=10 standardSpaceDesignPoint=class=Point name=Unnamed dimension=2 values=[0.0686039,1.14043]
result= class=SimulationReliabilityResult name=Unnamed simulationResult=probabilityEstimate=1.124246e-01 varianceEstimate=2.365465e-05 standard deviation=4.86e-03 coefficient of variation=4.33e-02 confidenceLength(0.95)=1.91e-02 outerSampling=100 blockSize=10
#!/usr/bin/env python
import openturns as ot
import persalys
myStudy = persalys.Study('myStudy')
persalys.Study.Add(myStudy)
dist_X0 = ot.Normal(1, 1)
X0 = persalys.Input('X0', 0, dist_X0, '')
dist_X1 = ot.Normal(1, 1)
X1 = persalys.Input('X1', 0, dist_X1, '')
X2 = persalys.Input('X2', 2, '')
Y0 = persalys.Output('Y0', '')
inputs = [X0, X1, X2]
outputs = [Y0]
formulas = ['sin(X0) + 8*X1 + X2']
aModelPhys = persalys.SymbolicPhysicalModel('aModelPhys', inputs, outputs, formulas)
myStudy.add(aModelPhys)
aLimitState = persalys.LimitState('aLimitState', aModelPhys, 'Y0', ot.Greater(), 20)
myStudy.add(aLimitState)
myIS = persalys.ImportanceSamplingAnalysis('myIS', aLimitState)
myIS.setMaximumCalls(1000)
myIS.setMaximumCoefficientOfVariation(0.1)
myIS.setMaximumElapsedTime(60)
myIS.setBlockSize(10)
myIS.setSeed(2)
myIS.setStandardSpaceDesignPoint([0,0])
myStudy.add(myIS)
myIS2 = persalys.ImportanceSamplingAnalysis('myIS2', aLimitState)
myIS2.setMaximumCoefficientOfVariation(0.02)
myIS2.setMaximumElapsedTime(100000)
myIS2.setBlockSize(100)
myIS2.setSeed(0)
myIS2.setStandardSpaceDesignPoint([0,0])
myStudy.add(myIS2)
myIS3 = persalys.ImportanceSamplingAnalysis('myIS3', aLimitState)
myIS3.setMaximumCalls(1000)
myIS3.setMaximumCoefficientOfVariation(0.1)
myIS3.setMaximumElapsedTime(60)
myIS3.setBlockSize(10)
myIS3.setSeed(0)
myIS3.setStandardSpaceDesignPoint([0,0])
myStudy.add(myIS3)
myIS4 = persalys.ImportanceSamplingAnalysis('myIS4', aLimitState)
myIS4.setMaximumCalls(1000)
myIS4.setMaximumCoefficientOfVariation(-1)
myIS4.setMaximumElapsedTime(60)
myIS4.setBlockSize(10)
myIS4.setSeed(0)
myIS4.setStandardSpaceDesignPoint([0.0686039,1.14043])
myStudy.add(myIS4)
|
