<div dir="ltr"><div dir="ltr"> Hi Aditya,<div><br></div><div>I wouldn't trust FEniCS' wrappers around PETSc. It might also depend on the version of DOLFIN you're working with.</div><div><br></div><div>Shameless plug: Try using pfibs </div><div><br></div><div><a href="https://github.com/NREL/pfibs" target="_blank">https://github.com/NREL/pfibs</a></div><div><br></div><div> We wrote a slightly better PETSc interface to FEniCS. It does require FEniCS 2018.1.0 or higher. This should at least give better error logs for your problem.</div><div><br></div><div>Also I know a lot of PETSc solvers rely on some sort of DM, which FEniCS (at the moment) does not automatically provide. We have a pFibs branch which attaches a DMShell to your functionspace. </div><div><br></div><div><a href="https://github.com/NREL/pfibs/tree/dm-nested-solver-feature" target="_blank">https://github.com/NREL/pfibs/tree/dm-nested-solver-feature</a></div><div><br></div><div>It's still in its early stages, but it should at least be able to handle the PETScOptions you've set above.</div><div><br></div><div>Hopefully this helps,</div><div>Justin</div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Feb 8, 2019 at 10:41 AM aditya kumar via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><p style="margin-top:0px;font-family:Helvetica,Arial,sans-serif;font-size:15px">Hello,</p><p style="font-family:Helvetica,Arial,sans-serif;font-size:15px">I am using PETSc with FEniCS project libraries to solve a nonlinear problem. I am using PETSc Krylov solver with the following configuration</p><p style="font-family:Helvetica,Arial,sans-serif;font-size:15px"><span style="font-family:Consolas,Menlo,Monaco,"Lucida Console","Liberation Mono","DejaVu Sans Mono","Bitstream Vera Sans Mono","Courier New",monospace;white-space:pre-wrap;background-color:rgb(249,249,249)">pc = PETScPreconditioner("petsc_amg") 

PETScOptions.set("mg_levels_ksp_type", "chebyshev")

PETScOptions.set("mg_levels_pc_type", "jacobi")

PETScOptions.set("mg_levels_esteig_ksp_type", "cg")

PETScOptions.set("mg_levels_ksp_chebyshev_esteig_steps", 50)

solver_u = PETScKrylovSolver("cg", pc)

<p style="font-size:15px"><span style="font-family:Helvetica,Arial,sans-serif">I am encountering this error sometimes during the solution of the linear system:<span style="font-family:Consolas,Menlo,Monaco,"Lucida Console","Liberation Mono","DejaVu Sans Mono","Bitstream Vera Sans Mono","Courier New",monospace;font-size:15px;white-space:pre-wrap;background-color:rgb(249,249,249)">Error:   Unable to successfully call PETSc function 'KSPSolve'.

*** Reason:  PETSc error code is: 77 (Petsc has generated inconsistent data).

*** Where:   This error was encountered inside /tmp/dolfin/dolfin/la/PETScKrylovSolver.cpp. </span><p style="font-size:15px"><span style="background-color:rgb(249,249,249);font-family:Consolas,Menlo,Monaco,"Lucida Console","Liberation Mono","DejaVu Sans Mono","Bitstream Vera Sans Mono","Courier New",monospace;white-space:pre-wrap">*** Process: 130</span></p><p style="font-size:15px">Can anybody help me in understanding the possible reasons for this error code? Or refer me to any available documentation? It will be really helpful. </p><p style="font-size:15px">Also, I will like to know if there is a way to <span style="font-family:Helvetica,Arial,sans-serif">prevent the code from stopping due to such PETSc errors, so that I can tweak some parameters to hopefully make it progress.</span>  </p><div><br></div><div>Thanks and Regards,</div><div>Aditya</div><div dir="ltr" class="gmail-m_-4213504756596182550gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><br></div></div></div></div></div></div>

</blockquote></div>