Mercurial > repos > pointsource_pde / file comparison
comparison: experiments/laser_and_mirrors_aux.py
experiments/laser_and_mirrors_aux.py
- changeset 10
- 0c9e19020217
- parent 3
- c3a4f4bb87f7
- child 11
- c4f923017349
equal
deleted
inserted
replaced
| 9:a5e6fac08015 | 10:0c9e19020217 |
|---|---|
| 150 for data in b | 150 for data in b |
| 151 ], | 151 ], |
| 152 ), | 152 ), |
| 153 pde, | 153 pde, |
| 154 ) | 154 ) |
| 155 dat.curvature_bound_components = lambda: (None, None, None) | |
| 155 | 156 |
| 156 # Save true and initial data | 157 # Save true and initial data |
| 157 plot_factory = PlotFactory(pde) | 158 plot_factory = PlotFactory(pde) |
| 158 plotter = plot_factory.produce(prefix) | 159 plotter = plot_factory.produce(prefix) |
| 159 ω, _ = dat.diff((DiscreteMeasure_2_f64([]), auxtrue)) | 160 ω, _ = dat.diff((DiscreteMeasure_2_f64([]), auxtrue)) |
| 207 ax = aux.apply(aux0) | 208 ax = aux.apply(aux0) |
| 208 inival = dat.apply((μ0, aux0)) | 209 inival = dat.apply((μ0, aux0)) |
| 209 print("Initial data term value:", inival + ax) | 210 print("Initial data term value:", inival + ax) |
| 210 print("Data term value at true μ:", dat.apply((μ_true, auxtrue)) + ax) | 211 print("Data term value at true μ:", dat.apply((μ_true, auxtrue)) + ax) |
| 211 | 212 |
| 212 # No curvature bound given: θ is aboslute. | |
| 213 dat.curvature_bound_components = lambda: (None, None) | |
| 214 | |
| 215 return Problem(dat, reg, aux, aux0, μ0, plot_factory=plot_factory) | 213 return Problem(dat, reg, aux, aux0, μ0, plot_factory=plot_factory) |
