/*
	This subroutine computes the Decimation in Time (DIT)
	Fast Fourier Transform (FFT) of its fixed point input.
	The integer value of the block exponent is returned
	by fft2048.

	The C Applications Library

	Created on 7/20/93

	#include <ffts.h>;
	int fft2048(int in_real[], int out_real[], int out_real[], int out_imag[]);
		NOTE: The input and output arrays must be different & the output
		      arrays must be one circular boundries
*/

.CONST          	N_=2048, N_div_2_=1024;    /*Const. for 2048 points*/
.CONST			nover2_=1024, nover4_=512;
.CONST			mod_value_=H#0008;
.CONST			log2Nminus2_=9, Nminus1_=2047;

.MODULE/IMAGE		__2048_DIT_FFT_Routine;
.VAR/PM/ANY/CIRC        twid_real [N_div_2_];
.VAR/PM/ANY/CIRC        twid_imag [N_div_2_];


#include		<lib_glob.h>;
#include 		<fft_glob.h>;

.EXTERNAL		__fftcore;

#include	<t2048_rl.h>;
#include	<t2048_im.h>;

.ENTRY		fft2048_;

fft2048_:	I1=I4;		/* push values for fft2048.  We haven't
				   done frame stuff, but we'll be okay */
		M3=-2;
		MODIFY(I1,M3);  /* pass over prev frame and return address */
		M3=-1;
		AX1=N_;
		DM(I1,M3)=AX1;
		AX1=mod_value_;
		DM(I1,M3)=AX1;
		AX1=Nminus1_;
		DM(I1,M3)=AX1;
		AX1=nover2_;
		DM(I1,M3)=AX1;
		M3=-3;
		MODIFY(I1,M3);  /* skip over out_* and node_space */
		M3=-1;
		AX1=^twid_real;
		DM(I1,M3)=AX1;
		AX1=^twid_imag;
		DM(I1,M3)=AX1;
		MODIFY(I1,M3);   /* skip over bflys_per_group and groups */
		MODIFY(I1,M3);
		AX1=nover4_;
		DM(I1,M3)=AX1;
		AX1=log2Nminus2_;
		DM(I1,M3)=AX1;
		JUMP __fftcore;

.ENDMOD;