[C] effet plasma C/SDL - Version imprimable

Code DIFF :
--- a/plasma.c  2013-09-13 20:23:06.017363199 +0200

+++ b/plasma.c  2013-09-13 20:23:12.537281232 +0200

@@ -13,6 +13,7 @@

 #define HAUTEUR        500

 #define PI 3.14159265359

 #define STEP 1000

+#define COEF 1000

 

 double distance (unsigned int x1, unsigned int y1, unsigned int x2, unsigned int y2) {

        // trigonometrie de ~15/16 ans : calcul de l'hypothenuse

@@ -29,6 +30,11 @@

        double f1, f2, f3, deg = PI/180;

        double *_sin = calloc(360*STEP,sizeof(double)), *_cos = calloc(360*STEP,sizeof(double));

        if (!_sin || !_cos) return -1; // au cas ou malloc se chie dessus, pas la peine d'aller plus loin

+       unsigned int step = 360 * STEP;

+       double *func1 = calloc(LARGEUR*HAUTEUR*16, sizeof(double));

+       double *func2 = calloc(LARGEUR*HAUTEUR*16, sizeof(double));

+       double *func3 = calloc(LARGEUR*HAUTEUR*16, sizeof(double));

+       if (!func1 || !func2 || !func3) return -2;

 

        // on initialise la sdl et le mode gfx

        SDL_Init(SDL_INIT_VIDEO);

@@ -47,6 +53,14 @@

                        ((unsigned char)(192+63*cos(x*0.1*deg)))

                );

        }

+       // precalcul des images en memoire (optimisation)

+       for (x=0; x<LARGEUR*4;x++) {

+               for (y=0;y<HAUTEUR*4;y++) {

+                       func1[y*LARGEUR*4+x] = sin(distance(x, y, LARGEUR*2, HAUTEUR*2) * deg * 0.3);

+                       func2[y*LARGEUR*4+x] = cos(distance(x, y, LARGEUR*2, HAUTEUR*2) * deg * 0.7);

+                       func3[y*LARGEUR*4+x] = cos(distance(x, y, LARGEUR*2, HAUTEUR*2) * deg * 1.1);

+               }

+       }

 

        // la boucle principale

        while (1) {

@@ -60,9 +74,9 @@

                for (y=0;y<HAUTEUR;y++) { // la routine principale, on balaye l'ecran, pour chaque pixel on calcule son index dans la palette

                        for (x=0;x<LARGEUR;x++) {

                                // nos fonctions trigo horribles et consomatrices de cycles

-                               f1 = sin(distance(x, y, LARGEUR + LARGEUR * _sin[(t*2)%(360*STEP)], HAUTEUR + (hhaut * _sin[(t)%(360*STEP)])) * deg * 0.3);

-                               f2 = cos(distance(x, y, LARGEUR + hlarg * _sin[(t*3)%(360*STEP)], HAUTEUR + (HAUTEUR * _cos[(t*4)%(360*STEP)])) * deg * 0.7);

-                               f3 = sin(distance(x, y, LARGEUR + hlarg * _cos[(t)%(360*STEP)], HAUTEUR + (hhaut * _cos[(t*2)%(360*STEP)])) * deg * 1.1);

+                               f1 = func1[((y+(int)(HAUTEUR + hhaut * _sin[t%step]))*(LARGEUR<<2))+(x+(int)(LARGEUR + LARGEUR * _sin[(t<<1)%step]))];

+                               f2 = func2[((y+(int)(HAUTEUR+HAUTEUR*_cos[(t<<2)%step]))*(LARGEUR<<2))+(x+(int)(LARGEUR+hlarg*_sin[(t*3)%step]))];

+                               f3 = func3[((y+(int)(HAUTEUR+hhaut*_cos[(t<<1)%step]))*(LARGEUR<<2))+(x+(int)(LARGEUR+hlarg*_cos[(t)%step]))];

                                // on somme et fait la moyenne des trois fonctions point a point, et on "anime" en ajoutant t, le tout donne l'indice dans la palette de couleurs

                                c = pal[((unsigned char)(512 + 512 * ((f1+f2+f3)/3))+t)%1024];

                                // putpixel

@@ -76,6 +90,7 @@

                t += STEP;

 

                SDL_UpdateRect(screen,0,0,0,0); // on update l'affichage

+               SDL_Delay(20);

        }

 

        // on quitte proprement