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

◆ 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().