fdurastante
/
psblas3
mirror of https://github.com/sfilippone/psblas3.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Fix GEPRT

Browse Source
repack-newsolve
Salvatore Filippone 4 months ago
parent 5430ba0e22
commit f10c6c1822
4 changed files with 16 additions and 12 deletions
Show Stats Download Patch File Download Diff File

5
base/serial/psb_cgeprt.f90
Unescape Escape View File

@ -29,7 +29,7 @@
! POSSIBILITY OF SUCH DAMAGE. ! POSSIBILITY OF SUCH DAMAGE.
! !
! !
! File: psb_scsprt.f90 ! File: psb_geprt.f90
! Subroutine: ! Subroutine:
! Arguments: ! Arguments:
@ -131,6 +131,7 @@ subroutine psb_cgeprt2(iout,a,head)
ncol = size(a,2) ncol = size(a,2)
write(iout,*) nrow,ncol write(iout,*) nrow,ncol
write(frmtv,'(a,i3.3,a)') '(',ncol,'2(es26.18,1x))' write(frmtv,'(a,i3.3,a)') '(',ncol,'2(es26.18,1x))'
do i=1,nrow do i=1,nrow
@ -161,7 +162,7 @@ subroutine psb_cgeprt1(iout,a,head)
write(iout,'(a)') '% ' write(iout,'(a)') '% '
nrow = size(a,1) nrow = size(a,1)
ncol = 1 ncol = 1
write(iout,*) nrow write(iout,*) nrow,ncol
write(frmtv,'(a,i3.3,a)') '(',ncol,'2(es26.18,1x))' write(frmtv,'(a,i3.3,a)') '(',ncol,'2(es26.18,1x))'

9
base/serial/psb_dgeprt.f90
Unescape Escape View File

@ -29,7 +29,7 @@
! POSSIBILITY OF SUCH DAMAGE. ! POSSIBILITY OF SUCH DAMAGE.
! !
! !
! File: psb_dcsprt.f90 ! File: psb_geprt.f90
! Subroutine: ! Subroutine:
! Arguments: ! Arguments:
@ -124,13 +124,14 @@ subroutine psb_dgeprt2(iout,a,head)
character(len=80) :: frmtv character(len=80) :: frmtv
integer(psb_ipk_) :: irs,ics,i,j, nmx, ni, nrow, ncol integer(psb_ipk_) :: irs,ics,i,j, nmx, ni, nrow, ncol
write(iout,'(a)') '%%MatrixMarket matrix array real general' write(iout,'(a)') '%%MatrixMarket matrix array complex general'
write(iout,'(a)') '% '//trim(head) write(iout,'(a)') '% '//trim(head)
write(iout,'(a)') '% ' write(iout,'(a)') '% '
nrow = size(a,1) nrow = size(a,1)
ncol = size(a,2) ncol = size(a,2)
write(iout,*) nrow,ncol write(iout,*) nrow,ncol
write(frmtv,'(a,i3.3,a)') '(',ncol,'(es26.18,1x))' write(frmtv,'(a,i3.3,a)') '(',ncol,'(es26.18,1x))'
do i=1,nrow do i=1,nrow
@ -156,12 +157,12 @@ subroutine psb_dgeprt1(iout,a,head)
character(len=80) :: frmtv character(len=80) :: frmtv
integer(psb_ipk_) :: irs,ics,i,j, nmx, ni, nrow, ncol integer(psb_ipk_) :: irs,ics,i,j, nmx, ni, nrow, ncol
write(iout,'(a)') '%%MatrixMarket matrix array real general' write(iout,'(a)') '%%MatrixMarket matrix array complex general'
write(iout,'(a)') '% '//trim(head) write(iout,'(a)') '% '//trim(head)
write(iout,'(a)') '% ' write(iout,'(a)') '% '
nrow = size(a,1) nrow = size(a,1)
ncol = 1 ncol = 1
write(iout,*) nrow write(iout,*) nrow,ncol
write(frmtv,'(a,i3.3,a)') '(',ncol,'(es26.18,1x))' write(frmtv,'(a,i3.3,a)') '(',ncol,'(es26.18,1x))'

