polycal 0.1.10

methods for fitting and using polynomial calibration functions following ISO/TS 28038
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

use super::{Covariance, Solution, SolveSystem, SolverError};
use ndarray::{Array1, Array2, ArrayView1, ArrayView2, ScalarOperand};
use ndarray_linalg::{Lapack, Scalar};

pub struct TotalLeastSquares<'a, E> {
    pub(crate) y: Array1<E>,
    pub(crate) covariance_x: Covariance<'a, E>,
    pub(crate) covariance_y: Covariance<'a, E>,
    pub(crate) h: Array2<E>,
}

impl<E: Lapack + Scalar<Real = E> + ScalarOperand> SolveSystem<E> for TotalLeastSquares<'_, E> {
    fn solve(&self) -> Result<Solution<E>, SolverError> {
        match (&self.covariance_x, &self.covariance_y) {
            (Covariance::Diagonal(ux), Covariance::Diagonal(uy)) => self.solve_diagonal(*ux, *uy),
            (Covariance::Matrix(vx), Covariance::Matrix(vy)) => self.solve_full_rank(*vx, *vy),
            _ => unreachable!("make this inaccessible via type."),
        }
    }
}

impl<'a, E: Lapack + Scalar<Real = E> + ScalarOperand> TotalLeastSquares<'a, E> {
    #[allow(clippy::unused_self)]
    fn solve_diagonal(
        &self,
        _ux: ArrayView1<'a, E>,
        _uy: ArrayView1<'a, E>,
    ) -> Result<Solution<E>, SolverError> {
        unimplemented!("no diagonal TLS impl");
    }

    #[allow(clippy::unused_self)]
    fn solve_full_rank(
        &self,
        _vx: ArrayView2<'a, E>,
        _vy: ArrayView2<'a, E>,
    ) -> Result<Solution<E>, SolverError> {
        unimplemented!("no full TLS impl");
    }
}