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! 70ccf82acSMartin Diehlpure subroutine FortranMultCRL(m, rmax, x, y, icols, acols) 8*fe66ebccSMartin Diehl use, intrinsic :: ISO_C_binding 90ccf82acSMartin Diehl implicit none(type, external) 100ccf82acSMartin Diehl PetscInt, intent(in) :: m, rmax, icols(m, rmax) 110ccf82acSMartin Diehl PetscScalar, intent(in) :: x(0:m - 1), acols(m, rmax) 120ccf82acSMartin Diehl PetscScalar, intent(out) :: y(m) 13c96caaccSSatish Balay 14d66e387eSMartin Diehl PetscInt :: i 15c96caaccSSatish Balay 16d66e387eSMartin Diehl y(1:m) = acols(1:m, 1)*x(icols(1:m, 1)) 170113e719SMartin Diehl do i = 2, rmax 18d66e387eSMartin Diehl y(1:m) = y(1:m) + acols(1:m, i)*x(icols(1:m, i)) 190113e719SMartin Diehl end do 200113e719SMartin Diehlend subroutine FortranMultCRL 21