solve_mass_system

scirs2_integrate::ode::utils::mass_matrix

Function solve_mass_system 

Source 
pub fn solve_mass_system<F>(
    mass: &MassMatrix<F>,
    t: F,
    y: ArrayView1<'_, F>,
    b: ArrayView1<'_, F>,
) -> IntegrateResult<Array1<F>>
where
    F: IntegrateFloat,
Expand description

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

For ODEs with mass matrices, we often need to solve M·x = f(t,y) to find x = y’ for the standard form y’ = g(t,y)

§Arguments

  • mass - The mass matrix structure
  • t - Current time
  • y - Current state
  • b - Right-hand side vector

§Returns

Solution vector x where M·x = b, or error if the system cannot be solved