D ddratio L m_D m_L m_max_system_size m_prefix m_rhs m_U m_work norm1 prefix reset RHS save_matrix save_system save_system_with_solution save_vector solve solve TridiagonalSystem U

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.