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

◆ solve()

void pism::AgeColumnSystem::solve ( std::vector< double > &  x)

First-order upwind scheme with implicit in the vertical: one column solve.

The PDE being solved is

\[ \frac{\partial \tau}{\partial t} + \frac{\partial}{\partial x}\left(u \tau\right) + \frac{\partial}{\partial y}\left(v \tau\right) + \frac{\partial}{\partial z}\left(w \tau\right) = 1. \]

Definition at line 68 of file AgeColumnSystem.cc.

References pism::RuntimeError::add_context(), pism::k, m_A, m_A_e, m_A_n, m_A_s, m_A_w, pism::columnSystemCtx::m_dt, pism::columnSystemCtx::m_dx, pism::columnSystemCtx::m_dy, pism::columnSystemCtx::m_i, pism::columnSystemCtx::m_j, pism::columnSystemCtx::m_ks, m_nu, pism::columnSystemCtx::m_solver, pism::columnSystemCtx::m_u, pism::columnSystemCtx::m_v, pism::columnSystemCtx::m_w, pism::columnSystemCtx::reportColumnZeroPivotErrorMFile(), and S().

Referenced by pism::AgeModel::update().