1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
|
#include "superlu_config.fh"
subroutine hbcode1(nrow, ncol, nnzero, values, rowind, colptr)
! ================================================================
! ... SAMPLE CODE FOR READING A SPARSE MATRIX IN STANDARD FORMAT
! ================================================================
CHARACTER TITLE*72 , KEY*8 , MXTYPE*3 , &
PTRFMT*16, INDFMT*16, VALFMT*20, RHSFMT*20
INTEGER TOTCRD, PTRCRD, INDCRD, VALCRD, RHSCRD, &
NROW , NCOL , NNZERO, NELTVL
#if (XSDK_INDEX_SIZE==64)
INTEGER*8 COLPTR (*), ROWIND (*)
#else
INTEGER COLPTR (*), ROWIND (*)
#endif
REAL*8 VALUES (*)
! ------------------------
! ... READ IN HEADER BLOCK
! ------------------------
READ ( *, 1000 ) TITLE , KEY , &
TOTCRD, PTRCRD, INDCRD, VALCRD, RHSCRD, &
MXTYPE, NROW , NCOL , NNZERO, NELTVL, &
PTRFMT, INDFMT, VALFMT, RHSFMT
1000 FORMAT ( A72, A8 / 5I14 / A3, 11X, 4I14 / 2A16, 2A20 )
! -------------------------
! ... READ MATRIX STRUCTURE
! -------------------------
READ ( *, PTRFMT ) ( COLPTR (I), I = 1, NCOL+1 )
do i = 1,ncol+1
write(I5, *) colptr(i)
end do
READ ( *, INDFMT ) ( ROWIND (I), I = 1, NNZERO )
IF ( VALCRD .GT. 0 ) THEN
! ----------------------
! ... READ MATRIX VALUES
! ----------------------
READ ( *, VALFMT ) ( VALUES (I), I = 1, NNZERO )
ENDIF
return
end
|
