Actual source code: fwaxpy.F

  1: !
  2: !  "$Id: fmaxpy.F,v 1.10 2001/08/07 03:05:24 balay Exp $";
  3: !
  4: !    Fortran kernel for the WAXPY() vector routine
  5: !
 6:  #include include/finclude/petscdef.h
  7: !
  8:       subroutine FortranWAXPY(n,a,x,y,w)
  9:       implicit none
 10:       PetscScalar  a
 11:       PetscScalar  x(*),y(*),w(*)
 12:       integer n

 14:       integer i

 16:       do 10,i=1,n
 17:         w(i) = a*x(i) + y(i)
 18:  10   continue

 20:       return
 21:       end