Well ... In the same run (other timesteps), I also see DIVERGED_INDEFINITE_MAT. In many timesteps it converges. It happens sometimes... I am using LU preconditioners.
I have a 2D circumference inside the domain (petroleum wellbore). I need to set a Neumann BC (total oil rate), but the pressure must be constant in the circumference. Similar to a distributed total force in a rigid plate - for mechanical solvers. My approach: 1) select a DOF for the circumference. 2) tie all other DOFs in the circumference to this one with "add_constrain_row". Can you see another way to do that? I am in try-and-try approach already... Renato On Wed, Mar 13, 2019 at 10:34 AM Stogner, Roy H <royst...@ices.utexas.edu> wrote: > > On Tue, 12 Mar 2019, Renato Poli wrote: > > > I see the message below, but I cannot identify the reason. > > I know it is related to constrain_rows" I use to tie together the > pressure > > in a boundary. > > > > 0 KSP Residual norm 8.180571722215e+07 > > 1 KSP Residual norm 4.475923747576e+03 > > 2 KSP Residual norm 4.475912845013e+03 > > # Debug[1:0]: Linear solver convergence/divergence reason: > > DIVERGED_INDEFINITE_PC > > > > What that means? > > Any suggestion on possible causes? > > IIRC that usually means you don't have an invertible preconditioner, > which unless you were setting up a separate preconditioner means you > don't have an invertible Jacobian. > --- > Roy > _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