9
base/serial/psb_sgeprt.f90
Unescape Escape View File

@ -29,7 +29,7 @@
! POSSIBILITY OF SUCH DAMAGE. ! POSSIBILITY OF SUCH DAMAGE.
! !
! !
! File: psb_scsprt.f90 ! File: psb_geprt.f90
! Subroutine: ! Subroutine:
! Arguments: ! Arguments:
@ -124,13 +124,14 @@ subroutine psb_sgeprt2(iout,a,head)
character(len=80) :: frmtv character(len=80) :: frmtv
integer(psb_ipk_) :: irs,ics,i,j, nmx, ni, nrow, ncol integer(psb_ipk_) :: irs,ics,i,j, nmx, ni, nrow, ncol
write(iout,'(a)') '%%MatrixMarket matrix array real general' write(iout,'(a)') '%%MatrixMarket matrix array complex general'
write(iout,'(a)') '% '//trim(head) write(iout,'(a)') '% '//trim(head)
write(iout,'(a)') '% ' write(iout,'(a)') '% '
nrow = size(a,1) nrow = size(a,1)
ncol = size(a,2) ncol = size(a,2)
write(iout,*) nrow,ncol write(iout,*) nrow,ncol
write(frmtv,'(a,i3.3,a)') '(',ncol,'(es26.18,1x))' write(frmtv,'(a,i3.3,a)') '(',ncol,'(es26.18,1x))'
do i=1,nrow do i=1,nrow
@ -156,12 +157,12 @@ subroutine psb_sgeprt1(iout,a,head)
character(len=80) :: frmtv character(len=80) :: frmtv
integer(psb_ipk_) :: irs,ics,i,j, nmx, ni, nrow, ncol integer(psb_ipk_) :: irs,ics,i,j, nmx, ni, nrow, ncol
write(iout,'(a)') '%%MatrixMarket matrix array real general' write(iout,'(a)') '%%MatrixMarket matrix array complex general'
write(iout,'(a)') '% '//trim(head) write(iout,'(a)') '% '//trim(head)
write(iout,'(a)') '% ' write(iout,'(a)') '% '
nrow = size(a,1) nrow = size(a,1)
ncol = 1 ncol = 1
write(iout,*) nrow write(iout,*) nrow,ncol
write(frmtv,'(a,i3.3,a)') '(',ncol,'(es26.18,1x))' write(frmtv,'(a,i3.3,a)') '(',ncol,'(es26.18,1x))'

5
base/serial/psb_zgeprt.f90
Unescape Escape View File

@ -29,7 +29,7 @@
! POSSIBILITY OF SUCH DAMAGE. ! POSSIBILITY OF SUCH DAMAGE.
! !
! !
! File: psb_scsprt.f90 ! File: psb_geprt.f90
! Subroutine: ! Subroutine:
! Arguments: ! Arguments:
@ -131,6 +131,7 @@ subroutine psb_zgeprt2(iout,a,head)
ncol = size(a,2) ncol = size(a,2)
write(iout,*) nrow,ncol write(iout,*) nrow,ncol
write(frmtv,'(a,i3.3,a)') '(',ncol,'2(es26.18,1x))' write(frmtv,'(a,i3.3,a)') '(',ncol,'2(es26.18,1x))'
do i=1,nrow do i=1,nrow
@ -161,7 +162,7 @@ subroutine psb_zgeprt1(iout,a,head)
write(iout,'(a)') '% ' write(iout,'(a)') '% '
nrow = size(a,1) nrow = size(a,1)
ncol = 1 ncol = 1
write(iout,*) nrow write(iout,*) nrow,ncol
write(frmtv,'(a,i3.3,a)') '(',ncol,'2(es26.18,1x))' write(frmtv,'(a,i3.3,a)') '(',ncol,'2(es26.18,1x))'

Write Preview
Loading…
Cancel
Save