File: daxpyf90-2.f90

package info (click to toggle)
blitz%2B%2B 1%3A0.9-10
  • links: PTS, VCS
  • area: main
  • in suites: squeeze
  • size: 12,684 kB
  • ctags: 12,926
  • sloc: cpp: 97,336; sh: 8,422; fortran: 1,208; makefile: 688; f90: 596
file content (15 lines) | stat: -rw-r--r-- 314 bytes parent folder | download | duplicates (10) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
! Fortran 90 DAXPY using arrays
SUBROUTINE f90daxpy(a, x, y, n, iters)
  IMPLICIT NONE
  INTEGER, INTENT( IN ) :: n, iters
  DOUBLE PRECISION, DIMENSION (n) :: x, y
  DOUBLE PRECISION, INTENT( IN ) :: a
  DOUBLE PRECISION :: b

  b = - a

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