PISM, A Parallel Ice Sheet Model  stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -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::Vector2d::u, and pism::Vector2d::v.

Referenced by compute_jacobian().