PISM, A Parallel Ice Sheet Model  stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800

◆ residual_dirichlet()

void pism::stressbalance::Blatter::residual_dirichlet ( const DMDALocalInfo &  info,
Parameters **  P,
const Vector2d ***  x,
Vector2d ***  R 
)
protected

Set the residual at Dirichlet locations

Compute the residual at Dirichlet locations and reset the residual to zero elsewhere.

Setting it to zero is necessary because we call DMDASNESSetFunctionLocal() with INSERT_VALUES.

Definition at line 262 of file residual.cc.

References pism::stressbalance::Blatter::Parameters::bed, dirichlet_node(), pism::grid_z(), pism::k, pism::Component::m_grid, pism::NODE_EXTERIOR, pism::stressbalance::Blatter::Parameters::thickness, u_bc(), and pism::stressbalance::u_exterior.

Referenced by compute_residual().