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