Mercurial > repos > pointsource_pde / file comparison
comparison: experiments/laser_and_mirrors_aux.py
experiments/laser_and_mirrors_aux.py
- changeset 18
- ddd24b5c2715
- parent 12
- a0d61ac3719e
equal
deleted
inserted
replaced
| 17:354747ed506c | 18:ddd24b5c2715 |
|---|---|
| 243 c1_noisy=c1_noisy, | 243 c1_noisy=c1_noisy, |
| 244 c2_noisy=c2_noisy, | 244 c2_noisy=c2_noisy, |
| 245 true_u_n_list_array=[u_n_list[i].x.array for i in plotter.plot_frames], | 245 true_u_n_list_array=[u_n_list[i].x.array for i in plotter.plot_frames], |
| 246 ) | 246 ) |
| 247 | 247 |
| 248 delta = 0.001 * k | 248 # delta = 0.001 * k |
| 249 dk_est = ( | 249 # dk_est = ( |
| 250 dat.apply((μ_true, (k + delta, (c1, c2)))) - dat.apply((μ_true, auxtrue)) | 250 # dat.apply((μ_true, (k + delta, (c1, c2)))) - dat.apply((μ_true, auxtrue)) |
| 251 ) / delta | 251 # ) / delta |
| 252 print("k-diff-test:", dk, ", ", dk_est) | 252 # print("k-diff-test:", dk, ", ", dk_est) |
| 253 delta = 0.001 * c1 | 253 # delta = 0.001 * c1 |
| 254 dc1_est = ( | 254 # dc1_est = ( |
| 255 dat.apply((μ_true, (k, (c1 + delta, c2)))) - dat.apply((μ_true, auxtrue)) | 255 # dat.apply((μ_true, (k, (c1 + delta, c2)))) - dat.apply((μ_true, auxtrue)) |
| 256 ) / delta | 256 # ) / delta |
| 257 print("c1-diff-test:", dc1, ", ", dc1_est) | 257 # print("c1-diff-test:", dc1, ", ", dc1_est) |
| 258 delta = 0.001 | 258 # delta = 0.001 |
| 259 μ_perturb = DiscreteMeasure_2_f64( | 259 # μ_perturb = DiscreteMeasure_2_f64( |
| 260 [([0.1, 0.3], 0.08 + delta), ([0.4, 0.25], 0.05), ([0.25, 0.13], 0.06)] | 260 # [([0.1, 0.3], 0.08 + delta), ([0.4, 0.25], 0.05), ([0.25, 0.13], 0.06)] |
| 261 ) | 261 # ) |
| 262 dcμ_est = (dat.apply((μ_perturb, auxtrue)) - dat.apply((μ_true, auxtrue))) / delta | 262 # dcμ_est = (dat.apply((μ_perturb, auxtrue)) - dat.apply((μ_true, auxtrue))) / delta |
| 263 print( | 263 # print( |
| 264 "μ-diff-test:", | 264 # "μ-diff-test:", |
| 265 eval_Function_f64(ω_true._cpp_object, [0.1, 0.3, 0.0]), | 265 # eval_Function_f64(ω_true._cpp_object, [0.1, 0.3, 0.0]), |
| 266 ", ", | 266 # ", ", |
| 267 dcμ_est, | 267 # dcμ_est, |
| 268 ) | 268 # ) |
| 269 | 269 |
| 270 return dat, auxtrue, aux_noisy, μ_bound, μ_true, plot_factory, pde | 270 return dat, auxtrue, aux_noisy, μ_bound, μ_true, plot_factory, pde |
| 271 | 271 |
| 272 | 272 |
| 273 def relnoise(v, σ, rng, level=None): | 273 def relnoise(v, σ, rng, level=None): |
