Mercurial > repos > pointsource_pde / file diff

--- a/experiments/laser_and_mirrors.py	Thu Feb 26 09:32:12 2026 -0500
+++ b/experiments/laser_and_mirrors.py	Wed Apr 22 22:32:00 2026 -0500
@@ -18,20 +18,18 @@
 
 # Setup routine
 def setup(prefix):
-    dat, auxtrue, μ_bound, μ_true, plot_factory = laser_and_mirrors_aux.generic_setup(
+    dat, auxtrue, _μ_bound, μ_true, plot_factory = laser_and_mirrors_aux.generic_setup(
         prefix
     )
 
     μ0 = DiscreteMeasure_2_f64([])
 
-    reg = RegTerm.NonnegRadon(
-        10e-7  # 5e-7  # 3e-7
-    )
+    reg = RegTerm.NonnegRadon(1.5e-6)
 
-    dat_simple = ComposeFnWithOperator(dat, InjectSecond(auxtrue), xbound=μ_bound)
+    dat_simple = ComposeFnWithOperator(dat, InjectSecond(auxtrue))
     dat_simple.curvature_bound_components = lambda: (None, None)
-    print("Initial value:", dat_simple.apply(μ0))
-    print("Value at true μ:", dat_simple.apply(μ_true))
+    print("Initial data term value:", dat_simple.apply(μ0))
+    print("Data term value at true μ:", dat_simple.apply(μ_true))
 
     return Problem(
         dataterm=dat_simple, regterm=reg, μinit=μ0, plot_factory=plot_factory