hecmw_ordering_metis.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
! Copyright (c) 2020 FrontISTR Commons
! This software is released under the MIT License, see LICENSE.txt
!-------------------------------------------------------------------------------
 
!----------------------------------------------------------------------
!         ordering provided by Metis
!----------------------------------------------------------------------
module hecmw_ordering_metis
  use hecmw_util
  implicit none
 
  private
  public :: hecmw_ordering_metis_nodend
 
contains
 
  !======================================================================!
  !======================================================================!
  subroutine hecmw_ordering_metis_nodend(Neqns,Xadj,Adj0,Perm,Invp)
    implicit none
    !------
    integer(kind=kint), intent(in):: neqns
    integer(kind=kint), intent(in):: adj0(:)
    integer(kind=kint), intent(in):: xadj(:)
    integer(kind=kint), intent(out):: perm(:)
    integer(kind=kint), intent(out):: invp(:)
    !------
#ifdef HECMW_WITH_METIS
 
#  if HECMW_METIS_VER == 5
 
    integer(kind=kint), allocatable:: vwght(:)
    integer(kind=kint):: options(40)
    integer(kind=kint):: ierror
 
    allocate(vwght(neqns),stat=ierror)
    if ( ierror/=0 ) stop "ALLOCATION ERROR, vwght: SUB. gennd"
 
    vwght(:)=1
 
    call metis_setdefaultoptions(options)
    ! set fortran numbering
    options(18)=1
 
    call metis_nodend(neqns,xadj,adj0,vwght,options,perm,invp)
 
    deallocate(vwght)
 
#  elif HECMW_METIS_VER == 4
 
    integer(kind=kint):: numflag
    integer(kind=kint):: options(8)
 
    numflag=1
    options(:)=0
 
    call metis_nodend(neqns,xadj,adj0,numflag,options,perm,invp)
 
#  else
#    error unknown HECMW_METIS_VER
#  endif
 
#else
    stop "METIS not available"
#endif
  end subroutine hecmw_ordering_metis_nodend
 
end module hecmw_ordering_metis