Actual source code: fwaxpy_bgl.F

  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