daxpyf90.f90

来自「A C++ class library for scientific compu」· F90 代码 · 共 13 行

      SUBROUTINE f90daxpy(a, x, y, n, iters)
      INTEGER n, iters
      DOUBLE PRECISION, DIMENSION (n) :: x, y
      DOUBLE PRECISION a, b

      b = - a

      DO i=1,iters
        y = y + a * x
        y = y + b * x
      END DO
      END SUBROUTINE