PISM, A Parallel Ice Sheet Model 2.2.2-d6b3a29ca committed by Constantine Khrulev on 2025-03-28
Loading...
Searching...
No Matches

◆ add_contribution() [1/2]

void pism::fem::Element::add_contribution ( const double K,
Mat  J 
) const

Add Jacobian contributions.

Add the contributions of an element-local Jacobian to the global Jacobian matrix.

The element-local Jacobian should be given as a row-major array of Nk*Nk values in the scalar case or (2Nk)*(2Nk) values in the vector valued case.

Note that MatSetValuesBlockedStencil ignores negative indexes, so values in K corresponding to locations marked using mark_row_invalid() and mark_col_invalid() are ignored. (Just as they should be.)

Definition at line 108 of file Element.cc.