hecmw_solve_error.f90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
! Copyright (c) 2019 FrontISTR Commons
! This software is released under the MIT License, see LICENSE.txt
!-------------------------------------------------------------------------------
 
module m_hecmw_solve_error
  use hecmw_util
 
  integer(kind=kint), parameter :: &
    hecmw_solver_error_incons_pc = 1001, &
    hecmw_solver_error_zero_diag = 2001, &
    hecmw_solver_error_zero_rhs = 2002, &
    hecmw_solver_error_noconv_maxit = 3001, &
    hecmw_solver_error_diverge_mat = 3002, &
    hecmw_solver_error_diverge_pc = 3003, &
    hecmw_solver_error_diverge_nan = 3004
 
contains
 
  subroutine hecmw_solve_error (hecMESH, IFLAG)
    use hecmw_util
    implicit none
    type (hecmwST_local_mesh) :: hecMESH
    integer(kind=kint) :: IFLAG
 
    if ( hecmesh%zero.eq.1 .and. &
        (iflag.eq.hecmw_solver_error_incons_pc .or. iflag.eq.hecmw_solver_error_zero_diag) ) then
      write (*,'(/a )')                                               &
        &           '###############################################'
      write (*,'( a )')                                               &
        &           '######## ERROR MESSAGE : LINEAR SOLVER ########'
      write (*,'( a/)')                                               &
        &           '###############################################'
    endif
 
    if (iflag.eq.hecmw_solver_error_incons_pc) then
      if (hecmesh%zero.eq.1) then
        write (*,'(/a )')'  #### HEC-MW-SOLVER-E-1001'
        write (*,'( a/)')'    inconsistent solver/preconditioning'
      endif
      !        call MPI_ABORT (hecMESH%MPI_COMM, ierr)
      call hecmw_abort( hecmw_comm_get_comm())
    endif
 
    if (iflag.eq.hecmw_solver_error_zero_diag) then
      if (hecmesh%zero.eq.1) then
        write (*,'(/a )')'  #### HEC-MW-SOLVER-E-2001: '
        write (*,'( a/)')'    ZERO component in diagonal block'
      endif
      !        call MPI_ABORT (hecMESH%MPI_COMM, ierr)
      call hecmw_abort( hecmw_comm_get_comm())
    endif
 
    if (iflag.eq.hecmw_solver_error_zero_rhs) then
      if (hecmesh%zero.eq.1) then
        write (*,'(/a )')'  #### HEC-MW-SOLVER-W-2002: '
        write (*,'( a/)')'    ZERO RHS norm'
      endif
    endif
 
    if (iflag.eq.hecmw_solver_error_noconv_maxit) then
      if (hecmesh%zero.eq.1) then
        write (*,'(/a )')'  #### HEC-MW-SOLVER-W-3001: '
        write (*,'( a/)')'    not converged within certain iterations'
      endif
    endif
 
    if (iflag.eq.hecmw_solver_error_diverge_mat) then
      if (hecmesh%zero.eq.1) then
        write (*,'(/a )')'  #### HEC-MW-SOLVER-W-3002: '
        write (*,'( a/)')'    diverged due to indefinite or negative definite matrix'
      endif
    endif
 
    if (iflag.eq.hecmw_solver_error_diverge_pc) then
      if (hecmesh%zero.eq.1) then
        write (*,'(/a )')'  #### HEC-MW-SOLVER-W-3003: '
        write (*,'( a/)')'    diverged due to indefinite preconditioner'
      endif
    endif
 
    if (iflag.eq.hecmw_solver_error_diverge_nan) then
      if (hecmesh%zero.eq.1) then
        write (*,'(/a )')'  #### HEC-MW-SOLVER-W-3004: '
        write (*,'( a/)')'    diverged due to NaN'
      endif
    endif
 
    !stop " PROGRAM STOP:"
 
  end subroutine hecmw_solve_error
 
end module m_hecmw_solve_error