[petsc-users] Solving NON-Diagonally dominant sparse system

[Kaushik Kulkarni]

Hello,
I am trying to solve a 2500x2500 sparse matrix. To get an idea about the
matrix structure I have added a file matrix.log which contains the output
of MatView() and also the output of Matview_draw in the image file.

>From the matrix structure it can be seen that Jacobi iteration won't work
and some of the diagonal entries being very low(of the order of 1E-16) LU
factorization would also fail.

C
​an someone please suggest what all could I try next, in order to make the
solution converge?

Thanks,
Kaushik
​

-- 
Kaushik Kulkarni
Fourth Year Undergraduate
Department of Mechanical Engineering
Indian Institute of Technology Bombay
Mumbai, India
https://kaushikcfd.github.io/About/
+91-9967687150
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20170410/0702de3a/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: matrix.log
Type: text/x-log
Size: 2091693 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20170410/0702de3a/attachment-0001.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: matrix_pattern.png
Type: image/png
Size: 25197 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20170410/0702de3a/attachment-0001.png>