Does the option to keep the diagonal for non-square matrices in assembly (i.e., no not let the linear algebra backend optimise it away) make much sense? It causes all sorts of problems with local (process-wise) indexing of dofs because the ownership ranges of the row and column maps differ. It would be handy to remove support for preserving the diagonal for non-square matrices.

Garth

_______________________________________________
fenics mailing list
[email protected]
http://fenicsproject.org/mailman/listinfo/fenics

Reply via email to