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

TridiagonalSystem() pism::TridiagonalSystem::TridiagonalSystem ( unsigned int  max_size, const std::string &  prefix  ) Allocate a tridiagonal system of maximum size max_system_size. Let N = max_system_size. Then allocated locations are like this: D[0] U[0] 0 0 0 ... L[1] D[1] U[1] 0 0 ... 0 L[2] D[2] U[2] 0 ... 0 0 L[3] D[3] U[3] ... with the last row 0 0 ... 0 L[N-1] D[N-1] Thus the index into the arrays L, D, U is always the row number. Definition at line 50 of file ColumnSystem.cc. References m_D, m_L, m_max_system_size, m_rhs, m_U, and m_work.