SimPEG.flow.richards.SimulationNDCellCentered.diagsJacobian#

SimulationNDCellCentered.diagsJacobian(m, hn, hn1, dt, bc)[source]#

Diagonals and rhs of the jacobian system

The matrix that we are computing has the form:

.. code::

.- -. .- -. .- -. | Adiag | | h1 | | b1 | | Asub Adiag | | h2 | | b2 | | Asub Adiag | | h3 | = | b3 | | … … | | .. | | .. | | Asub Adiag | | hn | | bn | ‘- -’ ‘- -’ ‘- -’