Actual source code: fwaxpy.F90

      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:       end