Blatter compute_averaged_velocity compute_jacobian compute_node_type compute_residual copy_solution define_model_state_impl 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 u_bc update velocity_u_sigma velocity_v_sigma write_model_state_impl ~Blatter

compute_jacobian() void pism::stressbalance::Blatter::compute_jacobian ( DMDALocalInfo *  petsc_info, const Vector2d ***  X, Mat  A, Mat  J  ) protected Compute the Jacobian matrix. Definition at line 199 of file jacobian.cc. References pism::fem::Element3::add_contribution(), dirichlet_node(), exterior_element(), pism::fem::q13d::FACE_BOTTOM, pism::GHOSTED, pism::grid_transpose(), pism::grid_z(), grounding_line(), I, J, jacobian_basal(), jacobian_dirichlet(), jacobian_f(), pism::hydrology::K(), pism::k, pism::fem::Element3::local_to_global(), m_face100, m_face4, pism::Component::m_grid, m_Nq, m_parameters, pism::fem::Q1Element3::mark_col_invalid(), pism::fem::Q1Element3::mark_row_invalid(), n, pism::fem::Element::n_chi(), pism::fem::q13d::n_chi, pism::fem::Element::n_pts(), nodal_parameter_values(), pism::fem::Element3::nodal_values(), PISM_CHK, pism::fem::Q1Element3Face::reset(), pism::fem::Q1Element3::reset(), u_bc(), pism::stressbalance::ShallowStressBalance::velocity(), pism::fem::Q1Element3::x(), pism::fem::Q1Element3::y(), and pism::fem::Q1Element3::z(). Referenced by jacobian_callback().