srank_revealing.F File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	smumps_get_ns_options_facto (n, keep, icntl, mpg)
subroutine	smumps_get_ns_options_solve (icntl, keep, nrhs, mpg, info)
subroutine	smumps_rr_init_pointers (id)
subroutine	smumps_rr_free_pointers (id)

Function/Subroutine Documentation

smumps_get_ns_options_facto()

subroutine smumps_get_ns_options_facto	(	integer	n,
		integer, dimension(500)	keep,
		integer, dimension(60)	icntl,
		integer	mpg )

Definition at line 14 of file srank_revealing.F.

      IMPLICIT NONE
          INTEGER N, KEEP(500), ICNTL(60), MPG
          keep(19)=0
          RETURN

smumps_get_ns_options_solve()

subroutine smumps_get_ns_options_solve	(	integer, dimension(60), intent(in)	icntl,
		integer, dimension(500), intent(in)	keep,
		integer, intent(in)	nrhs,
		integer, intent(in)	mpg,
		integer, dimension(80), intent(inout)	info )

Definition at line 20 of file srank_revealing.F.

      IMPLICIT NONE
      INTEGER, intent(in)   :: KEEP(500), NRHS, MPG, ICNTL(60)
      INTEGER, intent(inout):: INFO(80)
      IF (keep(19).EQ.0.AND.keep(110).EQ.0) THEN
        IF (keep(111).NE.0) THEN
         info(1) = -37
         info(2) = 56  
         IF (keep(110).EQ.0) info(2) = 24 
          IF(mpg.GT.0) THEN
           WRITE( mpg,'(A)')
     &'** ERROR  : Null space computation requirement'
          WRITE( mpg,'(A)')
     &'** not consistent with factorization options'
         ENDIF
         GOTO 333
        ENDIF
      ENDIF
       IF (icntl(9).NE.1) THEN
         IF (keep(111).NE.0) THEN
          info(1) = -37
          info(2) = 9  
          IF (mpg.GT.0) THEN
           WRITE(mpg,'(A)')
     &'** error  icntl(25) incompatible with '
           WRITE( MPG,'(a)')
     &'** option transposed system(icntl(9).ne.1) '
          ENDIF
         ENDIF
         GOTO 333
       ENDIF
.EQ.      IF (KEEP(19)2) THEN
.NE..AND..EQ.       IF ((KEEP(111)0)(KEEP(50)0)) THEN
         INFO(1) = -37
         INFO(2) = 0  
.GT.         IF (MPG0) THEN
          WRITE(MPG,'(a)')
     &'** error  icntl(25) incompatible with '
          WRITE( MPG,'(a)')
     &'** option rrqr(icnlt(56)=2) and unsym. matrices '
         ENDIF
        ENDIF
        GOTO 333
      ENDIF
.eq..AND..NE.      IF (KEEP(111)-1NRHSKEEP(112)+KEEP(17))THEN
            INFO(1)=-32
            INFO(2)=NRHS
            GOTO 333
      ENDIF
.gt..AND..NE.      IF (KEEP(111)0  NRHS  1) THEN
            INFO(1)=-32
            INFO(2)=NRHS
            GOTO 333
      ENDIF
.NE..AND..NE.      IF (KEEP(248) 0KEEP(111)0) THEN
.GT.         IF (MPG0) THEN
            WRITE(MPG,'(a)')
     &           ' error: icntl(20) and icntl(30) functionalities ',
     &           ' incompatible with null space'
         ENDIF
         INFO(1) = -37
.NE.         IF (KEEP(237)0) THEN
            INFO(2) = 30    
.GT.            IF (MPG0) THEN
               WRITE(MPG,'(a)')
     &           ' error: icntl(30) functionality ',
     &              ' incompatible with null space'
            ENDIF
         ELSE
.GT.            IF (MPG0) THEN
               WRITE(MPG,'(a)')
     &              ' error: icntl(20) functionality ',
     &              ' incompatible with null space'
            ENDIF
            INFO(2) = 20    
         ENDIF
         GOTO 333
      ENDIF
.LT..OR.      IF (( KEEP(111)  -1 ) 
.GT..OR.     &     (KEEP(111)KEEP(112)+KEEP(17)) 
.EQ..AND..EQ.     &     (KEEP(111) -1  KEEP(112)+KEEP(17)0))
     &     THEN
         INFO(1)=-36
         INFO(2)=KEEP(111)
         GOTO 333
      ENDIF
.NE..AND..NE.      IF (KEEP(221)0KEEP(111)0) THEN
         INFO(1)=-37
         INFO(2)=26
         GOTO 333
      ENDIF
 333  CONTINUE
      RETURN

smumps_rr_free_pointers()

subroutine smumps_rr_free_pointers

(

type (smumps_struc)

id

)

Definition at line 125 of file srank_revealing.F.

      USE smumps_struc_def
      IMPLICIT NONE
      TYPE (SMUMPS_STRUC) id
      IF (associated(id%root%QR_TAU))  THEN
        DEALLOCATE(id%root%QR_TAU)
        NULLIFY(id%root%QR_TAU)
      ENDIF
      IF (associated(id%root%SVD_U))  THEN
        DEALLOCATE(id%root%SVD_U)
        NULLIFY(id%root%SVD_U)
      ENDIF
      IF (associated(id%root%SVD_VT))  THEN
        DEALLOCATE(id%root%SVD_VT)
        NULLIFY(id%root%SVD_VT)
      ENDIF
      IF (associated(id%root%SINGULAR_VALUES))  THEN
        DEALLOCATE(id%root%SINGULAR_VALUES)
        NULLIFY(id%root%SINGULAR_VALUES)
      ENDIF
      RETURN

smumps_rr_init_pointers()

subroutine smumps_rr_init_pointers

(

type (smumps_struc)

id

)

Definition at line 115 of file srank_revealing.F.

      USE smumps_struc_def
      IMPLICIT NONE
      TYPE (SMUMPS_STRUC) id
      NULLIFY(id%root%QR_TAU)
      NULLIFY(id%root%SVD_U)
      NULLIFY(id%root%SVD_VT)
      NULLIFY(id%root%SINGULAR_VALUES)
      RETURN