PISM, A Parallel Ice Sheet Model  stable v2.0.5 committed by Constantine Khrulev on 2022-10-14 09:56:26 -0800

◆ jacobian_dirichlet()

void pism::stressbalance::Blatter::jacobian_dirichlet ( const DMDALocalInfo &  info,
Parameters **  P,
Mat  J 
)
protected

Set the Jacobian to identity at Dirichlet nodes.

Definition at line 168 of file jacobian.cc.

References dirichlet_node(), J, pism::k, pism::NODE_EXTERIOR, PISM_CHK, pism::fem::row(), pism::Vector2::u, and pism::Vector2::v.

Referenced by compute_jacobian().