ana_set_ordering.F

Go to the documentation of this file.

C
C  This file is part of MUMPS 5.5.1, released
C  on Tue Jul 12 13:17:24 UTC 2022
C
C
C  Copyright 1991-2022 CERFACS, CNRS, ENS Lyon, INP Toulouse, Inria,
C  Mumps Technologies, University of Bordeaux.
C
C  This version of MUMPS is provided to you free of charge. It is
C  released under the CeCILL-C license 
C  (see doc/CeCILL-C_V1-en.txt, doc/CeCILL-C_V1-fr.txt, and
C  https://cecill.info/licences/Licence_CeCILL-C_V1-en.html)
C
      SUBROUTINE mumps_set_ordering(N, KEEP, SYM, NPROCS, IORD,
     &                       NBQD, AvgDens,
     &                       PROK, MP)
      IMPLICIT NONE
      INTEGER, intent(in)    :: N, KEEP(500), NPROCS, SYM
      INTEGER, intent(in)    :: NBQD, AvgDens, MP
      LOGICAL, intent(in)    :: PROK
      INTEGER, intent(inout) :: IORD
      INTEGER MAXQD
      parameter(maxqd=2)
      INTEGER SMALLSYM, SMALLUNS
      parameter(smalluns=5000, smallsym=10000)
#if ! defined(metis) && ! defined(parmetis) && ! defined(metis4) && ! defined(parmetis3)
      IF ( iord .EQ. 5 ) THEN
        IF (prok) WRITE(mp,*)
     &  'WARNING: METIS not available. Ordering set to default.'
        iord = 7
      END IF
#endif
#if ! defined(pord)
      IF ( iord .EQ. 4 ) THEN
        IF (prok) WRITE(mp,*)
     &  'WARNING: PORD not available. Ordering set to default.'
        iord = 7
      END IF
#endif
#if ! defined(scotch) && !  defined(ptscotch)
      IF ( iord .EQ. 3 ) THEN
        IF (prok) WRITE(mp,*)
     &  'WARNING: SCOTCH not available. Ordering set to default.'
        iord = 7
      END IF
#endif
      IF (iord.EQ.7) THEN
        IF (sym.NE.0) THEN
          IF ( n.LE.smallsym ) THEN 
             IF (nbqd.GE.maxqd) THEN
               iord = 6         
             ELSE
               iord = 2         
             ENDIF
          ELSE
#if  defined(metis) || defined(parmetis) || defined(metis4) || defined(parmetis3)
             iord = 5
#else
#  if defined(scotch) || defined(ptscotch)
             iord = 3
#  else
#    if defined(pord)
               iord = 4
#    else
               iord = 6
#    endif
#  endif
#endif
          ENDIF
        ELSE
          IF ( n.LE.smalluns ) THEN
            IF (nbqd.GE.maxqd) THEN
              iord = 6  
            ELSE
              iord = 2   
            ENDIF
          ELSE
#if  defined(metis) || defined(parmetis) || defined(metis4) || defined(parmetis3)
            iord = 5
#else
#  if defined(scotch) || defined(ptscotch)
            iord = 3
#  else
#    if defined(pord)
              iord = 4
#    else
              iord = 6
#    endif
#  endif
#endif
          ENDIF
        ENDIF
      ENDIF
      RETURN
      END SUBROUTINE mumps_set_ordering