Category : Miscellaneous Language Source Code
Archive   : MATBYU.ZIP
Filename : PRINTS.FOR

C=======================================================================
C...THIS SUBROUTINE PRINTS A MATRIX IN SPARSE FORMAT.
C=======================================================================

SUBROUTINE PRINTS (MATNAM,DOUB,INTE,NREC,NPG,NLIN)

C...ARGUMENTS
DOUBLE PRECISION DOUB(*)
INTEGER NREC,NPG,NLIN,INTE(*)
CHARACTER MATNAM(*)*6

C...LOCAL VARIABLES
INTEGER N,M,IREC,I
CHARACTER NAME*6
LOGICAL DUP,INT

C-----------------------------------------------------------------------

C...PRINT A SPARSE MATRIX

C...READ THE MATRIX NAME
WRITE(*,*)'ENTER NAME OF THE SPARSE MATRIX TO BE PRINTED'
READ(*,1017)NAME

C...FIND THE LOCATION OF THE MATRIX IN THE LIBRARY
CALL DUPCHK(NAME,MATNAM,DUP,IREC,NREC)
IF (DUP) THEN

C...CHECK FOR AN INTEGER MATRIX
CALL INTCHK (NAME,INT)

C...READ AN INTEGER MATRIX
IF (INT) THEN
READ(1,REC=IREC+1)NAME,N,M,(INTE(I),I=1,N*M)
KROW=1+N*M
KCOL=KROW+N*M
KVAL=KCOL+N*M

C...PRINT AN INTEGER MATRIX
CALL ISPRNT (INTE(1),N,M,NAME,NPG,NLIN,INTE(KROW),
+ INTE(KCOL),INTE(KVAL))

C...READ A FLOATING POINT MATRIX
ELSE
READ(1,REC=IREC+1)MATNAM(IREC),N,M,(DOUB(I),I=1,N*M)
KROW=1+N*M
KCOL=KROW+N*M
KVAL=1+N*M

C...PRINT A FLOATING POINT MATRIX
CALL SPRINT (DOUB(1),N,M,NAME,NPG,NLIN,INTE(KROW),
+ INTE(KCOL),DOUB(KVAL))
END IF

C...REPORT THAT MATRIX WAS NOT IN LIBRARY FILE
ELSE
WRITE(*,*)'MATRIX NOT IN LIBRARY FILE'
RETURN
END IF

C...RETURN TO MATRIX I/O MENU
RETURN

C-----------------------------------------------------------------------
C...FORMATS
1017 FORMAT(A6)
C-----------------------------------------------------------------------

END