[petsc-users] Nonlinear solve did not converge due to DIVERGED_LINEAR_SOLVE

李瑞 lirui319 at hnu.edu.cn
Mon Jul 17 09:16:25 CDT 2017


Dear professor or engineer:
   I recently attempt to solve the system of equation that consist of two sub equations,which are respetivly spatially and timely dependant.The structure of equations are likely : F_1(x,y)=0;F_2(w(t),x)=0.and y can be explicitly obtained in each time step.We need to pay attention to the first equation is a linear equation and the second a  nonlinear equation dependent on time.Getting W(t0) by sequentially obtain y and x.
   When I solve the equation F1 in first time step,the result was convergent,but in the second step was not.The answer reveals that:
>>>In the first time step:
      ....
  126 KSP Residual norm 6.505554772167e-03 
  127 KSP Residual norm 2.563866095029e-03 
  128 KSP Residual norm 5.365504095451e-04 
  129 KSP Residual norm 9.699344087450e-05 
  130 KSP Residual norm 2.793539286585e-05 
  131 KSP Residual norm 8.413249813721e-06 
  132 KSP Residual norm 2.488394296066e-06 
  133 KSP Residual norm 1.055822777756e-06 
  134 KSP Residual norm 7.820629842889e-07
Nonlinear solve converged due to CONVERGED_ITS

>>>In the second time step :
       ....
  294 KSP Residual norm 2.632117164436e+03 
  295 KSP Residual norm 2.627315190869e+03 
  296 KSP Residual norm 2.626944159914e+03 
  297 KSP Residual norm 2.625073489387e+03 
  298 KSP Residual norm 2.624035038684e+03 
  299 KSP Residual norm 2.623956156495e+03 
  300 KSP Residual norm 2.623696336837e+03
Nonlinear solve did not converge due to DIVERGED_LINEAR_SOLVE

The following code is typing my makefile.
-linear_ksp_monitor -linear_snes_converged_reason -linear_snes_atol 1.e-4 -linear_snes_rtol 1.e-60 -linear_ksp_type gmres \
-linear_ksp_max_it 300 -linear_ksp_atol 1.e-6 -linear_ksp_gmres_restart 90 -linear_ksp_rtol 1.e-60 \
-linear_ksp_pc_side right -linear_pc_type asm -linear_pc_asm_type restrict -linear_pc_asm_overlap 1 \
-linear_sub_ksp_type preonly -linear_sub_pc_type lu \

 Is there anything wrong,Pc_type or Krylov method?Can you give me some detailed suggestions or other ways to solve my obsession?Thank for your heartful help!


More information about the petsc-users mailing list