diff -r c5d8bd1a7728 -r 0f59c0d02e13 src/run.rs
--- a/src/run.rs	Thu Jan 23 23:35:28 2025 +0100
+++ b/src/run.rs	Mon Jan 06 21:37:03 2025 -0500
@@ -327,7 +327,7 @@
 // }
 
 /// Plotting level configuration
-#[derive(Copy, Clone, Eq, PartialEq, Ord, PartialOrd, Serialize, ValueEnum, Debug)]
+#[derive(Copy, Clone, Eq, PartialEq, Ord, PartialOrd, Serialize, Deserialize, ValueEnum, Debug)]
 pub enum PlotLevel {
     /// Plot nothing
     #[clap(name = "none")]
@@ -357,7 +357,7 @@
 >;
 
 /// This is a dirty workaround to rust-csv not supporting struct flattening etc.
-#[derive(Serialize)]
+#[derive(Serialize, Deserialize)]
 struct CSVLog<F> {
     iter : usize,
     cpu_time : f64,
@@ -372,7 +372,7 @@
 }
 
 /// Collected experiment statistics
-#[derive(Clone, Debug, Serialize)]
+#[derive(Clone, Debug, Serialize, Deserialize)]
 struct ExperimentStats<F : Float> {
     /// Signal-to-noise ratio in decibels
     ssnr : F,
@@ -398,7 +398,7 @@
     }
 }
 /// Collected algorithm statistics
-#[derive(Clone, Debug, Serialize)]
+#[derive(Clone, Debug, Serialize, Deserialize)]
 struct AlgorithmStats<F : Float> {
     /// Overall CPU time spent
     cpu_time : F,
@@ -416,10 +416,29 @@
 
 
 /// Struct for experiment configurations
-#[derive(Debug, Clone, Serialize)]
+#[derive(Debug, Clone, Serialize, Deserialize)]
+#[serde(bound(
+    serialize = "Cube<F, N> : Serialize,
+                 NoiseDistr : Serialize,
+                 [usize; N] : Serialize,
+                 RNDM<F, N> : Serialize,
+                 Regularisation<F> : Serialize,
+                 F : Serialize,
+                 S : Serialize,
+                 P : Serialize,
+                 K : Serialize",
+    deserialize = "Cube<F, N> : for<'a> Deserialize<'a>,
+                   NoiseDistr : for<'a> Deserialize<'a>,
+                   [usize; N] : for<'a> Deserialize<'a>,
+                   RNDM<F, N> : for<'a> Deserialize<'a>,
+                   Regularisation<F> : for<'a> Deserialize<'a>,
+                   F : for<'a> Deserialize<'a>,
+                   S : for<'a> Deserialize<'a>,
+                   P : for<'a> Deserialize<'a>,
+                   K : for<'a> Deserialize<'a>,",
+))]
 pub struct ExperimentV2<F, NoiseDistr, S, K, P, const N : usize>
 where F : Float,
-      [usize; N] : Serialize,
       NoiseDistr : Distribution<F>,
       S : Sensor<F, N>,
       P : Spread<F, N>,
@@ -452,7 +471,15 @@
     pub algorithm_defaults : HashMap<DefaultAlgorithm, AlgorithmConfig<F>>,
 }
 
-#[derive(Debug, Clone, Serialize)]
+#[derive(Debug, Clone, Serialize, Deserialize)]
+#[serde(bound(
+    serialize = "ExperimentV2<F, NoiseDistr, S, K, P, N> : Serialize,
+                 B : Serialize,
+                 F : Serialize",
+    deserialize = "ExperimentV2<F, NoiseDistr, S, K, P, N >: for<'a> Deserialize<'a>,
+                   B : for<'a> Deserialize<'a>,
+                   F : for<'a> Deserialize<'a>",
+))]
 pub struct ExperimentBiased<F, NoiseDistr, S, K, P, B, const N : usize>
 where F : Float,
       [usize; N] : Serialize,
@@ -541,6 +568,7 @@
 ) ->  DynError
 where
     PlotLookup : Plotting<N>,
+    DeltaMeasure<Loc<F, N>, F> : Serialize,
 {
     let mut logs = Vec::new();
 
@@ -647,6 +675,8 @@
     BTNodeLookup: BTNode<F, usize, Bounds<F>, N>,
     RNDM<F, N> : SpikeMerging<F>,
     NoiseDistr : Distribution<F> + Serialize + std::fmt::Debug,
+    [F; N] : Serialize,
+    [[F; 2]; N] : Serialize,
     // DefaultSG<F, S, P, N> : ForwardModel<RNDM<F, N>, F, PreadjointCodomain = PreadjointCodomain, Observable=DVector<F::MixedType>>,
     // PreadjointCodomain : Space + Bounded<F> + DifferentiableRealMapping<F, N>,
     // DefaultSeminormOp<F, K, N> : ProxPenalty<F, PreadjointCodomain, RadonRegTerm<F>, N>,
@@ -912,6 +942,8 @@
     RNDM<F, N> : SpikeMerging<F>,
     NoiseDistr : Distribution<F> + Serialize + std::fmt::Debug,
     B : Mapping<Loc<F, N>, Codomain = F> + Serialize + std::fmt::Debug,
+    [F; N] : Serialize,
+    [[F; 2]; N] : Serialize,
     // DefaultSG<F, S, P, N> : ForwardModel<RNDM<F, N>, F, PreadjointCodomain = PreadjointCodomain, Observable=DVector<F::MixedType>>,
     // PreadjointCodomain : Bounded<F> + DifferentiableRealMapping<F, N>,
     // DefaultSeminormOp<F, K, N> :  ProxPenalty<F, PreadjointCodomain, RadonRegTerm<F>, N>,
@@ -1169,18 +1201,21 @@
     b : &A::Observable,
     kernel_plot_width : F,
 ) -> DynError
-where F : Float + ToNalgebraRealField,
-      Sensor : RealMapping<F, N> + Support<F, N> + Clone,
-      Spread : RealMapping<F, N> + Support<F, N> + Clone,
-      Kernel : RealMapping<F, N> + Support<F, N>,
-      Convolution<Sensor, Spread> : DifferentiableRealMapping<F, N> + Support<F, N>,
-      𝒟 : DiscreteMeasureOp<Loc<F, N>, F>,
-      𝒟::Codomain : RealMapping<F, N>,
-      A : ForwardModel<RNDM<F, N>, F>,
-      for<'a> &'a A::Observable : Instance<A::Observable>,
-      A::PreadjointCodomain : DifferentiableRealMapping<F, N> + Bounded<F>,
-      PlotLookup : Plotting<N>,
-      Cube<F, N> : SetOrd {
+where
+    F : Float + ToNalgebraRealField,
+    Sensor : RealMapping<F, N> + Support<F, N> + Clone,
+    Spread : RealMapping<F, N> + Support<F, N> + Clone,
+    Kernel : RealMapping<F, N> + Support<F, N>,
+    Convolution<Sensor, Spread> : DifferentiableRealMapping<F, N> + Support<F, N>,
+    𝒟 : DiscreteMeasureOp<Loc<F, N>, F>,
+    𝒟::Codomain : RealMapping<F, N>,
+    A : ForwardModel<RNDM<F, N>, F>,
+    for<'a> &'a A::Observable : Instance<A::Observable>,
+    A::PreadjointCodomain : DifferentiableRealMapping<F, N> + Bounded<F>,
+    PlotLookup : Plotting<N>,
+    Cube<F, N> : SetOrd,
+    DeltaMeasure<Loc<F, N>, F> : Serialize,
+{
 
     if cli.plot < PlotLevel::Data {
         return Ok(())