Fixed usage of psb_spcnv.

psblas3-type-indexed
Salvatore Filippone 18 years ago
parent 6fdc3c2c0e
commit 031c0a6007

@ -124,7 +124,7 @@ contains
if (info == 0) read (infile,fmt=rhsfmt) (b(i),i=1,nrow) if (info == 0) read (infile,fmt=rhsfmt) (b(i),i=1,nrow)
endif endif
call psb_ipcsr2coo(a,ircode) call psb_spcnv(a,ircode,afmt='csr')
if (ircode /= 0) goto 993 if (ircode /= 0) goto 993
call psb_sp_reall(a,2*nnzero,ircode) call psb_sp_reall(a,2*nnzero,ircode)
@ -139,7 +139,7 @@ contains
end if end if
end do end do
a%infoa(psb_nnz_) = nzr a%infoa(psb_nnz_) = nzr
call psb_ipcoo2csr(a,ircode) call psb_spcnv(a,ircode,afmt='csr')
if (ircode /= 0) goto 993 if (ircode /= 0) goto 993
else else
@ -384,7 +384,7 @@ contains
endif endif
call psb_ipcsr2coo(a,ircode) call psb_spcnv(a,ircode,afmt='coo')
if (ircode /= 0) then if (ircode /= 0) then
write(0,*) 'ipcsr2coo ',ircode write(0,*) 'ipcsr2coo ',ircode
goto 993 goto 993
@ -408,7 +408,7 @@ contains
end if end if
end do end do
a%infoa(psb_nnz_) = nzr a%infoa(psb_nnz_) = nzr
call psb_ipcoo2csr(a,ircode) call psb_spcnv(a,ircode,afmt='csr')
if (ircode /= 0) then if (ircode /= 0) then
write(0,*) 'ipcoo2csr ',ircode write(0,*) 'ipcoo2csr ',ircode
goto 993 goto 993
@ -439,7 +439,7 @@ contains
if (info == 0) read (infile,fmt=rhsfmt) (b(i),i=1,nrow) if (info == 0) read (infile,fmt=rhsfmt) (b(i),i=1,nrow)
endif endif
call psb_ipcsr2coo(a,ircode) call psb_spcnv(a,ircode,afmt='coo')
if (ircode /= 0) then if (ircode /= 0) then
write(0,*) 'ipcsr2coo ',ircode write(0,*) 'ipcsr2coo ',ircode
goto 993 goto 993
@ -462,7 +462,7 @@ contains
end if end if
end do end do
a%infoa(psb_nnz_) = nzr a%infoa(psb_nnz_) = nzr
call psb_ipcoo2csr(a,ircode) call psb_spcnv(a,ircode,afmt='csr')
if (ircode /= 0) then if (ircode /= 0) then
write(0,*) 'ipcoo2csr ',ircode write(0,*) 'ipcoo2csr ',ircode
goto 993 goto 993

@ -100,7 +100,7 @@ contains
read(infile,fmt=*,end=902) a%ia1(i),a%ia2(i),a%aspk(i) read(infile,fmt=*,end=902) a%ia1(i),a%ia2(i),a%aspk(i)
end do end do
a%infoa(psb_nnz_) = nnzero a%infoa(psb_nnz_) = nnzero
call psb_ipcoo2csr(a,ircode) call psb_spcnv(a,ircode,afmt='csr')
else if ((tolower(type) == 'real').and.(tolower(sym) == 'symmetric')) then else if ((tolower(type) == 'real').and.(tolower(sym) == 'symmetric')) then
! we are generally working with non-symmetric matrices, so ! we are generally working with non-symmetric matrices, so
@ -123,7 +123,7 @@ contains
end if end if
end do end do
a%infoa(psb_nnz_) = nzr a%infoa(psb_nnz_) = nzr
call psb_ipcoo2csr(a,ircode) call psb_spcnv(a,ircode,afmt='csr')
else else
write(0,*) 'read_matrix: matrix type not yet supported' write(0,*) 'read_matrix: matrix type not yet supported'
@ -259,7 +259,7 @@ contains
end do end do
a%infoa(psb_nnz_) = nnzero a%infoa(psb_nnz_) = nnzero
call psb_ipcoo2csr(a,ircode) call psb_spcnv(a,ircode,afmt='csr')
else if ((tolower(type) == 'complex').and.(tolower(sym) == 'symmetric')) then else if ((tolower(type) == 'complex').and.(tolower(sym) == 'symmetric')) then
! we are generally working with non-symmetric matrices, so ! we are generally working with non-symmetric matrices, so
@ -283,7 +283,7 @@ contains
end if end if
end do end do
a%infoa(psb_nnz_) = nzr a%infoa(psb_nnz_) = nzr
call psb_ipcoo2csr(a,ircode) call psb_spcnv(a,ircode,afmt='csr')
else if ((tolower(type) == 'complex').and.(tolower(sym) == 'hermitian')) then else if ((tolower(type) == 'complex').and.(tolower(sym) == 'hermitian')) then
! we are generally working with non-symmetric matrices, so ! we are generally working with non-symmetric matrices, so
@ -307,7 +307,7 @@ contains
end if end if
end do end do
a%infoa(psb_nnz_) = nzr a%infoa(psb_nnz_) = nzr
call psb_ipcoo2csr(a,ircode) call psb_spcnv(a,ircode,afmt='csr')
else else
write(0,*) 'read_matrix: matrix type not yet supported' write(0,*) 'read_matrix: matrix type not yet supported'

Loading…
Cancel
Save