#include <manip.h>
#include "compress.h"


#define lg	72

void exec_fen(int fen,int item,double valeur,...);

long *entree,*comp,*rec;


main()
{
nouveauD(1,0,"compress",exec_fen);
return(1);
}


//--------------------  fonction  exec  de  la  fenetre : copie  -------------------------------


void exec_fen(int fen,int item,double valeur,...) 
{
int i;
long t;
double x,y;
double bruit=0;

switch(item)
    {
    case ouverture               :


random(0.2);
entree=malloc(lg*4);
comp=malloc(lg*4);
rec=malloc(lg*4);

for(i=0;i<lg;i++) 
	{
	bruit=bruit*0.8 + 0.2*random(0);
	entree[i]=200000.*cos(i*0.005) + 100000.*bruit;
//	entree[i]=200000.*cos(i*0.005) + 100000.*bruit + (i%2)*80000;
	if(i>40)	entree[i]+=500000.*exp(- 0.05*(i-40));
//	if(i>40)	entree[i]+=200000.*exp(- 0.05*(i-40));
	entree[i]/=16;	//  entre 0 et 50 000  soit 16 bits
//	entree[i]*=2;	// de 0 a 15 e6  soit 21 bits
	}


for(i=0;i<lg;i++) trace(1,i,entree[i]);
graph->couleur_trace=bleu;
style(1);
retrace(1);

alerte(1,"");

printf("comprime 4.1 avec lg=%d : ",lg);
t=TickCount();
for(i=0;i<1000;i++)	  compress_4_1((unsigned long*)entree,(unsigned long*)comp,lg,1);
			decompress_4_1(comp,rec,lg);
y=0;for(i=0;i<lg;i++) {x=entree[i]-rec[i];y+=x*x;} y=sqrt(y/(double)lg);
printf("fin de comprime 4_1  : t=%d   ecart total = %g  \n\n",TickCount()-t,y);


/*
printf("comprime : ");
t=TickCount();
for(i=0;i<1000;i++)	  compress_4_2((unsigned long*)entree,(unsigned long*)comp,lg,1);
			decompress_4_2(comp,rec,lg);
y=0;for(i=0;i<lg;i++) {x=entree[i]-rec[i];y+=x*x;} y=sqrt(y/(double)lg);
printf("fin de comprime 4_2  : t=%d   ecart total = %g  \n\n",TickCount()-t,y);

*/
/*
printf("comprime : ");
t=TickCount();
for(i=0;i<1000;i++)	  compress_7_1((unsigned long*)entree,(unsigned long*)comp,lg,1);
			decompress_7_1(comp,rec,lg);
y=0;for(i=0;i<lg;i++) {x=entree[i]-rec[i];y+=x*x;} y=sqrt(y/(double)lg);
printf("fin de comprime 7_1  : t=%d   ecart total = %g  \n\n",TickCount()-t,y);



printf("comprime : ");
t=TickCount();
for(i=0;i<1;i++)	  compress_7_2((unsigned long*)entree,(unsigned long*)comp,lg,1);
			decompress_7_2(comp,rec,lg);
y=0;for(i=0;i<lg;i++) {x=entree[i]-rec[i];y+=x*x;} y=sqrt(y/(double)lg);
printf("fin de comprime 7_2  : t=%d   ecart total = %g  \n\n",TickCount()-t,y);
*/

for(i=0;i<lg;i++) trace(1,i,rec[i]);
graph->couleur_trace=vert;
style(1);
for(i=0;i<lg;i++) trace(1,i,entree[i]-rec[i]);


                 break;
    case fermeture               :
                 break;
    default  :   break;
    }
}