AgeColumnSystem init m_A m_A_e m_A_n m_A_s m_A_w m_age3 m_nu solve

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().