MINICERN/mathlib/gen/d/s512.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 S512(A,B)
  11       REAL A(256),B(256),C(256)
  12       COMMON /TAB/D(126),E(127)
  13       CALL D701BD
  14       DO 10 I=3,127,2
  15       A1 = A(I-1)-A(I+1)
  16       A2 = A(257-I)-A(259-I)
  17       C1R = A1+A2
  18       C2R = A1-A2
  19       C1I = A(I)-A(258-I)
  20       C2I = A(I)+A(258-I)
  21       C1 = D(I-2)
  22       C2 = D(I-1)
  23       C3R = C1*C2R-C2*C2I
  24       C3I = C1*C2I+C2*C2R
  25       C(I) = C1R+C3R
  26       C(I+1) = C1I+C3I
  27       C(258-I) = C1R-C3R
  28    10 C(259-I) = C3I-C1I
  29       C(1) = 2.*(A(256)-A(2))
  30       C(2) = -2.*(A(256)+A(2))
  31       C(129) = 2.*(A(128)-A(130))
  32       C(130) = 2.*A(129)
  33       CALL I128(C,C(2),2,C,C(2),2)
  34       B(129) = 0.5*C(129)
  35       DO 20 I=2,128
  36       E1 = 0.5*(C(I)-C(258-I))
  37       E2 = E(I-1)*(C(I)+C(258-I))
  38       B(I) = E2-E1
  39    20 B(258-I) = E2+E1
  40       RETURN
  41       END