ana_set_ordering.F File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	mumps_set_ordering (n, keep, sym, nprocs, iord, nbqd, avgdens, prok, mp)

Function/Subroutine Documentation

mumps_set_ordering()

subroutine mumps_set_ordering	(	integer, intent(in)	n,
		integer, dimension(500), intent(in)	keep,
		integer, intent(in)	sym,
		integer, intent(in)	nprocs,
		integer, intent(inout)	iord,
		integer, intent(in)	nbqd,
		integer, intent(in)	avgdens,
		logical, intent(in)	prok,
		integer, intent(in)	mp )

Definition at line 14 of file ana_set_ordering.F.

      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