fluid_aux.f90

Go to the documentation of this file.

 
module fluid_aux
  use, intrinsic :: ieee_arithmetic, only: ieee_is_nan
  use krylov, only : ksp_monitor_t
  use logger, only : log_size
  use utils, only : neko_error, neko_warning
  use time_state, only : time_state_t
  implicit none
  private
 
  public :: fluid_step_info
 
contains
 
  subroutine fluid_step_info(time, ksp_results, full_stress_formulation, &
       strict_convergence)
    type(ksp_monitor_t), dimension(:), intent(in) :: ksp_results
    type(time_state_t), intent(in) :: time
    logical, intent(in) :: full_stress_formulation
    logical, intent(in), optional :: strict_convergence
    character(len=LOG_SIZE) :: log_buf
    integer :: i, n
 
    n = size(ksp_results)
    if (full_stress_formulation) n = 2
 
    ! Do the printing
    call ksp_results(1)%print_header()
    do i = 1, n
       call ksp_results(i)%print_result(time%tstep)
    end do
 
    ! Check for convergence
    do i = 1, n
       if (ieee_is_nan(ksp_results(i)%res_final)) then
          call neko_error("Fluid solver diverged for " // &
               trim(ksp_results(i)%name))
       end if
 
       if (present(strict_convergence)) then
 
          if (.not. ksp_results(i)%converged) then
             log_buf = 'Fluid solver did not converge for ' &
                  // trim(ksp_results(i)%name)
 
             if (strict_convergence) then
                call neko_error(log_buf)
             else
                call neko_warning(log_buf)
             end if
          end if
       end if
    end do
 
  end subroutine fluid_step_info
 
end module fluid_aux