Module mass_matrix 

Utilities for working with mass matrices in ODE systems

This module provides functions for handling mass matrices in ODEs of the form: M(t,y)·y’ = f(t,y), where M is a mass matrix that may depend on time t and state y.

Functions

apply_mass

Apply mass matrix to a vector: result = M·v

check_mass_compatibility

Check if a mass matrix is compatible with an ODE state

is_singular

Check if a matrix is singular or ill-conditioned

solve_mass_system

Solve a linear system with mass matrix: M·x = b

transform_to_standard_form

Transform standard ODE to form with identity mass matrix