assemble_matrix assemble_rhs compute_hardav_staggered compute_nuH_norm compute_nuH_staggered compute_nuH_staggered_cfbc diagnostics_impl fracture_induced_softening init_impl integrated_viscosity is_marginal m_A m_b m_default_pc_failure_count m_default_pc_failure_max_count m_hardness m_KSP m_nuH m_nuH_old m_nuh_viewer m_nuh_viewer_size m_scaling m_velocity_old m_view_nuh m_work 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 ) 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_b) and iteratively solve linear systems \[ A x = b \] where \(x\) (= Vec SSAX). A PETSc SNES object is never created. Definition at line 61 of file SSAFD.cc. References pism::VariableMetadata::long_name(), m_A, pism::stressbalance::SSA::m_da, pism::stressbalance::ShallowStressBalance::m_flow_law, pism::Component::m_grid, m_hardness, m_KSP, m_nuH, m_nuH_old, m_nuh_viewer_size, m_velocity_old, m_view_nuh, m_work, pism::array::Array::metadata(), PISM_CHK, pism::printf(), pism::Wrapper< T >::rawptr(), pism::VariableMetadata::set_units_without_validation(), and pism::VariableMetadata::units(). Referenced by pism::stressbalance::SSAFDFactory().