assemble_matrix compute_nuH_norm init_impl m_A m_default_pc_failure_count m_default_pc_failure_max_count m_KSP m_nuH_old m_nuh_viewer m_regional_mode m_velocity_old pc_setup_asm pc_setup_bjacobi picard_iteration picard_manager picard_strategy_regularization solve SSAFD update_nuH_viewers write_system_petsc ~SSAFD

SSAFD() pism::stressbalance::SSAFD::SSAFD ( std::shared_ptr< const Grid >  grid, bool  regional_mode  ) Because the FD implementation of the SSA uses Picard iteration, a PETSc KSP and Mat are used directly. In particular we set up \(A\) (Mat m_A) and a \(b\) (= Vec m_rhs) and iteratively solve linear systems \[ A x = b \] where \(x\) (= Vec m_velocity_global). A PETSc SNES object is never created. Definition at line 54 of file SSAFD.cc. References pism::array::Array::dm(), pism::VariableMetadata::long_name(), m_A, pism::Component::m_grid, m_KSP, m_nuH_old, m_nuh_viewer, pism::stressbalance::SSA::m_velocity_global, m_velocity_old, pism::array::Array::metadata(), PISM_CHK, pism::Wrapper< T >::rawptr(), and pism::VariableMetadata::units().