hecmw_jadm.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_jad_type
  use hecmw_util
  use m_hecmw_comm_f
  use hecmw_jad_type_33
  use hecmw_jad_type_44
  use hecmw_jad_type_nn
  implicit none
 
  private
 
  public :: hecmw_jad_init
  public :: hecmw_jad_finalize
  public :: hecmw_jad_is_initialized
  public :: hecmw_jad_matvec
 
  !C---------------------- AU&AL
  integer(kind=kint) :: INITIALIZED = 0
 
contains
 
  subroutine hecmw_jad_init(hecMAT)
    type(hecmwst_matrix) :: hecmat
    select case(hecmat%NDOF)
      case(3)
        call hecmw_jad_init_33(hecmat)
      case(4)
        call hecmw_jad_init_44(hecmat)
      case default
        call hecmw_jad_init_nn(hecmat)
    end select
    initialized = 1
  end subroutine hecmw_jad_init
 
  subroutine hecmw_jad_finalize(hecMAT)
    type(hecmwst_matrix) :: hecmat
 
    select case(hecmat%NDOF)
      case(3)
        call hecmw_jad_finalize_33()
      case(4)
        call hecmw_jad_finalize_44()
      case default
        call hecmw_jad_finalize_nn()
    end select
    initialized = 0
  end subroutine hecmw_jad_finalize
 
  function hecmw_jad_is_initialized()
    integer(kind=kint) :: hecmw_jad_is_initialized
    hecmw_jad_is_initialized = initialized
  end function hecmw_jad_is_initialized
 
  subroutine hecmw_jad_matvec(hecMESH, hecMAT, X, Y, COMMtime)
 
    type(hecmwst_local_mesh), intent(in) :: hecmesh
    type(hecmwst_matrix), intent(in), target :: hecmat
    real(kind=kreal), intent(in) :: x(:)
    real(kind=kreal), intent(out) :: y(:)
    real(kind=kreal), intent(inout) :: commtime
    select case(hecmat%NDOF)
      case(3)
        call hecmw_jad_matvec_33(hecmesh, hecmat, x, y, commtime)
      case(4)
        call hecmw_jad_matvec_44(hecmesh, hecmat, x, y, commtime)
      case default
        call hecmw_jad_matvec_nn(hecmesh, hecmat, x, y, commtime)
    end select
  end subroutine hecmw_jad_matvec
 
end module hecmw_jad_type