MINICERN/mathlib/gen/d/c64.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:02:23  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 #include "gen/pilot.h"
  10       SUBROUTINE C64(A,B)
  11       REAL A(33),B(33),C(32)
  12       COMMON /TAB/D(126),E(127)
  13       CALL D701BD
  14       S2 = A(16)+A(18)
  15       DO 10 I=3,15,2
  16       S2 = S2+A(I-1)+A(35-I)
  17       A1 = A(I+1)-A(I-1)
  18       A2 = A(35-I)-A(33-I)
  19       C1R = A(I)+A(34-I)
  20       C2R = A(I)-A(34-I)
  21       C1I = A1-A2
  22       C2I = A1+A2
  23       C1 = -D(8*I-9)
  24       C2 = D(8*I-8)
  25       C3R = C1*C2R-C2*C2I
  26       C3I = C1*C2I+C2*C2R
  27       C(I) = C1R+C3R
  28       C(I+1) = C1I+C3I
  29       C(34-I) = C1R-C3R
  30    10 C(35-I) = C3I-C1I
  31       C(1) = A(1)+A(33)
  32       C(2) = A(1)-A(33)
  33       C(17) = 2.*A(17)
  34       C(18) = 2.*(A(16)-A(18))
  35       CALL D16(C,C(2),2,C,C(2),2)
  36       S2 = 2.*S2
  37       B(1) = C(1)+S2
  38       B(33) = C(1)-S2
  39       B(17) = C(17)
  40       DO 20 I=2,16
  41       E1 = 0.5*(C(I)+C(34-I))
  42       E2 = E(8*I-8)*(C(I)-C(34-I))
  43       B(I) = E1+E2
  44    20 B(34-I) = E1-E2
  45       RETURN
  46       END