1c96caaccSSatish Balay! 2c96caaccSSatish Balay! 3c96caaccSSatish Balay! Fortran kernel for sparse matrix-vector product in the AIJ/CRL format 4c96caaccSSatish Balay! 5c96caaccSSatish Balay#include <petsc/finclude/petscsys.h> 6c96caaccSSatish Balay! 7c96caaccSSatish Balaysubroutine FortranMultCRL(m,rmax,x,y,icols,acols) 8c96caaccSSatish Balay implicit none 9c96caaccSSatish Balay PetscInt m,rmax,icols(m,rmax) 10c96caaccSSatish Balay PetscScalar x(0:m-1),y(m) 11c96caaccSSatish Balay PetscScalar acols(m,rmax) 12c96caaccSSatish Balay 13c96caaccSSatish Balay PetscInt i,j 14c96caaccSSatish Balay 15*0113e719SMartin Diehl do j=1,m 16c96caaccSSatish Balay y(j) = acols(j,1)*x(icols(j,1)) 17*0113e719SMartin Diehl end do 18c96caaccSSatish Balay 19*0113e719SMartin Diehl do i=2,rmax 20*0113e719SMartin Diehl do j=1,m 21c96caaccSSatish Balay y(j) = y(j) + acols(j,i)*x(icols(j,i)) 22*0113e719SMartin Diehl end do 23*0113e719SMartin Diehl end do 24c96caaccSSatish Balay 25*0113e719SMartin Diehlend subroutine FortranMultCRL 26