hecmw_precond_SSOR.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_ssor
  use hecmw_util
  use hecmw_precond_ssor_11
  use hecmw_precond_ssor_22
  use hecmw_precond_ssor_33
  use hecmw_precond_ssor_44
  use hecmw_precond_ssor_66
  use hecmw_precond_ssor_nn
  implicit none
 
  private
  public :: hecmw_precond_ssor_setup
  public :: hecmw_precond_ssor_clear
  public :: hecmw_precond_ssor_apply
 
contains
 
  subroutine hecmw_precond_ssor_setup(hecMAT)
    implicit none
    type (hecmwst_matrix), intent(inout) :: hecmat
 
    select case(hecmat%NDOF)
      case(3)
        call hecmw_precond_ssor_33_setup(hecmat)
      case(4)
        call hecmw_precond_ssor_44_setup(hecmat)
      case(6)
        call hecmw_precond_ssor_66_setup(hecmat)
      case(1)
        call hecmw_precond_ssor_11_setup(hecmat)
      case(2)
        call hecmw_precond_ssor_22_setup(hecmat)
      case default
        call hecmw_precond_ssor_nn_setup(hecmat)
    end select
  end subroutine hecmw_precond_ssor_setup
 
  subroutine hecmw_precond_ssor_clear(hecMAT)
    implicit none
    type (hecmwst_matrix), intent(inout) :: hecmat
 
    select case(hecmat%NDOF)
      case(3)
        call hecmw_precond_ssor_33_clear(hecmat)
      case(4)
        call hecmw_precond_ssor_44_clear(hecmat)
      case(6)
        call hecmw_precond_ssor_66_clear(hecmat)
      case(1)
        call hecmw_precond_ssor_11_clear(hecmat)
      case(2)
        call hecmw_precond_ssor_22_clear(hecmat)
      case default
        call hecmw_precond_ssor_nn_clear(hecmat)
    end select
  end subroutine hecmw_precond_ssor_clear
 
  subroutine hecmw_precond_ssor_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_ssor_33_apply(zp)
      case(4)
        call hecmw_precond_ssor_44_apply(zp)
      case(6)
        call hecmw_precond_ssor_66_apply(zp)
      case(1)
        call hecmw_precond_ssor_11_apply(zp)
      case(2)
        call hecmw_precond_ssor_22_apply(zp)
      case default
        call hecmw_precond_ssor_nn_apply(zp,ndof)
    end select
  end subroutine hecmw_precond_ssor_apply
 
end module hecmw_precond_ssor