Renamed pargen sample.

psblas3-type-indexed
Salvatore Filippone 17 years ago
parent e099c19498
commit d2ad1ed196

@ -37,7 +37,7 @@ program df_sample
implicit none implicit none
! input parameters ! input parameters
character(len=40) :: cmethd, ptype, mtrx_file, rhs_file character(len=40) :: kmethd, ptype, mtrx_file, rhs_file
! sparse matrices ! sparse matrices
type(psb_dspmat_type) :: a, aux_a type(psb_dspmat_type) :: a, aux_a
@ -95,7 +95,7 @@ program df_sample
! !
! get parameters ! get parameters
! !
call get_parms(ictxt,mtrx_file,rhs_file,cmethd,ptype,& call get_parms(ictxt,mtrx_file,rhs_file,kmethd,ptype,&
& ipart,afmt,istopc,itmax,itrace,irst,eps) & ipart,afmt,istopc,itmax,itrace,irst,eps)
call psb_barrier(ictxt) call psb_barrier(ictxt)
@ -214,7 +214,7 @@ program df_sample
iparm = 0 iparm = 0
call psb_barrier(ictxt) call psb_barrier(ictxt)
t1 = psb_wtime() t1 = psb_wtime()
call psb_krylov(cmethd,a,prec,b_col,x_col,eps,desc_a,info,& call psb_krylov(kmethd,a,prec,b_col,x_col,eps,desc_a,info,&
& itmax=itmax,iter=iter,err=err,itrace=itrace,istop=istopc,irst=irst) & itmax=itmax,iter=iter,err=err,itrace=itrace,istop=istopc,irst=irst)
call psb_barrier(ictxt) call psb_barrier(ictxt)
t2 = psb_wtime() - t1 t2 = psb_wtime() - t1

@ -35,11 +35,11 @@ contains
! !
! Get iteration parameters from the command line ! Get iteration parameters from the command line
! !
subroutine get_parms(ictxt,mtrx_file,rhs_file,cmethd,ptype,ipart,& subroutine get_parms(ictxt,mtrx_file,rhs_file,kmethd,ptype,ipart,&
& afmt,istopc,itmax,itrace,irst,eps) & afmt,istopc,itmax,itrace,irst,eps)
use psb_base_mod use psb_base_mod
integer :: ictxt integer :: ictxt
character(len=40) :: cmethd, mtrx_file, rhs_file, ptype character(len=40) :: kmethd, mtrx_file, rhs_file, ptype
integer :: iret, istopc,itmax,itrace,ipart,irst integer :: iret, istopc,itmax,itrace,ipart,irst
character(len=40) :: charbuf character(len=40) :: charbuf
real(kind(1.d0)) :: eps real(kind(1.d0)) :: eps
@ -54,14 +54,14 @@ contains
if (ip >= 5) then if (ip >= 5) then
read(*,*) mtrx_file read(*,*) mtrx_file
read(*,*) rhs_file read(*,*) rhs_file
read(*,*) cmethd read(*,*) kmethd
read(*,*) ptype read(*,*) ptype
read(*,*) afmt read(*,*) afmt
call psb_bcast(ictxt,mtrx_file) call psb_bcast(ictxt,mtrx_file)
call psb_bcast(ictxt,rhs_file) call psb_bcast(ictxt,rhs_file)
call psb_bcast(ictxt,cmethd) call psb_bcast(ictxt,kmethd)
call psb_bcast(ictxt,ptype) call psb_bcast(ictxt,ptype)
call psb_bcast(ictxt,afmt) call psb_bcast(ictxt,afmt)
@ -102,7 +102,7 @@ contains
write(*,'("Solving matrix : ",a)') mtrx_file write(*,'("Solving matrix : ",a)') mtrx_file
write(*,'("Number of processors : ",i3)') np write(*,'("Number of processors : ",i3)') np
write(*,'("Data distribution : ",i2)') ipart write(*,'("Data distribution : ",i2)') ipart
write(*,'("Iterative method : ",a)') cmethd write(*,'("Iterative method : ",a)') kmethd
write(*,'("Preconditioner : ",a)') ptype write(*,'("Preconditioner : ",a)') ptype
write(*,'("Restart parameter : ",i2)') irst write(*,'("Restart parameter : ",i2)') irst
write(*,'("Storage format : ",a)') afmt(1:3) write(*,'("Storage format : ",a)') afmt(1:3)
@ -116,7 +116,7 @@ contains
! Receive Parameters ! Receive Parameters
call psb_bcast(ictxt,mtrx_file) call psb_bcast(ictxt,mtrx_file)
call psb_bcast(ictxt,rhs_file) call psb_bcast(ictxt,rhs_file)
call psb_bcast(ictxt,cmethd) call psb_bcast(ictxt,kmethd)
call psb_bcast(ictxt,ptype) call psb_bcast(ictxt,ptype)
call psb_bcast(ictxt,afmt) call psb_bcast(ictxt,afmt)

@ -227,8 +227,7 @@ contains
integer :: ictxt integer :: ictxt
character(len=*) :: kmethd, ptype, afmt character(len=*) :: kmethd, ptype, afmt
integer :: idim, istopc,itmax,itrace,irst integer :: idim, istopc,itmax,itrace,irst
integer :: iargc, np, iam integer :: np, iam
external iargc
integer :: intbuf(10), ip integer :: intbuf(10), ip
call psb_info(ictxt, iam, np) call psb_info(ictxt, iam, np)

Loading…
Cancel
Save