daxpyf90.f

来自「数值计算工具库,C语言编写的,可以直接调用.」· F 代码 · 共 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