matvec.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	matvec (n, x, y, nelt, ia, ja, a, isym)

Function/Subroutine Documentation

matvec()

subroutine matvec	(	integer	n,
		real*8, dimension(*)	x,
		real*8, dimension(*)	y,
		integer	nelt,
		integer, dimension(*)	ia,
		integer, dimension(*)	ja,
		real*8, dimension(*)	a,
		integer	isym
	)

      use omp_lib
!
      implicit none
!
      integer ia(*),ja(*),i,j,n,nelt,isym,nflnei
      real*8 y(*),x(*),a(*)
!
!     number of off-diagonal terms
!
      nflnei=nelt-n
!
c$omp parallel default(none)
c$omp& shared(n,x,a,y,ja,ia,nflnei)
c$omp& private(i,j)
c$omp do
      do i=1,n
         y(i)=a(nflnei+i)*x(i)
         do j=ja(i)+1,ja(i+1)
c            if(ia(j).ne.0) then
               y(i)=y(i)+a(j)*x(ia(j))
c            endif
         enddo
      enddo
c$omp end do
c$omp end parallel
!
      return