PISM, A Parallel Ice Sheet Model  stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800

◆ solve() [1/2]

void pism::TridiagonalSystem::solve ( unsigned int  system_size,
double *  result 
)

The actual code for solving a tridiagonal system.

This is modified slightly from a Numerical Recipes version.

Input size n is size of instance. Requires n <= TridiagonalSystem::m_max_system_size.

Solution of system in x.

Definition at line 215 of file ColumnSystem.cc.

References pism::RuntimeError::formatted(), pism::k, m_D, m_L, m_max_system_size, m_rhs, m_U, m_work, and PISM_ERROR_LOCATION.