<div dir="ltr"><div class="gmail_quote"><div dir="ltr">On Tue, Oct 30, 2018 at 5:21 PM Appel, Thibaut 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear users,<br>
<br>
Following a suggestion from Matthew Knepley I’ve been trying to apply fieldsplit/gamg for my set of PDEs but I’m still encountering issues despite various tests. pc_gamg simply won’t start.<br>
Note that direct solvers always yield the correct, physical result.<br>
Removing the fieldsplit to focus on the gamg bit </blockquote><div><br></div><div>This is a mistake I think because you want a simpler operator to try GAMG on.</div><div><br></div><div>Can you go back to splitting the system, and try using GMRES as the smoother. Its important to see</div><div>whether the smoother makes no progress, or the coarse correction stinks.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">and trying to solve the linear system on a modest size problem still gives, with<br>
<br>
'-ksp_monitor -ksp_rtol 1.0e-10 -ksp_gmres_restart 300 -ksp_type gmres -pc_type gamg'<br>
<br>
[3]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<br>
[3]PETSC ERROR: Petsc has generated inconsistent data<br>
[3]PETSC ERROR: Have un-symmetric graph (apparently). Use '-(null)pc_gamg_sym_graph true' to symetrize the graph or '-(null)pc_gamg_threshold -1' if the matrix is structurally symmetric.<br>
<br>
And since then, after adding '-pc_gamg_sym_graph true' I have been getting<br>
[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<br>
[0]PETSC ERROR: Petsc has generated inconsistent data<br>
[0]PETSC ERROR: Eigen estimator failed: DIVERGED_NANORINF at iteration<br>
<br>
-ksp_chebyshev_esteig_noisy 0/1 does not change anything<br>
<br>
Knowing that Chebyshev eigen estimator needs a positive spectrum I tried ‘-mg_levels_ksp_type gmres’ but iterations would just go on endlessly.<br>
<br>
It seems that I have indeed eigenvalues of rather high magnitude in the spectrum of my operator without being able to determine the reason.<br>
The eigenvectors look like small artifacts at the wall-inflow or wall-outflow corners with zero anywhere else but I do not know how to interpret this.<br>
Equations are time-harmonic linearized Navier-Stokes to which a forcing is applied, there’s no time-marching.<br>
<br>
Matrix is formed with a MPIAIJ type. The formulation is incompressible, in complex arithmetic and the 2D physical domain is mapped to a logically rectangular, regular collocated grid with a high-order finite difference method.<br>
I determine the ownership of the rows/degrees of freedom of the matrix with PetscSplitOwnership and I’m not using DMDA.<br>
<br>
The Fortran application code is memory-leak free and has undergone a strict verification/validation procedure for different variations of the PDEs.<br>
<br>
If there’s any problem with the matrix what could help for the diagnostic? At this point I’m running out of ideas so I would really appreciate additional suggestions and discussions.<br>
<br>
Thanks for your continued support,<br>
<br>
<br>
Thibaut</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>