File testcases/testcases_2D.py changed (mode: 100644) (index dfece1e..e70060e) |
... |
... |
h = f_models.UnCorD((stats.norm(-12,5), stats.norm(12,5))) |
16 |
16 |
#l = f_models.UnCorD([stats.lognorm(s=1), stats.lognorm(s=1)]) |
#l = f_models.UnCorD([stats.lognorm(s=1), stats.lognorm(s=1)]) |
17 |
17 |
#č pro různorodost |
#č pro různorodost |
18 |
18 |
gew = f_models.Nataf((stats.gumbel_r, stats.weibull_min(c=1.5)), [[1,0.8], [0.8,1]]) |
gew = f_models.Nataf((stats.gumbel_r, stats.weibull_min(c=1.5)), [[1,0.8], [0.8,1]]) |
|
19 |
|
gewm = f_models.Nataf((stats.gumbel_r, stats.weibull_min(c=1.5)), [[1,-0.8], [-0.8,1]]) |
19 |
20 |
#gew = f_models.Nataf((stats.norm, stats.norm), [[1,0.8], [0.8,1]]) |
#gew = f_models.Nataf((stats.norm, stats.norm), [[1,0.8], [0.8,1]]) |
20 |
21 |
#gew = f_models.Nataf((stats.gumbel_r, stats.gumbel_l), [[1,0.8], [0.8,1]]) |
#gew = f_models.Nataf((stats.gumbel_r, stats.gumbel_l), [[1,0.8], [0.8,1]]) |
21 |
22 |
#č uniforma je na [0, 1], takže v cajku |
#č uniforma je na [0, 1], takže v cajku |
|
... |
... |
The same as Linear_nD, but defined via |
83 |
84 |
beta in sense of reliability index""" |
beta in sense of reliability index""" |
84 |
85 |
# Z_sum s divokým rozdělením |
# Z_sum s divokým rozdělením |
85 |
86 |
# pro obyčejný norm dist viz. HyperPlane |
# pro obyčejný norm dist viz. HyperPlane |
86 |
|
add('sum_2D') |
|
87 |
|
def sum_2D(): |
|
88 |
|
return whitebox.WhiteBox(gew, gm.Z_sum(nvar=2, beta_exact=2)) |
|
89 |
|
|
|
|
87 |
|
add('nataf_sum_2D') |
|
88 |
|
def nataf_sum_2D(): |
|
89 |
|
wt = whitebox.WhiteBox(gew, gm.Z_sum(nvar=2, beta_exact=np.sqrt(2))) |
|
90 |
|
wt.pf_exact = 0.000495 |
|
91 |
|
wt.Nsim = int(10e6) |
|
92 |
|
wt.pf_exact_method = 'MC' |
|
93 |
|
return wt |
90 |
94 |
|
|
91 |
95 |
|
|
92 |
96 |
# soucin velicin plus nějaká konstanta |
# soucin velicin plus nějaká konstanta |
|
... |
... |
def min_2D(): |
144 |
148 |
|
|
145 |
149 |
|
|
146 |
150 |
# g = np.sum(np.exp(-(sample**2)), axis=1) + self._const |
# g = np.sum(np.exp(-(sample**2)), axis=1) + self._const |
147 |
|
add('sumexp_2D') |
|
148 |
|
def sumexp_2D(): |
|
149 |
|
return whitebox.WhiteBox(f, gm.Z_sumexp(const=-0.1)) |
|
150 |
|
#whitebox.Gaussian_Z_sumexp(nvar=2, co=-2) |
|
151 |
|
|
|
|
151 |
|
add('nataf_supexp_2D') |
|
152 |
|
def nataf_supexp_2D(): |
|
153 |
|
return whitebox.WhiteBox(gew, gm.Z_sumexp(const=-0.003)) |
|
154 |
|
add('nataf_sumexp_2D') |
|
155 |
|
def nataf_sumexp_2D(): |
|
156 |
|
wt = whitebox.WhiteBox(gewm, gm.Z_sumexp(const=-0.003)) |
|
157 |
|
wt.pf_exact = 1.2e-05 |
|
158 |
|
wt.Nsim = int(10e6) |
|
159 |
|
wt.pf_exact_method = 'MC' |
|
160 |
|
return wt |
152 |
161 |
|
|
153 |
162 |
|
|
154 |
163 |
""" |
""" |
|
... |
... |
def nataf_sumsq_2D(): |
175 |
184 |
# obecný S_ball bez výpočtu pf |
# obecný S_ball bez výpočtu pf |
176 |
185 |
add('nataf_ball_2D') |
add('nataf_ball_2D') |
177 |
186 |
def nataf_ball_2D(): |
def nataf_ball_2D(): |
178 |
|
return whitebox.WhiteBox(gew, gm.S_ball(r=5)) |
|
179 |
|
|
|
|
187 |
|
wt = whitebox.WhiteBox(gew, gm.S_ball(r=5)) |
|
188 |
|
wt.pf_exact = 0.012403 |
|
189 |
|
wt.Nsim = int(10e6) |
|
190 |
|
wt.pf_exact_method = 'MC' |
|
191 |
|
return wt |
|
192 |
|
|
|
193 |
|
add('natafm_ball_2D') |
|
194 |
|
def natafm_ball_2D(): |
|
195 |
|
wt = whitebox.WhiteBox(gewm, gm.S_ball(r=5)) |
|
196 |
|
wt.pf_exact = 0.006812 |
|
197 |
|
wt.Nsim = int(10e6) |
|
198 |
|
wt.pf_exact_method = 'MC' |
|
199 |
|
return wt |
|
200 |
|
|
180 |
201 |
|
|
181 |
202 |
|
|
182 |
203 |
|
|
|
... |
... |
add('exp_P_2D') |
202 |
223 |
def exp_P_2D(): |
def exp_P_2D(): |
203 |
224 |
return whitebox.WhiteBox(u, gm.Exp_P(k=1., pow=5)) |
return whitebox.WhiteBox(u, gm.Exp_P(k=1., pow=5)) |
204 |
225 |
|
|
205 |
|
|
|
|
226 |
|
add('nataf_exp_2D') |
|
227 |
|
def nataf_exp_2D(): |
|
228 |
|
return whitebox.WhiteBox(gew, gm.Exp_P(k=1., pow=-0.5)) |
206 |
229 |
|
|
207 |
230 |
|
|
208 |
231 |
# Sin2D |
# Sin2D |