Dear Prof. Wolfgang, 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)?
Regards Phani On Mon, Apr 2, 2018 at 2:28 PM, Wolfgang Bangerth <[email protected]> wrote: > On 04/02/2018 11:48 AM, Phani Motamarri wrote: > >> We have an ongoing research project with Caltech where we are trying to >> explore ways of approximating the subspace spanned by k lowest eigen values >> of a given positive definite matrix(FEM discretized matrix) in a sparse >> format using the incomplete Cholesky factor with reduced computational >> complexity. So I need access to the incomplete Cholesky factor. >> > > Interesting. > > The "cheap and dirty" way to get what you want is to modify the > declaration of the class so that the derivation of SparseILU from > SparseMatrix is 'public', not 'private' or 'protected'. This way, you can > access all of the usual member functions of the sparse matrix. We don't > want to do that for the actual sources we distribute for the reasons stated > previously, but there is nothing inherently wrong with doing this yourself > in your local copy. > > > 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/fo > rum/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. > -- 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.
