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

◆ add_contribution()

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

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 232 of file Element.cc.

References J, pism::hydrology::K(), m_block_size, m_col, m_row, and PISM_CHK.