#![allow(non_snake_case)]
use crate::algebra::*;
pub struct DefaultEquilibrationData<T> {
pub d: Vec<T>,
pub dinv: Vec<T>,
pub e: Vec<T>,
pub einv: Vec<T>,
pub c: T,
}
impl<T> DefaultEquilibrationData<T>
where
T: FloatT,
{
pub fn new(n: usize, m: usize) -> Self {
let d = vec![T::one(); n];
let dinv = vec![T::one(); n];
let e = vec![T::one(); m];
let einv = vec![T::one(); m];
let c = T::one();
Self {
d,
dinv,
e,
einv,
c,
}
}
}