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
56
57
58
59
60
61
62
63
64
65
66
67
68
|
!
! Copyright (c) 2003, The Regents of the University of California, through
! Lawrence Berkeley National Laboratory (subject to receipt of any required
! approvals from U.S. Dept. of Energy)
!
! All rights reserved.
!
! The source code is distributed under BSD license, see the file License.txt
! at the top-level directory.
!
!#include "superlu_config.fh"
program f_main
integer maxn, maxnz
parameter ( maxn = 10000, maxnz = 100000 )
#if (XSDK_INDEX_SIZE==64)
integer*8 rowind(maxnz), colptr(maxn+1)
#else
integer rowind(maxnz), colptr(maxn+1)
#endif
real*8 values(maxnz), b(maxn)
integer n, nnz, nrhs, ldb, info, iopt
integer*8 factors
!
call hbcode1(n, n, nnz, values, rowind, colptr)
! write(I6, *) (rowind(i), i=1,nnz)
! write(I6, *) (colptr(i), i=1,n+1)
!
nrhs = 1
ldb = n
do i = 1, n
b(i) = 1
enddo
!
! First, factorize the matrix. The factors are stored in *factors* handle.
iopt = 1
call c_fortran_dgssv( iopt, n, nnz, nrhs, values, rowind, colptr, &
b, ldb, factors, info )
!
if (info .eq. 0) then
write (*,*) 'Factorization succeeded'
else
write(*,*) 'INFO from factorization = ', info
endif
!
! Second, solve the system using the existing factors.
iopt = 2
call c_fortran_dgssv( iopt, n, nnz, nrhs, values, rowind, colptr, &
b, ldb, factors, info )
!
if (info .eq. 0) then
write (*,*) 'Solve succeeded'
write (*,*) (b(i), i=1, 10)
else
write(*,*) 'INFO from triangular solve = ', info
endif
! Last, free the storage allocated inside SuperLU
iopt = 3
call c_fortran_dgssv( iopt, n, nnz, nrhs, values, rowind, colptr, &
b, ldb, factors, info )
!
stop
end
|
