#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "complex.h"

#define FS 32.78 *pow(10,6)


void gen_sig(float *sig, int size);

void g_sig(struct complex *data, struct complex32 *data32, struct complex16 *data16,int size);

int main(){
	const int size = 307200;
	printf("%d\n", size);
	float sig[size];
    struct complex *data;
    struct complex32 *data32;
    struct complex16 *data16;





    //data   = malloc(sizeof(struct complex));
   // data32 = malloc(sizeof(struct complex32));
    //data16 = malloc(sizeof(struct complex16));
 
	printf("FS: %f\n" ,FS);
    if ((data=malloc(sizeof(struct complex)*(size_t)size))==NULL)
    {
      fprintf(stderr, "Out of memory!\n");
      exit(1);
    }

  if ((data16=malloc(sizeof(struct complex16)*(size_t)size))==NULL)
    {
      fprintf(stderr, "Out of memory!\n");
      exit(1);
    }

  if ((data32=malloc(sizeof(struct complex32)*(size_t)size))==NULL)
    {
      fprintf(stderr, "Out of memory!\n");
      exit(1);
    }


    g_sig(data,data32, data16, size);
    
return 0;
}

void g_sig(struct complex *data, struct complex32 *data32, struct complex16 *data16,int size){

	int write = 0;
	int j;
    // Generating a table got 100 MHz channel
 	int cp0 = 160; // First cyclic prefix
	int cp  = 144; // Normal cyclic prefix
    int cpe = 512; // Extended cyclic prefix
    int deltaF = 200; //Hz
    int ofdm_size = 2048;
	for (j=0;j<size;j++){
		data[j].r = cos(20*pow(10,6)*j/FS);//Is is Mhz on this line;
 		data[j].i = sin(20*pow(10,6)*j/FS);
    }
    //printf("%d",(1<<15)); 
    //Truncate data to fixed-point representation
    // I've used 16 and 32 bits representation, 
	for (j=0; j<size;j++){
		data16[j].r = (short)data[j].r*(1<<15)-1;
		data16[j].i = (short)data[j].i*(1<<15)-1;
		data32[j].r = (int)data[j].r*(1<<15)-1;
		data32[j].i = (int)data[j].i*(1<<15)-1;
        printf("%f\n %f\n", data16[j].r, data16[j].i);
						
	}
}





void gen_sig(float *sig, int size){
	int write = 0;
    // Generating a table got 100 MHz channel
 	int cp0 = 160; // First cyclic prefix
	int cp  = 144; // Normal cyclic prefix
    int cpe = 512; // Extended cyclic prefix
    int deltaF = 200; //Hz
    int ofdm_size = 2048;
    
    if (write){ //write to file.txt
    	FILE *f;
		f = fopen("file.txt", "w"); 
		if(f==NULL){
    		exit(1); 
    	}
		fprintf(f, "x[%d]={ ", size);
	      	
		for (int i=1;i<size; i+=2){
			sig[i-1] = sin(20*i/FS);
        	sig[i]   = cos(20*i/FS);
			fprintf(f, "%f, ", sig[i-1]);
		}
    	fprintf(f,"}; \n");
	}   
	for (int i=1;i<size; i+=2){
		sig[i-1] = sin(20*i/FS);
        sig[i]   = cos(20*i/FS);
	}
}