L2 distance can now compute best model fit

specific/data/objective/l2distance.cpp

@@ -45,6 +45,20 @@ namespace l2distance
         }
         return _v;
     }
+
+    template <Physic T_Physic>
+    std::vector < std::vector < std::complex<double>>> Objective<T_Physic>::improvedShotIntensity(const Data<T_Physic> &d1) const
+    {
+        throw Exception("Shot fitting only works in Physic::acoustic so far");
+        //return leastSquareShotIntensity(d1, d0);
+    }
+
+    template <>
+    std::vector<std::vector<std::complex<double>>> Objective<Physic::acoustic>::improvedShotIntensity(const Data<Physic::acoustic> &d1) const
+    {
+        return leastSquareShotIntensity(d1, d0);
+    }
+
     template class Objective<Physic::acoustic>;
     template class Objective<Physic::electromagnetic>;
     template class Objective<Physic::elastodynamic>;

specific/data/objective/l2distance.h

@@ -32,6 +32,8 @@ namespace l2distance
 
         virtual double performance(const Data<T_Physic>& d);
         virtual const Data<T_Physic>& update(Type type, const DataStateEvaluator<T_Physic>& ds);
+
+        std::vector<std::vector<std::complex<double>>> improvedShotIntensity(const Data<T_Physic>& d) const;
     };
 };