Minor changes to algorithm dispatch and choice

src/algorithms/map.rs

@@ -15,6 +15,7 @@ use super::{initialization, output::CycleLog, NonParametricAlgorithm};
 /// Maximum a posteriori (MAP) estimation
 ///
 /// Calculate the MAP estimate of the parameters of the model given the data.
+#[derive(Debug, Clone)]
 pub struct MAP<E: Equation> {
     equation: E,
     psi: Array2<f64>,

src/algorithms/mod.rs

@@ -3,7 +3,7 @@ use std::path::Path;
 
 use crate::prelude::{self, settings::Settings};
 
-use anyhow::{bail, Result};
+use anyhow::Result;
 use anyhow::{Context, Error};
 use map::MAP;
 use ndarray::Array2;
@@ -23,26 +23,14 @@ pub mod routines;
 /// Supported algorithms by `PMcore`
 #[derive(Debug, Clone, Serialize, Deserialize, PartialEq, PartialOrd)]
 pub enum Algorithm {
-    NonParametric(NonParametric),
-    Parametric(Parametric),
-}
-
-/// Supported non-parametric algorithms
-#[derive(Debug, Clone, Serialize, Deserialize, PartialEq, PartialOrd)]
-pub enum NonParametric {
+    // Non-parametric algorithms
     NPAG,
     NPOD,
     MAP,
+    // Parametric algorithms
 }
 
-/// Supported parametric algorithms
-#[derive(Debug, Clone, Serialize, Deserialize, PartialEq, PartialOrd)]
-pub enum Parametric {
-    FOCE,
-    NPSA,
-}
-
-/// This traint defines the methods for non-parametric (NP) algorithms
+/// This trait defines the methods for non-parametric (NP) algorithms
 pub trait NonParametricAlgorithm<E: Equation> {
     fn new(config: Settings, equation: E, data: Data) -> Result<Box<Self>, Error>
     where
@@ -122,15 +110,20 @@ pub trait NonParametricAlgorithm<E: Equation> {
     fn into_npresult(&self) -> NPResult<E>;
 }
 
+pub trait ParametricAlgorithm<E: Equation> {
+    fn fit(&mut self) -> Result<()> {
+        unimplemented!()
+    }
+}
+
 pub fn dispatch_algorithm<E: Equation>(
     settings: Settings,
     equation: E,
     data: Data,
 ) -> Result<Box<dyn NonParametricAlgorithm<E>>> {
     match settings.config().algorithm {
-        Algorithm::NonParametric(NonParametric::NPAG) => Ok(NPAG::new(settings, equation, data)?),
-        Algorithm::NonParametric(NonParametric::NPOD) => Ok(NPOD::new(settings, equation, data)?),
-        Algorithm::NonParametric(NonParametric::MAP) => Ok(MAP::new(settings, equation, data)?),
-        _ => bail!("Unsupported algorithm"),
+        Algorithm::NPAG => Ok(NPAG::new(settings, equation, data)?),
+        Algorithm::NPOD => Ok(NPOD::new(settings, equation, data)?),
+        Algorithm::MAP => Ok(MAP::new(settings, equation, data)?),
     }
 }

src/algorithms/npag.rs

@@ -24,7 +24,7 @@ const THETA_G: f64 = 1e-4; // Objective function convergence criteria
 const THETA_F: f64 = 1e-2;
 const THETA_D: f64 = 1e-4;
 
-#[derive(Debug)]
+#[derive(Debug, Clone)]
 pub struct NPAG<E: Equation> {
     equation: E,
     psi: Array2<f64>,

src/algorithms/npod.rs

@@ -28,6 +28,7 @@ use super::{
 const THETA_F: f64 = 1e-2;
 const THETA_D: f64 = 1e-4;
 
+#[derive(Debug, Clone)]
 pub struct NPOD<E: Equation> {
     equation: E,
     psi: Array2<f64>,

src/algorithms/routines/settings.rs

@@ -227,7 +227,7 @@ impl Default for Config {
     fn default() -> Self {
         Config {
             cycles: 100,
-            algorithm: Algorithm::NonParametric(crate::algorithms::NonParametric::NPAG),
+            algorithm: Algorithm::NPAG,
             cache: true,
         }
     }
@@ -793,7 +793,7 @@ impl SettingsBuilder<ErrorSet> {
 
 mod tests {
     use super::*;
-    use crate::algorithms::{Algorithm, NonParametric};
+    use crate::algorithms::Algorithm;
     use pharmsol::prelude::data::ErrorType;
 
     #[test]
@@ -805,7 +805,7 @@ mod tests {
             .unwrap();
 
         let settings = SettingsBuilder::new()
-            .set_algorithm(Algorithm::NonParametric(NonParametric::NPAG)) // Step 1: Define algorithm
+            .set_algorithm(Algorithm::NPAG) // Step 1: Define algorithm
             .set_parameters(parameters) // Step 2: Define parameters
             .set_error_model(Error {
                 value: 0.1,
@@ -814,9 +814,6 @@ mod tests {
             }) // Step 3: Define error model
             .build(); // Final step
 
-        assert_eq!(
-            settings.config.algorithm,
-            Algorithm::NonParametric(NonParametric::NPAG,)
-        );
+        assert_eq!(settings.config.algorithm, Algorithm::NPAG);
     }
 }