<div dir="ltr">Thanks for the input. I was doing a few tests. I will put quadruple precision choice on hold for now. Since the condition number is big at the beginning, scale of 10^2, I will try to smooth the input of the matrix by introducing heaviside function.  I do have another question. the incompressible simulation I'm simulating is using all reflection(free slip neumann) boundary condition, what will be a good preconditioner? Shall I start a new thread? so, when other users search, they could use it as a reference<br></div>