1: ! 2: ! 3: ! Fortran kernel for the WAXPY() vector routine 4: ! 5: #include finclude/petscdef.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: do 10,i=1,n 16: w(i) = a*x(i) + y(i) 17: 10 continue 19: return 20: end