hecmw_precond_DIAG.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 hecmw_precond_diag
  use hecmw_util
  use hecmw_precond_diag_11
  use hecmw_precond_diag_22
  use hecmw_precond_diag_33
  use hecmw_precond_diag_44
  use hecmw_precond_diag_66
  use hecmw_precond_diag_nn
  implicit none
 
  private
  public :: hecmw_precond_diag_setup
  public :: hecmw_precond_diag_clear
  public :: hecmw_precond_diag_apply
 
contains
 
  subroutine hecmw_precond_diag_setup(hecMAT)
    implicit none
    type (hecmwst_matrix), intent(inout) :: hecmat
 
    select case(hecmat%NDOF)
      case(3)
        call hecmw_precond_diag_33_setup(hecmat)
      case(4)
        call hecmw_precond_diag_44_setup(hecmat)
      case(6)
        call hecmw_precond_diag_66_setup(hecmat)
      case(1)
        call hecmw_precond_diag_11_setup(hecmat)
      case(2)
        call hecmw_precond_diag_22_setup(hecmat)
      case default
        call hecmw_precond_diag_nn_setup(hecmat)
    end select
  end subroutine hecmw_precond_diag_setup
 
  subroutine hecmw_precond_diag_clear(NDOF)
    implicit none
    integer(kind=kint), intent(in) :: ndof
 
    select case(ndof)
      case(3)
        call hecmw_precond_diag_33_clear()
      case(4)
        call hecmw_precond_diag_44_clear()
      case(6)
        call hecmw_precond_diag_66_clear()
      case(1)
        call hecmw_precond_diag_11_clear()
      case(2)
        call hecmw_precond_diag_22_clear()
      case default
        call hecmw_precond_diag_nn_clear()
    end select
  end subroutine hecmw_precond_diag_clear
 
  subroutine hecmw_precond_diag_apply(ZP, NDOF)
    implicit none
    real(kind=kreal), intent(inout) :: zp(:)
    integer(kind=kint), intent(in) :: ndof
 
    select case(ndof)
      case(3)
        call hecmw_precond_diag_33_apply(zp)
      case(4)
        call hecmw_precond_diag_44_apply(zp)
      case(6)
        call hecmw_precond_diag_66_apply(zp)
      case(1)
        call hecmw_precond_diag_11_apply(zp)
      case(2)
        call hecmw_precond_diag_22_apply(zp)
      case default
        call hecmw_precond_diag_nn_apply(zp,ndof)
    end select
  end subroutine hecmw_precond_diag_apply
 
end module hecmw_precond_diag