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