hecmw_precond_BILU.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_bilu
  use hecmw_util
  use hecmw_precond_bilu_33
  use hecmw_precond_bilu_44
  use hecmw_precond_bilu_66
  use hecmw_precond_bilu_nn
  implicit none
 
  private
  public :: hecmw_precond_bilu_setup
  public :: hecmw_precond_bilu_clear
  public :: hecmw_precond_bilu_apply
 
contains
 
  subroutine hecmw_precond_bilu_setup(hecMAT)
    implicit none
    type (hecmwst_matrix), intent(inout) :: hecmat
 
    select case(hecmat%NDOF)
      case(3)
        call hecmw_precond_bilu_33_setup(hecmat)
      case(4)
        call hecmw_precond_bilu_44_setup(hecmat)
      case(6)
        call hecmw_precond_bilu_66_setup(hecmat)
      case default
        call hecmw_precond_bilu_nn_setup(hecmat)
    end select
  end subroutine hecmw_precond_bilu_setup
 
  subroutine hecmw_precond_bilu_clear(NDOF)
    implicit none
    integer(kind=kint), intent(in) :: ndof
 
    select case(ndof)
      case(3)
        call hecmw_precond_bilu_33_clear()
      case(4)
        call hecmw_precond_bilu_44_clear()
      case(6)
        call hecmw_precond_bilu_66_clear()
      case default
        call hecmw_precond_bilu_nn_clear()
    end select
  end subroutine hecmw_precond_bilu_clear
 
  subroutine hecmw_precond_bilu_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_bilu_33_apply(zp)
      case(4)
        call hecmw_precond_bilu_44_apply(zp)
      case(6)
        call hecmw_precond_bilu_66_apply(zp)
      case default
        call hecmw_precond_bilu_nn_apply(zp,ndof)
    end select
  end subroutine hecmw_precond_bilu_apply
 
end module hecmw_precond_bilu