program Mat_Vec_Mult

!--This does matrix vector multiplication--!

implicit none
REAL(KIND=8):: sum
REAL(KIND=8),allocatable,dimension(:):: b,x
REAL(KIND=8),allocatable,dimension(:,:):: A
INTEGER:: i,j,size

size = 3

ALLOCATE(b(1:size),x(1:size))
ALLOCATE(A(1:size,1:size))


do i=1,size
  do j=1,size
    A(i,j) = (i-1)*(j-2)+1.0;
  end do
  x(i) = 1.0;
end do

WRITE(*,*) A
! Matrix Multiply !
do i=1,size
  sum = 0.0;
  do j=1,size
    sum = sum + A(i,j)*x(j)
  end do
  b(i) = sum;
end do

WRITE(*,*) b

stop
end program Mat_Vec_Mult