predgmres.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	predgmres (n, b, x, nelt, ia, ja, a, isym, itol, tol, itmax, iter, err, ierr, iunit, sb, sx, rgwk, lrgw, igwk, ligw, rwork, iwork)

Function/Subroutine Documentation

predgmres()

subroutine predgmres	(	integer	n,
		real*8, dimension(*)	b,
		real*8, dimension(*)	x,
		integer	nelt,
		integer, dimension(*)	ia,
		integer, dimension(*)	ja,
		real*8, dimension(*)	a,
		integer	isym,
		integer	itol,
		real*8	tol,
		integer	itmax,
		integer	iter,
		real*8	err,
		integer	ierr,
		integer	iunit,
		real*8, dimension(*)	sb,
		real*8, dimension(*)	sx,
		real*8, dimension(*)	rgwk,
		integer	lrgw,
		integer, dimension(*)	igwk,
		integer	ligw,
		real*8, dimension(*)	rwork,
		integer, dimension(*)	iwork
	)

!
      implicit none
!
      integer n,nelt,ia(*),ja(*),isym,itol,itmax,iter,ierr,
     &  iunit,lrgw,igwk(*),ligw,iwork(*)
!
      real*8 b(*),x(*),a(*),tol,err,sb(*),sx(*),rgwk(*),
     &  rwork(*)
!
      external matvec,msolve
!
      itol=0
      tol=1.e-6
      itmax=0
      iunit=0
!
      igwk(1)=10
      igwk(2)=10
      igwk(3)=0
      igwk(4)=1
      igwk(5)=10
      ligw=20
!
      call dgmres(n,b,x,nelt,ia,ja,a,isym,matvec,msolve,itol,tol,itmax,
     &  iter,err,ierr,iunit,sb,sx,rgwk,lrgw,igwk,ligw,rwork,iwork)
!
      return