File: cxloop.f

package info (click to toggle)
nastran 0.1.95-2
  • links: PTS, VCS
  • area: non-free
  • in suites: bookworm, bullseye, sid
  • size: 122,540 kB
  • sloc: fortran: 284,409; sh: 771; makefile: 324
file content (14 lines) | stat: -rw-r--r-- 398 bytes parent folder | download | duplicates (2) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
      SUBROUTINE CXLOOP (X,Y,N)
      DOUBLE PRECISION X(1), Y(1)
      DOUBLE PRECISION XX(2) , YY(2), MPY(2)
      NN = N + N
      DO 10 I = 1,NN
   10 X(I) = Y(I)
      RETURN
      ENTRY CLOOP( XX, YY, MPY, M)
      MM = M+M
      DO 20 I = 1,MM,2
      XX(I) = XX(I) - MPY(1)*YY(I) + MPY(2) * YY(I+1)
   20 XX(I+1) = XX(I+1) -MPY(2)*YY(I) -MPY(1)*YY(I+1)
      RETURN
      END