Actual source code: fmultcrl.F90

  1: !
  2: !
  3: !    Fortran kernel for sparse matrix-vector product in the AIJ/CRL format
  4: !
  5: #include <petsc/finclude/petscsys.h>
  6: !
  7: pure subroutine FortranMultCRL(m, rmax, x, y, icols, acols)
  8:   use, intrinsic :: ISO_C_binding
  9:   implicit none(type, external)
 10:   PetscInt, intent(in) :: m, rmax, icols(m, rmax)
 11:   PetscScalar, intent(in) :: x(0:m - 1), acols(m, rmax)
 12:   PetscScalar, intent(out) :: y(m)

 14:   PetscInt :: i

 16:   y(1:m) = acols(1:m, 1)*x(icols(1:m, 1))
 17:   do i = 2, rmax
 18:     y(1:m) = y(1:m) + acols(1:m, i)*x(icols(1:m, i))
 19:   end do
 20: end subroutine FortranMultCRL