Phani, On Monday, April 2, 2018 at 4:59:02 PM UTC-4, 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)? > This is already done for UMFPACK. Take a look here https://github.com/dealii/dealii/blob/master/source/lac/sparse_direct.cc#L106
Best, Bruno -- 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.
