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

is_marginal() bool pism::stressbalance::SSAFD::is_marginal ( int  i, int  j, bool  ssa_dirichlet_bc  ) protectedvirtual Checks if a cell is near or at the ice front. You need to create array::AccessScope object and add mask to it. Note that a cell is a CFBC location of one of four direct neighbors is ice-free. If one of the diagonal neighbors is ice-free we don't use the CFBC, but we do need to compute weights used in the SSA discretization (see assemble_matrix()) to avoid differencing across interfaces between icy and ice-free cells. This method ensures that checks in assemble_rhs() and assemble_matrix() are consistent. Definition at line 1698 of file SSAFD.cc. References pism::array::CellType1::box_int(), pism::mask::ice_free(), pism::mask::ice_free_ocean(), pism::mask::icy(), and pism::stressbalance::SSA::m_mask. Referenced by assemble_matrix(), and assemble_rhs().