Struct ParamEstProblem 

pub struct ParamEstProblem<S: Scalar> { /* private fields */ }

Builder for ODE parameter estimation problems.

Implementations

impl<S: Scalar> ParamEstProblem<S>

pub fn new(n_params: usize, n_states: usize) -> Self

Create a new parameter estimation problem.

• n_params: number of parameters to estimate.
• n_states: dimension of the ODE state vector.

pub fn model<F>(self, f: F) -> Self

where F: Fn(S, &[S], &mut [S], &[S]) + Send + Sync + 'static,

Set the ODE right-hand side: f(t, y, dydt, params).

pub fn initial_state(self, y0: Vec<S>) -> Self

Set the initial state y(t0).

pub fn params(self, p0: Vec<S>) -> Self

Set the initial parameter guess.

pub fn param_bounds(self, i: usize, bounds: (S, S)) -> Self

Set bounds for parameter i.

pub fn all_param_bounds(self, bounds: Vec<Option<(S, S)>>) -> Self

Set bounds for all parameters at once.

pub fn data(self, t_data: Vec<S>, y_data: Vec<S>) -> Self

Set measurement data.

y_data is flat row-major: y_data[i * n_observed + j] for time index i and observed state j.

pub fn observed(self, indices: Vec<usize>) -> Self

Specify which state indices are observed.

If not called, all states are assumed observed.

pub fn ode_solver(self, choice: OdeSolverChoice) -> Self

Choose the ODE solver for forward integrations.

pub fn ode_tolerances(self, rtol: S, atol: S) -> Self

Set ODE solver tolerances.

pub fn max_iter(self, n: usize) -> Self

Set maximum optimizer iterations.

pub fn solve(self) -> Result<ParamEstResult<S>, OcpError>

where S: SimpleEntity + Conjugate<Canonical = S> + ComplexField,

Run parameter estimation.