apply_jacobian_design apply_jacobian_design apply_jacobian_design apply_jacobian_design_transpose apply_jacobian_design_transpose apply_jacobian_design_transpose apply_linearization apply_linearization_transpose assemble_jacobian_state assemble_residual assemble_residual DesignVec DesignVecGhosted get_da init IP_SSATaucForwardProblem linearize_at m_du_global m_du_local m_dzeta_local m_element m_element_index m_fixed_tauc_locations m_J_state m_ksp m_reason m_rebuild_J_state m_tauc_copy m_tauc_param m_velocity_shared m_zeta set_design set_tauc_fixed_locations solution StateVec StateVec1 tauc_param ~IP_SSATaucForwardProblem

apply_jacobian_design() [3/3] void pism::inverse::IP_SSATaucForwardProblem::apply_jacobian_design ( array::Vector &  u, array::Scalar &  dzeta, Vector2d **  du_a  ) virtual Applies the design Jacobian matrix to a perturbation of the design variable. The matrix depends on the current value of the design variable \(\zeta\) and the current value of the state variable \(u\). The specification of \(\zeta\) is done earlier with set_design or linearize_at. The value of \(u\) is specified explicitly as an argument to this method. Parameters [in] u Current state variable value. [in] dzeta Perturbation of the design variable. Prefers vectors with ghosts; will copy to a ghosted vector if needed. [out] du_a Computed corresponding perturbation of the state variable. The array du_a should be extracted first from a Vec or an array::Array. Typically this method is called via one of its overloads. Definition at line 220 of file IP_SSATaucForwardProblem.cc. References pism::array::Array2D< T >::add(), pism::fem::Element2::add_contribution(), pism::array::Array2D< T >::array(), pism::ParallelSection::check(), pism::fem::Element::chi(), pism::fem::DirichletData::constrain(), pism::array::Array2D< T >::copy_from(), pism::IceBasalResistancePlasticLaw::drag(), pism::fem::DirichletData_Vector::enforce(), pism::fem::DirichletData_Scalar::enforce_homogeneous(), pism::fem::Element2::evaluate(), pism::ParallelSection::failed(), pism::fem::DirichletData_Vector::fix_residual_homogeneous(), pism::mask::grounded_ice(), pism::k, pism::stressbalance::ShallowStressBalance::m_basal_sliding_law, pism::stressbalance::SSAFEM::m_bc_mask, pism::stressbalance::SSAFEM::m_bc_values, pism::stressbalance::SSAFEM::m_coefficients, pism::stressbalance::SSAFEM::m_dirichletScale, m_dzeta_local, m_element, m_element_index, m_fixed_tauc_locations, pism::Component::m_grid, m_tauc_param, m_zeta, pism::fem::MAX_QUADRATURE_SIZE, pism::fem::q1::n_chi, pism::fem::Element::n_pts(), pism::fem::Element2::nodal_values(), pism::stressbalance::SSAFEM::quad_point_values(), pism::fem::Element2::reset(), pism::array::Array::stencil_width(), pism::inverse::IPDesignVariableParameterization::toDesignVariable(), pism::Vector2d::u, pism::Vector2d::v, pism::fem::Germ::val, pism::fem::Element::weight(), pism::fem::ElementIterator::xm, pism::fem::ElementIterator::xs, pism::fem::ElementIterator::ym, and pism::fem::ElementIterator::ys.