File: arrdaxpyf.f

package info (click to toggle)
blitz%2B%2B 1%3A0.10-3.2
  • links: PTS, VCS
  • area: main
  • in suites: jessie, jessie-kfreebsd
  • size: 13,276 kB
  • ctags: 12,037
  • sloc: cpp: 70,465; sh: 11,116; fortran: 1,510; python: 1,246; f90: 852; makefile: 701
file content (26 lines) | stat: -rw-r--r-- 546 bytes parent folder | download | duplicates (5) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
 
       subroutine arrdaxpyf(A, B, N, c1)
       integer N, iters
       double precision A(N,N,N), B(N,N,N)
       double precision c1
       double precision c2
       c2 = - c1

         do k=2,N-1
           do j=2,N-1
             do i=2,N-1
               A(i,j,k) = A(i,j,k) + c1 * B(i,j,k);
             enddo
           enddo
         enddo

         do k=2,N-1
           do j=2,N-1
             do i=2,N-1
               A(i,j,k) = A(i,j,k) + c2 * B(i,j,k);
             enddo
           enddo
         enddo

       return
       end