hecmw_precond_SAINV.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_sainv
  use hecmw_util
  use hecmw_precond_sainv_33
  use hecmw_precond_sainv_nn
  implicit none
 
  private
  public :: hecmw_precond_sainv_setup
  public :: hecmw_precond_sainv_clear
  public :: hecmw_precond_sainv_apply
 
contains
 
  subroutine hecmw_precond_sainv_setup(hecMAT)
    implicit none
    type (hecmwst_matrix), intent(inout) :: hecmat
 
    select case(hecmat%NDOF)
      case(3)
        call hecmw_precond_sainv_33_setup(hecmat)
      case default
        call hecmw_precond_sainv_nn_setup(hecmat)
    end select
  end subroutine hecmw_precond_sainv_setup
 
  subroutine hecmw_precond_sainv_clear(NDOF)
    implicit none
    integer(kind=kint), intent(in) :: ndof
 
    select case(ndof)
      case(3)
        call hecmw_precond_sainv_33_clear()
      case default
        call hecmw_precond_sainv_nn_clear()
    end select
  end subroutine hecmw_precond_sainv_clear
 
  subroutine hecmw_precond_sainv_apply(R, ZP, NDOF)
    implicit none
    real(kind=kreal), intent(in) :: r(:)
    real(kind=kreal), intent(inout) :: zp(:)
    integer(kind=kint), intent(in) :: ndof
 
    select case(ndof)
      case(3)
        call hecmw_precond_sainv_33_apply(r,zp)
      case default
        call hecmw_precond_sainv_nn_apply(r,zp)
    end select
  end subroutine hecmw_precond_sainv_apply
 
end module hecmw_precond_sainv