Added rnorm output in NewtonKrylov bypass.

nekStab · Dec 11, 2024 · d4a742b · d4a742b
1 parent 3e2039d
commit d4a742b
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 0 deletions.
diff --git a/src/NewtonKrylov.f90 b/src/NewtonKrylov.f90
@@ -189,6 +189,8 @@ subroutine newton_rsp(sys, X, solver, info, rtol, atol, options, linear_solver_o
 
       ! Check for lucky convergence.
       if (rnorm < target_tol) then
+         write(msg,'(2(A,E9.2))') 'nrorm= ', rnorm, ', target= ', target_tol
+         call logger%log_warning(msg, module=this_module, procedure='newton_rsp')
          write(msg,'(A)') 'Initial guess is a fixed point to tolerance!'
          call logger%log_warning(msg, module=this_module, procedure='newton_rsp')
          newton_meta%converged = .true.
@@ -352,6 +354,8 @@ subroutine newton_rdp(sys, X, solver, info, rtol, atol, options, linear_solver_o
 
       ! Check for lucky convergence.
       if (rnorm < target_tol) then
+         write(msg,'(2(A,E9.2))') 'nrorm= ', rnorm, ', target= ', target_tol
+         call logger%log_warning(msg, module=this_module, procedure='newton_rdp')
          write(msg,'(A)') 'Initial guess is a fixed point to tolerance!'
          call logger%log_warning(msg, module=this_module, procedure='newton_rdp')
          newton_meta%converged = .true.
@@ -515,6 +519,8 @@ subroutine newton_csp(sys, X, solver, info, rtol, atol, options, linear_solver_o
 
       ! Check for lucky convergence.
       if (rnorm < target_tol) then
+         write(msg,'(2(A,E9.2))') 'nrorm= ', rnorm, ', target= ', target_tol
+         call logger%log_warning(msg, module=this_module, procedure='newton_csp')
          write(msg,'(A)') 'Initial guess is a fixed point to tolerance!'
          call logger%log_warning(msg, module=this_module, procedure='newton_csp')
          newton_meta%converged = .true.
@@ -678,6 +684,8 @@ subroutine newton_cdp(sys, X, solver, info, rtol, atol, options, linear_solver_o
 
       ! Check for lucky convergence.
       if (rnorm < target_tol) then
+         write(msg,'(2(A,E9.2))') 'nrorm= ', rnorm, ', target= ', target_tol
+         call logger%log_warning(msg, module=this_module, procedure='newton_cdp')
          write(msg,'(A)') 'Initial guess is a fixed point to tolerance!'
          call logger%log_warning(msg, module=this_module, procedure='newton_cdp')
          newton_meta%converged = .true.

diff --git a/src/NewtonKrylov.fypp b/src/NewtonKrylov.fypp
@@ -178,6 +178,8 @@ contains
 
       ! Check for lucky convergence.
       if (rnorm < target_tol) then
+         write(msg,'(2(A,E9.2))') 'nrorm= ', rnorm, ', target= ', target_tol
+         call logger%log_warning(msg, module=this_module, procedure='newton_${type[0]}$${kind}$')
          write(msg,'(A)') 'Initial guess is a fixed point to tolerance!'
          call logger%log_warning(msg, module=this_module, procedure='newton_${type[0]}$${kind}$')
          newton_meta%converged = .true.