Quantcast
Channel: Math.NET Numerics
Viewing all articles
Browse latest Browse all 971

New Post: ILU0 using and sparse matrix operations

$
0
0
Hello. I have some problems with solving system of linear equations. Really I have matrix about 20000*20000. I create Sparse Matrix A, Initialized it and created right hand side - vector B.

There is my code:
        var solver = new BiCgStab();

        var preconditioner = new ILU0Preconditioner();
        preconditioner.Initialize(A_sparse);
        var X = A_sparse.SolveIterative(BB, solver,monitor, preconditioner).ToArray();
But it work more than 30 minutes, I have never waited for it finish. When I use MILU0Preconditioner, solver do it very fast (5 -10 seconds) but the result is not that I need. I tryied in Matlab same algorithm and same matrix (I used bicg and LU preconditioner). Without preconditioner bicg in Matlab get similar result. But with LU preconditioner Matlab get very good results, unlike Math Net. I can't understand why, maybe somebody have experience with ILU0 preconditioner in Math Met Numerics

Viewing all articles
Browse latest Browse all 971

Latest Images

Trending Articles



Latest Images

<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>