| 25 μ0 = DiscreteMeasure_2_f64([]) |
25 μ0 = DiscreteMeasure_2_f64([]) |
| 26 |
26 |
| 27 reg = RegTerm.NonnegRadon(1.5e-6) |
27 reg = RegTerm.NonnegRadon(1.5e-6) |
| 28 |
28 |
| 29 dat_simple = ComposeFnWithOperator(dat, InjectSecond(auxtrue)) |
29 dat_simple = ComposeFnWithOperator(dat, InjectSecond(auxtrue)) |
| 30 dat_simple.curvature_bound_components = lambda: (None, None) |
30 |
| 31 print("Initial data term value:", dat_simple.apply(μ0)) |
31 print("Initial data term value:", dat_simple.apply(μ0)) |
| 32 print("Data term value at true μ:", dat_simple.apply(μ_true)) |
32 print("Data term value at true μ:", dat_simple.apply(μ_true)) |
| 33 |
33 |
| 34 return Problem( |
34 return Problem( |
| 35 dataterm=dat_simple, regterm=reg, μinit=μ0, plot_factory=plot_factory |
35 dataterm=dat_simple, regterm=reg, μinit=μ0, plot_factory=plot_factory |