source: Sophya/trunk/Poubelle/archediab.old/prog_compression_unsigned/test_compress.c

#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;
    }
}