Struct ConstantDerivativeMatrixProof
pub struct ConstantDerivativeMatrixProof<C = f64> { /* private fields */ }Description
A constant derivative Jacobian whose rank is computed without a numerical tolerance.
Every finite f64 coefficient is interpreted as the exact binary rational
number represented by its bits. Rank is then recomputed with arbitrary-
precision rational elimination. The stored rank is therefore evidence
about the lowered matrix itself, not a sample-state estimate or a
backend-dependent floating-point classification. The binary64 constructor
retains binary64 coefficients for numerical time backends; from_exact
retains rational coefficients for symbolic compatibility projection. Both
representations share the same exact elimination and rank owner.
Implementations§
§impl ConstantDerivativeMatrixProof
impl ConstantDerivativeMatrixProof
pub fn new(
dimension: usize,
coefficients: Vec<f64>,
) -> Result<ConstantDerivativeMatrixProof, Diagnostic>
pub fn new( dimension: usize, coefficients: Vec<f64>, ) -> Result<ConstantDerivativeMatrixProof, Diagnostic>
Construct and exactly classify one square constant derivative matrix.
§Errors
Returns EQ0705 for an empty/non-square matrix or a non-finite
coefficient.
§impl ConstantDerivativeMatrixProof<Ratio<BigInt>>
impl ConstantDerivativeMatrixProof<Ratio<BigInt>>
pub fn from_exact(
dimension: usize,
coefficients: Vec<Ratio<BigInt>>,
) -> Result<ConstantDerivativeMatrixProof<Ratio<BigInt>>, Diagnostic>
pub fn from_exact( dimension: usize, coefficients: Vec<Ratio<BigInt>>, ) -> Result<ConstantDerivativeMatrixProof<Ratio<BigInt>>, Diagnostic>
Construct a matrix from exact symbolic coefficients without first rounding their sums or products to binary64.
§Errors
Rejects empty/non-square matrices and raw rationals with zero denominators.
§impl<C> ConstantDerivativeMatrixProof<C>
impl<C> ConstantDerivativeMatrixProof<C>
pub const fn dimension(&self) -> usize
pub const fn dimension(&self) -> usize
Number of state coordinates and residual rows.
pub fn coefficients(&self) -> &[C]
pub fn coefficients(&self) -> &[C]
Complete row-major coefficient storage.
pub fn row(&self, row: usize) -> Option<&[C]>
pub fn row(&self, row: usize) -> Option<&[C]>
One residual row in state-coordinate order.
pub const fn exact_rank(&self) -> usize
pub const fn exact_rank(&self) -> usize
Exact rank over the binary rational values represented by the matrix.
pub fn compatibility_residual(
&self,
additional_columns: &[Vec<Ratio<BigInt>>],
rhs: &[Ratio<BigInt>],
) -> Result<Vec<f64>, Diagnostic>
pub fn compatibility_residual( &self, additional_columns: &[Vec<Ratio<BigInt>>], rhs: &[Ratio<BigInt>], ) -> Result<Vec<f64>, Diagnostic>
Orthogonal residual of M a + C u = rhs in the original row coordinates.
additional_columns supplies C column by column. Elimination and the
least-squares projection use exact binary rationals; only the returned
residual is rounded to binary64. The RHS is already exact so coefficient
accumulation cannot erase a constrained direction before projection.
Callers own residual scaling/tolerance.
This does not decide infinity-norm tolerance feasibility.
§Errors
Returns EQ0705 for invalid dimensions, nonfinite inputs, or a residual
that cannot be represented as a finite binary64 value.
pub fn is_compatible(
&self,
additional_columns: &[Vec<Ratio<BigInt>>],
column: &[Ratio<BigInt>],
) -> Result<bool, Diagnostic>
pub fn is_compatible( &self, additional_columns: &[Vec<Ratio<BigInt>>], column: &[Ratio<BigInt>], ) -> Result<bool, Diagnostic>
Whether a column belongs exactly to the span of M and the additional columns, before rounding any residual. This is not a tolerance test.
§Errors
Returns EQ0705 for nonfinite inputs or mismatched row dimensions.
pub fn require_index_one_regularity(
&self,
state_jacobian: &[f64],
) -> Result<(), Diagnostic>
pub fn require_index_one_regularity( &self, state_jacobian: &[f64], ) -> Result<(), Diagnostic>
Require local index-one regularity of the constant-mass residual.
state_jacobian is the row-major Jacobian of the regular residuals
with respect to values, excluding initial equations and their tangents.
This checks one point; it performs no index reduction or state selection.
§Errors
Returns EQ0705 for an invalid Jacobian or a singular/high-index block.
§impl ConstantDerivativeMatrixProof
impl ConstantDerivativeMatrixProof
pub fn monomial_rows(&self) -> Option<Vec<MonomialDerivativeRow>>
pub fn monomial_rows(&self) -> Option<Vec<MonomialDerivativeRow>>
Derive the monomial row view used only for explicit-ODE normalization.
Returns None unless every row has exactly one non-zero coefficient
and every state coordinate occurs exactly once.
Trait Implementations§
§impl<C> Clone for ConstantDerivativeMatrixProof<C>where
C: Clone,
impl<C> Clone for ConstantDerivativeMatrixProof<C>where
C: Clone,
§impl<C> Debug for ConstantDerivativeMatrixProof<C>where
C: Debug,
impl<C> Debug for ConstantDerivativeMatrixProof<C>where
C: Debug,
§impl<C> PartialEq for ConstantDerivativeMatrixProof<C>where
C: PartialEq,
impl<C> PartialEq for ConstantDerivativeMatrixProof<C>where
C: PartialEq,
§fn eq(&self, other: &ConstantDerivativeMatrixProof<C>) -> bool
fn eq(&self, other: &ConstantDerivativeMatrixProof<C>) -> bool
self and other values to be equal, and is used by ==.impl<C> StructuralPartialEq for ConstantDerivativeMatrixProof<C>
Auto Trait Implementations§
impl<C> Freeze for ConstantDerivativeMatrixProof<C>
impl<C> RefUnwindSafe for ConstantDerivativeMatrixProof<C>where
C: RefUnwindSafe,
impl<C> Send for ConstantDerivativeMatrixProof<C>where
C: Send,
impl<C> Sync for ConstantDerivativeMatrixProof<C>where
C: Sync,
impl<C> Unpin for ConstantDerivativeMatrixProof<C>where
C: Unpin,
impl<C> UnsafeUnpin for ConstantDerivativeMatrixProof<C>
impl<C> UnwindSafe for ConstantDerivativeMatrixProof<C>where
C: UnwindSafe,
Blanket Implementations§
Source§impl<T> Any for Twhere
T: 'static + ?Sized,
impl<T> Any for Twhere
T: 'static + ?Sized,
§impl<Src, Scheme> ApproxFrom<Src, Scheme> for Srcwhere
Scheme: ApproxScheme,
impl<Src, Scheme> ApproxFrom<Src, Scheme> for Srcwhere
Scheme: ApproxScheme,
§fn approx_from(src: Src) -> Result<Src, <Src as ApproxFrom<Src, Scheme>>::Err>
fn approx_from(src: Src) -> Result<Src, <Src as ApproxFrom<Src, Scheme>>::Err>
§impl<Dst, Src, Scheme> ApproxInto<Dst, Scheme> for Srcwhere
Dst: ApproxFrom<Src, Scheme>,
Scheme: ApproxScheme,
impl<Dst, Src, Scheme> ApproxInto<Dst, Scheme> for Srcwhere
Dst: ApproxFrom<Src, Scheme>,
Scheme: ApproxScheme,
§type Err = <Dst as ApproxFrom<Src, Scheme>>::Err
type Err = <Dst as ApproxFrom<Src, Scheme>>::Err
§fn approx_into(self) -> Result<Dst, <Src as ApproxInto<Dst, Scheme>>::Err>
fn approx_into(self) -> Result<Dst, <Src as ApproxInto<Dst, Scheme>>::Err>
Source§impl<T> Borrow<T> for Twhere
T: ?Sized,
impl<T> Borrow<T> for Twhere
T: ?Sized,
Source§impl<T> BorrowMut<T> for Twhere
T: ?Sized,
impl<T> BorrowMut<T> for Twhere
T: ?Sized,
Source§fn borrow_mut(&mut self) -> &mut T
fn borrow_mut(&mut self) -> &mut T
Source§impl<T> CloneToUninit for Twhere
T: Clone,
impl<T> CloneToUninit for Twhere
T: Clone,
§impl<T, Dst> ConvAsUtil<Dst> for T
impl<T, Dst> ConvAsUtil<Dst> for T
§fn approx(self) -> Result<Dst, Self::Err>where
Self: Sized + ApproxInto<Dst>,
fn approx(self) -> Result<Dst, Self::Err>where
Self: Sized + ApproxInto<Dst>,
§impl<T> ConvUtil for T
impl<T> ConvUtil for T
§fn approx_as<Dst>(self) -> Result<Dst, Self::Err>where
Self: Sized + ApproxInto<Dst>,
fn approx_as<Dst>(self) -> Result<Dst, Self::Err>where
Self: Sized + ApproxInto<Dst>,
§fn approx_as_by<Dst, Scheme>(self) -> Result<Dst, Self::Err>where
Self: Sized + ApproxInto<Dst, Scheme>,
Scheme: ApproxScheme,
fn approx_as_by<Dst, Scheme>(self) -> Result<Dst, Self::Err>where
Self: Sized + ApproxInto<Dst, Scheme>,
Scheme: ApproxScheme,
§fn into_as<Dst>(self) -> Dstwhere
Self: Sized + Into<Dst>,
fn into_as<Dst>(self) -> Dstwhere
Self: Sized + Into<Dst>,
§fn try_as<Dst>(self) -> Result<Dst, Self::Err>where
Self: Sized + TryInto<Dst>,
fn try_as<Dst>(self) -> Result<Dst, Self::Err>where
Self: Sized + TryInto<Dst>,
§impl<T> DistributionExt for Twhere
T: ?Sized,
impl<T> DistributionExt for Twhere
T: ?Sized,
fn rand<T>(&self, rng: &mut (impl Rng + ?Sized)) -> Twhere
Self: Distribution<T>,
Source§impl<T> From<T> for T
impl<T> From<T> for T
§impl<T> Instrument for T
impl<T> Instrument for T
§fn instrument(self, span: Span) -> Instrumented<Self>
fn instrument(self, span: Span) -> Instrumented<Self>
§fn in_current_span(self) -> Instrumented<Self>
fn in_current_span(self) -> Instrumented<Self>
Source§impl<T, U> Into<U> for Twhere
U: From<T>,
impl<T, U> Into<U> for Twhere
U: From<T>,
Source§impl<T> IntoEither for T
impl<T> IntoEither for T
Source§fn into_either(self, into_left: bool) -> Either<Self, Self>
fn into_either(self, into_left: bool) -> Either<Self, Self>
self into a Left variant of Either<Self, Self>
if into_left is true.
Converts self into a Right variant of Either<Self, Self>
otherwise. Read moreSource§fn into_either_with<F>(self, into_left: F) -> Either<Self, Self>where
F: FnOnce(&Self) -> bool,
fn into_either_with<F>(self, into_left: F) -> Either<Self, Self>where
F: FnOnce(&Self) -> bool,
self into a Left variant of Either<Self, Self>
if into_left(&self) returns true.
Converts self into a Right variant of Either<Self, Self>
otherwise. Read more§impl<T> Pointable for T
impl<T> Pointable for T
Source§impl<T> Same for T
impl<T> Same for T
§impl<SS, SP> SupersetOf<SS> for SPwhere
SS: SubsetOf<SP>,
impl<SS, SP> SupersetOf<SS> for SPwhere
SS: SubsetOf<SP>,
§fn to_subset(&self) -> Option<SS>
fn to_subset(&self) -> Option<SS>
self from the equivalent element of its
superset. Read more§fn is_in_subset(&self) -> bool
fn is_in_subset(&self) -> bool
self is actually part of its subset T (and can be converted to it).§fn to_subset_unchecked(&self) -> SS
fn to_subset_unchecked(&self) -> SS
self.to_subset but without any property checks. Always succeeds.§fn from_subset(element: &SS) -> SP
fn from_subset(element: &SS) -> SP
self to the equivalent element of its superset.