Blatter compute_averaged_velocity compute_jacobian compute_node_type compute_residual copy_solution dirichlet_node exterior_element function_callback get_basal_velocity grounding_line init_2d_parameters init_ice_hardness init_impl jacobian_basal jacobian_callback jacobian_dirichlet jacobian_f m_da m_E_viscosity m_eta_transform m_face100 m_face4 m_glen_exponent m_ksp_use_ew m_n_work m_Nq m_parameters m_rho_ice_g m_rho_ocean_g m_scaling m_snes m_u_sigma m_v_sigma m_viscosity_eps m_work m_work2 m_x m_x_old marine_boundary nodal_parameter_values parameter_continuation partially_submerged_face report_mesh_info residual_basal residual_dirichlet residual_f residual_lateral residual_source_term residual_surface set_initial_guess setup solve state_impl u_bc update velocity_u_sigma velocity_v_sigma write_state_impl ~Blatter

residual_dirichlet() void pism::stressbalance::Blatter::residual_dirichlet ( const DMDALocalInfo &  info, Parameters **  P, const Vector2d ***  x, Vector2d ***  R  ) protected Set the residual at Dirichlet locations Compute the residual at Dirichlet locations and reset the residual to zero elsewhere. Setting it to zero is necessary because we call DMDASNESSetFunctionLocal() with INSERT_VALUES. Definition at line 263 of file residual.cc. References pism::stressbalance::Blatter::Parameters::bed, dirichlet_node(), pism::grid_z(), pism::k, pism::Component::m_grid, pism::NODE_EXTERIOR, pism::stressbalance::Blatter::Parameters::node_type, pism::stressbalance::Blatter::Parameters::thickness, u_bc(), and pism::stressbalance::u_exterior. Referenced by compute_residual().