1: !
  2: !
  3: !    Fortran kernel for the AYPX() vector routine
  4: !
  5: #include <petsc/finclude/petscsys.h>
  6: !
  7: pure subroutine FortranAYPX(n,a,x,y)
  8:   use, intrinsic :: ISO_C_binding
  9:   implicit none (type, external)
 10:   PetscScalar, intent(in) :: a
 11:   PetscScalar, intent(in) :: x(*)
 12:   PetscScalar, intent(inout) :: y(*)
 13:   PetscInt, intent(in) :: n

 15:   y(1:n) = x(1:n) + a*y(1:n)
 16: end subroutine FortranAYPX