PISM, A Parallel Ice Sheet Model  stable v2.1.1 committed by Constantine Khrulev on 2024-12-04 13:36:58 -0900

◆ 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::Vector2d::u, and pism::Vector2d::v.

Referenced by compute_jacobian().