aphreco 0.1.3

A small-scale physiological modeling library
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62

use super::genetic_algorithm::GeneticAlgorithm;
use super::neldermead::NelderMead;
use crate::model::OptModelTrait;
use crate::objective::Objective;
use crate::optimizer::OptResult;

pub enum OptOptions {
  Default,

  NelderMead {
    max_iter: u64,
    adaptive: bool,
    x_abstol: f64,
    f_abstol: f64,
    verbose: bool,
  },

  GeneticAlgorithm {
    max_gen: u64,
    n_pop: usize,
    mutation_rate: f64,
    verbose: bool,
  },
}

pub enum Optimizer {
  NelderMead(OptOptions),
  GeneticAlgorithm(OptOptions),
}

impl Optimizer {
  pub fn run<M, const LEN_Y: usize, const LEN_P: usize, const LEN_B: usize, const LEN_X: usize>(
    &self,
    objective: &mut Objective<M, LEN_Y, LEN_P, LEN_B, LEN_X>,
  ) -> OptResult
  where
    M: OptModelTrait<LEN_Y, LEN_P, LEN_B, LEN_X>,
  {
    match self {
      Optimizer::NelderMead(options) => {
        let opt = NelderMead::new(objective.len_x, options);
        opt.run(objective)
      }

      Optimizer::GeneticAlgorithm(options) => {
        let opt = GeneticAlgorithm::new(objective.len_x, options);
        opt.run(objective)
      }
    }
  }
}

pub trait ConcreteOptimizer {
  fn new(len_x: usize, options: &OptOptions) -> Self;

  fn run<M, const LEN_Y: usize, const LEN_P: usize, const LEN_B: usize, const LEN_X: usize>(
    &self,
    objective: &mut Objective<M, LEN_Y, LEN_P, LEN_B, LEN_X>,
  ) -> OptResult
  where
    M: OptModelTrait<LEN_Y, LEN_P, LEN_B, LEN_X>;
}