compress_nnz.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	compress_c_nnz (array, len)
subroutine	compress_i_nnz (array, len)
subroutine	compress_r_nnz (array, len)

Function/Subroutine Documentation

compress_c_nnz()

subroutine compress_c_nnz	(	integer, dimension(len), intent(in)	array,
		integer, intent(in)	len )

Definition at line 28 of file compress_nnz.F.

        IMPLICIT NONE
C---------------------------------------------
C         Arguments
C---------------------------------------------
        INTEGER, INTENT(IN) :: LEN
        INTEGER, INTENT(IN) :: ARRAY(LEN)
C      INTEGER, INTENT(INOUT) :: COMPR(*)
C---------------------------------------------
C        Local Variables
C---------------------------------------------
        INTEGER I,J,IBEGIN, NNZ
        LOGICAL continue_loop
 
 
        ibegin = 1
        j = 1
        DO WHILE (ibegin <= len)
          i = ibegin
          nnz = 0
          continue_loop = (array(i) /= ichar(' '))
          DO WHILE(continue_loop)
            IF(i <= len) THEN
              IF(array(i) /= ichar(' ')) THEN
                continue_loop = .true.
                nnz = nnz + 1
                i = i + 1
              ELSE
                continue_loop = .false.
              ENDIF
            ELSE
              continue_loop = .false.
            ENDIF
          ENDDO
          IF(nnz > 0) THEN
C            COMPR(J) = IBEGIN
C            COMPR(J+1) = NNZ
            CALL write_i_c(ibegin,1)
            CALL write_i_c(nnz,1)
            j = j + 2
            CALL write_i_c(array(ibegin:ibegin+nnz-1),nnz)
            j = j + nnz
          ENDIF
          ibegin = i + 1
        ENDDO
C        COMPR(J) = -1
        i = -1
        CALL write_i_c(i,1)
C        WRITE(6,*) LEN,"compressed into ",J, dble(J)/dble(LEN)

compress_i_nnz()

subroutine compress_i_nnz	(	integer, dimension(len), intent(in)	array,
		integer, intent(in)	len )

Definition at line 90 of file compress_nnz.F.

        IMPLICIT NONE
C---------------------------------------------
C         Arguments
C---------------------------------------------
        INTEGER, INTENT(IN) :: LEN
        INTEGER, INTENT(IN) :: ARRAY(LEN)
C      INTEGER, INTENT(INOUT) :: COMPR(*)
C---------------------------------------------
C        Local Variables
C---------------------------------------------
        INTEGER I,J,IBEGIN, NNZ
        LOGICAL continue_loop
 
        ibegin = 1
        j = 1
        DO WHILE (ibegin <= len)
          i = ibegin
          nnz = 0
          continue_loop = (array(i) /= 0)
          DO WHILE(continue_loop)
            IF(i <= len) THEN
              IF(array(i) /= 0) THEN
                continue_loop = .true.
                nnz = nnz + 1
                i = i + 1
              ELSE
                continue_loop = .false.
              ENDIF
            ELSE
              continue_loop = .false.
            ENDIF
          ENDDO
          IF(nnz > 0) THEN
C            COMPR(J) = IBEGIN
C            COMPR(J+1) = NNZ
            CALL write_i_c(ibegin,1)
            CALL write_i_c(nnz,1)
            j = j + 2
            CALL write_i_c(array(ibegin:ibegin+nnz-1),nnz)
            j = j + nnz
          ENDIF
          ibegin = i + 1
        ENDDO
C        COMPR(J) = -1
        i = -1
        CALL write_i_c(i,1)
C        WRITE(6,*) LEN,"I compressed into ",J, dble(J)/dble(LEN)

compress_r_nnz()

subroutine compress_r_nnz	(	dimension(len), intent(in)	array,
		integer, intent(in)	len )

Definition at line 153 of file compress_nnz.F.

#include      "implicit_f.inc"
C---------------------------------------------
C         Arguments
C---------------------------------------------
        INTEGER, INTENT(IN) :: LEN
        my_real , INTENT(IN) :: array(len)
C      INTEGER, INTENT(INOUT) :: COMPR(*)
C---------------------------------------------
C        Local Variables
C---------------------------------------------
        INTEGER I,J,IBEGIN, NNZ
        LOGICAL continue_loop
 
        ibegin = 1
        j = 1
        DO WHILE (ibegin <= len)
          i = ibegin
          nnz = 0
          continue_loop = (array(i) /= 0)
          DO WHILE(continue_loop)
            IF(i <= len) THEN
              IF(array(i) /= 0) THEN
                continue_loop = .true.
                nnz = nnz + 1
                i = i + 1
              ELSE
                continue_loop = .false.
              ENDIF
            ELSE
              continue_loop = .false.
            ENDIF
          ENDDO
          IF(nnz > 0) THEN
C            COMPR(J) = IBEGIN
C            COMPR(J+1) = NNZ
            CALL write_i_c(ibegin,1)
            CALL write_i_c(nnz,1)
            j = j + 2
            CALL write_db(array(ibegin:ibegin+nnz-1),nnz)
            j = j + nnz
          ENDIF
          ibegin = i + 1
        ENDDO
C        COMPR(J) = -1
        i = -1
        CALL write_i_c(i,1)
C        WRITE(6,*) LEN,"R compressed into ",J, dble(J)/dble(LEN)