"Kong, Fande" <fande.k...@inl.gov> writes:

> Hi Developers,
> Any reason to force users provide a matrix which does not miss any diagonal
> entries when using a LU-type solver?

Automatically adding the entries for incomplete factorization is a
slight increase in complexity and means that "ILU(0)" actually has some
fill.  Note that a literal ILU(0) cannot work if there are missing
diagonal entries.

It's probably not a significant code complexity difference for LU.

> Sometime, it is impossible to have all diagonal entries in a matrix, that
> is, the matrix has to miss some diagonal entries. For example, there is a
> saddle-point matrix from the discretization of incomprehensible equations,
> and the lower part of the matrix is a zero block. The matrix usually looks
> like:
> | A   B^T  |
> | B    0     |

Just insert explicit zeros on the diagonal.

Attachment: signature.asc
Description: PGP signature

Reply via email to