hecmw_precond_RIF.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_rif
  use hecmw_util
  use hecmw_precond_rif_33
  use hecmw_precond_rif_nn
  implicit none
 
  private
  public :: hecmw_precond_rif_setup
  public :: hecmw_precond_rif_clear
  public :: hecmw_precond_rif_apply
 
contains
 
  subroutine hecmw_precond_rif_setup(hecMAT)
    implicit none
    type (hecmwst_matrix), intent(inout) :: hecmat
 
    select case(hecmat%NDOF)
      case(3)
        call hecmw_precond_rif_33_setup(hecmat)
      case default
        call hecmw_precond_rif_nn_setup(hecmat)
    end select
  end subroutine hecmw_precond_rif_setup
 
  subroutine hecmw_precond_rif_clear(NDOF)
    implicit none
    integer(kind=kint), intent(in) :: ndof
 
    select case(ndof)
      case(3)
        call hecmw_precond_rif_33_clear()
      case default
        call hecmw_precond_rif_nn_clear()
    end select
  end subroutine hecmw_precond_rif_clear
 
  subroutine hecmw_precond_rif_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_rif_33_apply(zp)
      case default
        call hecmw_precond_rif_nn_apply(zp,ndof)
    end select
  end subroutine hecmw_precond_rif_apply
 
end module hecmw_precond_rif