<div dir="ltr">Since I can't send binary files, I hardcoded the values on main.py. You only need main.py file to run the code.<div><br></div><div>Thanks,</div></div><br><div class="gmail_quote"><div dir="ltr">On Mon, Nov 12, 2018 at 5:17 PM HeeHo Park <<a href="mailto:heeho.park@gmail.com">heeho.park@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi PETSc, <div><br></div><div>These Python looks complicated but it really isn't. please run main.py</div><div><br></div><div>I am writing up a script that would read petsc binary files extracted from one of PFLOTRAN newton iteration step and solve with petsc4py to investigate linear solvers.</div><div><br></div><div>But I am having trouble solving a matrix that converges PFLOTRAN, let alone running a trimmed size (15x15) matrix of A. You'll see in the script.  </div><div><br></div><div>I don't think there is any problem with loading sparse.csr_matrix onto PETSc and same with array as I tested with very simple matrix which it solves.</div><div><br></div><div>Can you figure out why I'm getting inf solutions?</div><div><br></div><div>If I use the sparse matrix package to solve it, it solves it very quickly.</div><div><br></div><div>x = sla.spsolve(A,b)   </div><div><br></div><div>Thanks,<br><div><div><br></div>-- <br><div dir="ltr" class="m_-6457269117685951849gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">HeeHo Daniel Park</div></div></div></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">HeeHo Daniel Park</div></div>