rearrange.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	rearrange (au, irow, icol, ndim, neq)

Function/Subroutine Documentation

rearrange()

subroutine rearrange	(	real*8, dimension(*)	au,
		integer, dimension(*)	irow,
		integer, dimension(*)	icol,
		integer	ndim,
		integer	neq
	)

!
!     modifies the sparse storage mode for the iterative solver of
!     Ernst Rank (pcgsolver)
!
      implicit none
!
      integer irow(*),icol(*),ndim,i,neq,k,icr,istart,idiag,kflag
      real*8 au(*)
!
      kflag=2
!
      call isortiid(irow,icol,au,ndim,kflag)
!
      istart=1
      k=irow(1)
      icr=0
      idiag=0
!
      do i=1,ndim
         if(irow(i).eq.k) then
            icr=icr+1
            cycle
         else
            call isortid(icol(istart),au(istart),icr,kflag)
            icr=1
            istart=i
            k=irow(i)
            idiag=idiag+1
            irow(idiag)=i-1
         endif
      enddo
!
!     last row
!
      call isortid(icol(istart),au(istart),icr,kflag)
      idiag=idiag+1
      irow(idiag)=ndim
!
      return