1: ! 2: ! 3: ! Fortran kernel for the WAXPY() vector routine 4: ! 5: #include <petsc/finclude/petscsys.h> 6: ! 7: subroutine FortranWAXPY(n,a,x,y,w) 8: implicit none 9: PetscScalar a 10: PetscScalar x(*),y(*),w(*) 11: PetscInt n 13: PetscInt i 15: PETSC_AssertAlignx(16,x(1)) 16: PETSC_AssertAlignx(16,y(1)) 17: PETSC_AssertAlignx(16,w(1)) 19: do 10,i=1,n 20: w(i) = a*x(i) + y(i) 21: 10 continue 23: return 24: end