#include "LDPC_head.h"
#include <stdlib.h>


code *Gaussian_channel(code *codewords,double n0)
{
	unsigned int i;
	double u1,u0;
	double gaussnoise;

	for(i=0;i<codewords->length;i++)
	{
		u0 = (double)rand()/(double)RAND_MAX;
		u1 = (double)rand()/(double)RAND_MAX;
	
		if( u0 < limit) 
		{
			u0 = limit;
		}
		gaussnoise = sqrt( n0*log(1./u0) )*cos( 2*pi*u1 );
//		gaussnoise=0;
		codewords->codeword[i]+=gaussnoise;
	}
	return (codewords);
}