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