On 04/02/2018 02:59 PM, Phani Motamarri wrote:
I am trying to modify the function "SparseILU<number>::initialize"
inside sparse_ilu.templates.h to call the mkl function "dcsrilu0" which
requires the use of some of the data like ia, ja accessed inside
"initialize" function. The tricky part is that the sparsity pattern
stores the diagonal element at the first position of every row and the
remaining non-zero entries are stored next . Do you have any private
members of the SparsityPattern/SparseMatrix function which gives an
array of non-zero entries directly in CSR format (A,IA,JA)?
No -- you'll have to do the sorting yourself I'm afraid.
It's possible that you can find code that does this sort of thing in
sparse_direct.cc where we already interface with sparse direct solvers
that all want their own format. I don't recall whether they have the
sortedness restriction, but either way, that's the place to look for
related code.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
