cache_inputs cache_residual_cfbc compute_local_function compute_local_jacobian driving_stress explicit_driving_stress function_callback init_impl jacobian_callback m_alpha m_bc_mask m_bc_values m_beta_ice_free_bedrock m_boundary_integral m_callback_data m_coefficients m_dirichletScale m_driving_stress_x m_driving_stress_y m_element_index m_epsilon_ssa m_gc m_node_type m_q1_element m_rho_g m_snes monitor_function monitor_jacobian PointwiseNuHAndBeta quad_point_values solve solve_nocache solve_with_reason SSAFEM ~SSAFEM

solve() void pism::stressbalance::SSAFEM::solve ( const Inputs &  inputs ) protectedvirtual Solve the SSA system of equations. The FEM solver computes values of the various coefficients (most notably: the vertically-averaged ice hardness) at each of the element nodes only once. When running in an ice-model context, at each time step, SSA::update() is called, which calls SSAFEM::solve(). Since coefficients have generally changed between time steps, we need to recompute coefficients. On the other hand, in the context of inversion, coefficients will not change between iteration and there is no need to recompute their values. So there are two different solve methods, SSAFEM::solve() and SSAFEM::solve_nocache(). The only difference is that SSAFEM::solve() recomputes the cached values of the coefficients before calling SSAFEM::solve_nocache(). Implements pism::stressbalance::SSA. Definition at line 171 of file SSAFEM.cc. References pism::RuntimeError::formatted(), pism::Component::m_log, pism::stressbalance::SSA::m_stdout_ssa, PISM_ERROR_LOCATION, and solve_with_reason().