fluid_aux.f90

Go to the documentation of this file.

 
module fluid_aux
  use logger, only : neko_log, log_size
  use num_types, only : rp
  use krylov, only : ksp_monitor_t
  use, intrinsic :: ieee_arithmetic, only: ieee_is_nan
  use utils, only : neko_error, neko_warning
  implicit none
  private
 
  public :: fluid_step_info
 
contains
 
  subroutine fluid_step_info(step, t, dt, ksp_results, strict_convergence)
    type(ksp_monitor_t), intent(in) :: ksp_results(4)
    integer, intent(in) :: step
    real(kind=rp), intent(in) :: t, dt
    logical, intent(in) :: strict_convergence
    character(len=LOG_SIZE) :: log_buf
    integer :: i
 
    call neko_log%message('Pressure')
 
    write(log_buf, '(A,A,A)') 'Iterations:   ',&
         'Start residual:     ', 'Final residual:'
    call neko_log%message(log_buf)
    write(log_buf, '(I11,3x, E15.7,5x, E15.7)') ksp_results(1)%iter, &
         ksp_results(1)%res_start, ksp_results(1)%res_final
    call neko_log%message(log_buf)
 
    call neko_log%message('X-Velocity')
    write(log_buf, '(A,A,A)') 'Iterations:   ',&
         'Start residual:     ', 'Final residual:'
    call neko_log%message(log_buf)
    write(log_buf, '(I11,3x, E15.7,5x, E15.7)') ksp_results(2)%iter, &
         ksp_results(2)%res_start, ksp_results(2)%res_final
    call neko_log%message(log_buf)
 
    call neko_log%message('Y-Velocity')
    write(log_buf, '(A,A,A)') 'Iterations:   ',&
         'Start residual:     ', 'Final residual:'
    call neko_log%message(log_buf)
    write(log_buf, '(I11,3x, E15.7,5x, E15.7)') ksp_results(3)%iter, &
         ksp_results(3)%res_start, ksp_results(3)%res_final
    call neko_log%message(log_buf)
 
    call neko_log%message('Z-Velocity')
    write(log_buf, '(A,A,A)') 'Iterations:   ', &
         'Start residual:     ', 'Final residual:'
    call neko_log%message(log_buf)
    write(log_buf, '(I11,3x, E15.7,5x, E15.7)') ksp_results(4)%iter, &
         ksp_results(4)%res_start, ksp_results(4)%res_final
    call neko_log%message(log_buf)
 
    ! Check for convergence
    do i = 1, 4
       if (ieee_is_nan(ksp_results(i)%res_final)) then
          call neko_error("Fluid solver diverged")
       end if
 
       if (.not. ksp_results(i)%converged) then
          log_buf = 'Fluid solver did not converge for '
          select case(i)
            case(1)
             log_buf = trim(log_buf) // 'pressure'
            case(2)
             log_buf = trim(log_buf) // 'x-velocity'
            case(3)
             log_buf = trim(log_buf) // 'y-velocity'
            case(4)
             log_buf = trim(log_buf) // 'z-velocity'
          end select
 
          if (strict_convergence) then
             call neko_error(log_buf)
          else
             call neko_warning(log_buf)
          end if
       end if
    end do
 
  end subroutine fluid_step_info
 
end module fluid_aux