[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / d / s512.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1996/04/01 15:02:23  mclareni
* Mathlib gen
*
*
#include "gen/pilot.h"
      SUBROUTINE S512(A,B)
      REAL A(256),B(256),C(256)
      COMMON /TAB/D(126),E(127)
      CALL D701BD
      DO 10 I=3,127,2
      A1 = A(I-1)-A(I+1)
      A2 = A(257-I)-A(259-I)
      C1R = A1+A2
      C2R = A1-A2
      C1I = A(I)-A(258-I)
      C2I = A(I)+A(258-I)
      C1 = D(I-2)
      C2 = D(I-1)
      C3R = C1*C2R-C2*C2I
      C3I = C1*C2I+C2*C2R
      C(I) = C1R+C3R
      C(I+1) = C1I+C3I
      C(258-I) = C1R-C3R
   10 C(259-I) = C3I-C1I
      C(1) = 2.*(A(256)-A(2))
      C(2) = -2.*(A(256)+A(2))
      C(129) = 2.*(A(128)-A(130))
      C(130) = 2.*A(129)
      CALL I128(C,C(2),2,C,C(2),2)
      B(129) = 0.5*C(129)
      DO 20 I=2,128
      E1 = 0.5*(C(I)-C(258-I))
      E2 = E(I-1)*(C(I)+C(258-I))
      B(I) = E2-E1
   20 B(258-I) = E2+E1
      RETURN
      END
Commit	Line	Data
fe4da5cc	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 = C1C2R-C2C2I
	24	C3I = C1C2I+C2C2R
	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