#include <stdio.h>
#include <iostream>
#include <fstream>

#include "GenFading.h"

const int BRANCHNO=13; // number of independent fading processes

void main()
{
	CFadingEnergyGenerator* pGen = new CFadingEnergyGenerator();


	double fd = 100*18.46/1000/4; // 100Hz Doppler frequency for GSM burst interval
	pGen->Initialize(BRANCHNO, fd);
	
	std::ofstream coeffRealFile("coeffReal.dat");
	std::ofstream coeffImagFile("coeffImag.dat");

	double coeffReal[BRANCHNO], coeffImag[BRANCHNO];

//	unsigned long sampleNo = 80000/(18.46/4)*1000;
	unsigned long sampleNo = 5000;
	for(unsigned int ii=0;ii<sampleNo;ii++)
	{
		pGen->Update(1);
		pGen->GetComplexEnvelope(coeffReal, coeffImag);

		for(int jj=0;jj<BRANCHNO;jj++)
		{
			 coeffRealFile<<coeffReal[jj]<<"  ";
			 coeffImagFile<<coeffImag[jj]<<"  ";
		}

		coeffRealFile<<"\n";
		coeffImagFile<<"\n";
	}

	coeffRealFile.close();
	coeffImagFile.close();
}