CodeS-SourceS est hébergé par Frontier.fr

Ce site au démarrage

begin process at 2012 05 27 19:18:36

Accueil Codes Tutoriels Forum Livres Emploi Services Connexion

Trouver un code source :

dans

[ Dernières recherches ]

Accueil >

SIMULATION DE LA DIFFUSION THERMIQUE

Information sur la source

Note :
Aucune note

Catégorie :Graphique Classé sous :diffusion, chaleur, SDL, simulation Niveau :Débutant Date de création :05/08/2009 Date de mise à jour :07/08/2009 11:03:05 Vu / téléchargé :6 556 / 403

Auteur : ncoder

Ecrire un message privé

Site perso

Ce membre participe au partage de revenus publicitaires

Commentaire sur cette source (5)

Ajouter un commentaire et/ou une note

Description

Cliquez pour voir la capture en taille normale

Ce code simule la diffusion de la chaleur dans une plaque métallique (on peut l'appliquer à n'importe quelle autre diffusion d'ailleurs).

Le principe de la simulation repose sur la "discrétisation" du phénomène : la diffusion est régie par l'équation aux dérivées partielles : d(T)/dt = D*Laplacien(T) où T est la température, t est le temps, et D le coefficient de diffusion. Ensuite, on la "discrétise" (c'est à dire on fait un calcul point par point sachant qu'un point de la plaque ne peut influer que sur son plus proche voisin).

Les températures mises en jeu, pour une bonne simplification, sont comprises entre 0°C et 255°C.

La plaque est initialement à 50°C. Les températures des bords de la plaque sont constantes, pour voir l'évolution du phénomène.
On peut choisir, la température des bords de la plaque (en dur dans le code) : j'ai choisi 255 pour les bords à gauche et en haut.
Lorsque vous lancez le programme, vous pouvez ajouter des "points chauds" manuellement n'importe où sur la plaque :
clic gauche => point à 255°C
clic droit => point à 0°C

Si vous les créez pendant la simulation, leur température restera constante tout au long de la simulation.

Pour lancer la simulation, appuyez sur ENTREE.

Source

#include <stdlib.h>
#include <stdio.h>
#include <windows.h>
#include <SDL/SDL.h>
#include <SDL/SDL_getenv.h>
//Constantes de la fenètre
#define LARGEUR 200
#define HAUTEUR 200
//Constantes de calcul de la diffusion
#define DELAI_RAFRAICHISSEMENT 25 // En millisecondes
#define COEFF_DIFFUSION 0.15
#define TYPE float
//Température initiale de toute la plaque
#define TEMPERATURE_INIT 50
//Nombre maximum de points chauds autorisés
#define NB_PTS_CHAUDS_MAX 100
struct pts_chauds
{
SDL_Rect position;
TYPE temperature;
} pts_chauds[NB_PTS_CHAUDS_MAX];
TYPE temperature[2][LARGEUR][HAUTEUR]; //Deux matrices pour le calcul de la température
SDL_Surface *ecran = NULL, *point = NULL; // Les pointeurs qui vont stocker la surface de l'écran et le pixel
int total = 0;
int nb_pts_chauds = 0;
void Initialisation (void); //On y fixe la température de la plaque et des bords à t = 0
void Simulation (void); //Calcul de la diffusion
void setPixel(SDL_Surface *surface, int x, int y, Uint32 pixel); //Pour afficher un pixel
void ActualiserFenetre (void); //Pour afficher les points chauds après création
void AttentePtsChauds (void); //Pour créer des points chauds temporaires avant de lancer la simulation
void RecreerPtsChauds (void); //Pour permettre aux points chauds permanents de conserver leur température
int main(int argc, char *argv[])
{
SDL_Init(SDL_INIT_VIDEO);
putenv("SDL_VIDEO_WINDOW_POS=center"); //Pour centrer la fenêtre
ecran = SDL_SetVideoMode(LARGEUR, HAUTEUR, 8, SDL_HWSURFACE|SDL_DOUBLEBUF); // On tente d'ouvrir une fenêtre
if (ecran == NULL) // Si l'ouverture a échoué, on écrit l'erreur et on arrête
{
fprintf(stderr, "Impossible de charger le mode vidéo : %s\n", SDL_GetError());
exit(EXIT_FAILURE);
}
SDL_WM_SetCaption("Diffusion", NULL); //On change le nom de la fenêtre
SDL_FillRect(ecran, NULL, SDL_MapRGB(ecran->format, TEMPERATURE_INIT, 0, 255-TEMPERATURE_INIT)); //On remplit l'écran
SDL_Flip(ecran); /* Mise à jour de l'écran */
Initialisation(); //On crée la plaque
AttentePtsChauds(); //Ici l'utilisateur crée les points chauds
Simulation(); //La simulation commence !
SDL_Quit(); //On quitte
return EXIT_SUCCESS;
}
// On crée la plaque
void Initialisation (void)
{
int i,j;
for (i=0; i<LARGEUR; i++)
{
for (j=0;j<HAUTEUR; j++)
{
temperature[0][i][j] = TEMPERATURE_INIT;
temperature[1][i][j] = TEMPERATURE_INIT;
}
}
for(j=0;j<HAUTEUR; j++){
temperature[0][0][j]= 255; temperature[1][0][j]= 255; //Ici le coté gauche de la plaque est fixé à 255°C
//temperature[0][HAUTEUR-1][j]= 255; temperature[1][HAUTEUR-1][j]= 255; //Ici le coté droit de la plaque est fixé à 255°C
}
for(i=0;i<LARGEUR; i++){
temperature[0][i][0]= 255; temperature[1][i][0]= 255; //Ici le haut de la plaque est fixé à 255°C
//temperature[0][i][LARGEUR-1]= 255; temperature[1][i][LARGEUR-1]= 255; //Ici le bas de la plaque est fixé à 255°C
}
}
// Si l'utilisateur veut rajouter des points chauds temporaires à la main
void AttentePtsChauds (void)
{
SDL_Event event;
int continuer = 1;
while (continuer)
{
SDL_WaitEvent(&event);
switch(event.type)
{
case SDL_MOUSEBUTTONUP: //Souris
if (event.button.button == SDL_BUTTON_RIGHT && nb_pts_chauds<NB_PTS_CHAUDS_MAX) /* Si on a fait un clic droit, on rajoute un point temporaire à 0°C */
{
//Ajout du point froid temporaire (0°C)
temperature[0][event.button.x][event.button.y] = 0;
temperature[0][event.button.x-1][event.button.y] = 0;
temperature[0][event.button.x+1][event.button.y] = 0;
temperature[0][event.button.x][event.button.y-1] = 0;
temperature[0][event.button.x][event.button.y+1] = 0;
ActualiserFenetre(); //On affiche les points chauds
}
if (event.button.button == SDL_BUTTON_LEFT && nb_pts_chauds<NB_PTS_CHAUDS_MAX) /* Si on a fait un clic gauche, on rajoute un point temporaire à 255°C */
{
//Ajout du point chaud temporaire (255°C)
temperature[0][event.button.x][event.button.y] = 255;
temperature[0][event.button.x-1][event.button.y] = 255;
temperature[0][event.button.x+1][event.button.y] = 255;
temperature[0][event.button.x][event.button.y-1] = 255;
temperature[0][event.button.x][event.button.y+1] = 255;
ActualiserFenetre(); //On affiche les points chauds
}
break;
case SDL_QUIT:
continuer = 0;
break;
case SDL_KEYDOWN: //Clavier
switch(event.key.keysym.sym)
{
case SDLK_RETURN: //Si on appuie sur entrée, la simulation commence
continuer = 0;
}
break;
}
}
}
void Simulation (void)
{
int i = 0, j = 0, alt = 0, tempsActuel = 0, tempsPrecedent = 0;
int continuer = 1;
SDL_Rect position;
SDL_Event event;
// Ajout des pts chauds
RecreerPtsChauds();
// Boucle
while (continuer)
{
SDL_PollEvent(&event);
switch (event.type)
{
case SDL_QUIT:
continuer = 0;
break;
case SDL_MOUSEBUTTONUP: //Souris
if (event.button.button == SDL_BUTTON_RIGHT && nb_pts_chauds<NB_PTS_CHAUDS_MAX) /* Si on a fait un clic droit, on rajoute un point à 0°C */
{
//Ajout du point froid permanent ((0°C)
pts_chauds[nb_pts_chauds].temperature = 0;
pts_chauds[nb_pts_chauds].position.x = event.button.x;
pts_chauds[nb_pts_chauds].position.y = event.button.y;
temperature[0][event.button.x][event.button.y] = 0;
temperature[0][event.button.x-1][event.button.y] = 0;
temperature[0][event.button.x+1][event.button.y] = 0;
temperature[0][event.button.x][event.button.y-1] = 0;
temperature[0][event.button.x][event.button.y+1] = 0;
nb_pts_chauds++;
}
if (event.button.button == SDL_BUTTON_LEFT && nb_pts_chauds<NB_PTS_CHAUDS_MAX) /* Si on a fait un clic gauche, on rajoute un point à 255°C */
{
//Ajout du point chaud permanent (255°C)
pts_chauds[nb_pts_chauds].temperature = 255;
pts_chauds[nb_pts_chauds].position.x = event.button.x;
pts_chauds[nb_pts_chauds].position.y = event.button.y;
temperature[0][event.button.x][event.button.y] = 255;
temperature[0][event.button.x-1][event.button.y] = 255;
temperature[0][event.button.x+1][event.button.y] = 255;
temperature[0][event.button.x][event.button.y-1] = 255;
temperature[0][event.button.x][event.button.y+1] = 255;
nb_pts_chauds++;
}
break;
case SDL_KEYDOWN: //Clavier
switch(event.key.keysym.sym)
{
case SDLK_ESCAPE: //Echap
continuer = 0;
break;
}
break;
}
// Gestion du temps
tempsActuel = SDL_GetTicks();
if (tempsActuel - tempsPrecedent < DELAI_RAFRAICHISSEMENT)
{
SDL_Delay(DELAI_RAFRAICHISSEMENT - (tempsActuel - tempsPrecedent));
}
tempsPrecedent = SDL_GetTicks();
// Dessin de la surface
SDL_LockSurface(ecran); /* On bloque la surface ecran car on va directement modifier ses pixels */
//Calcul des nouvelles températures
for (i = 1 ; i < LARGEUR -1 ; i++)
{
for (j = 1 ; j < HAUTEUR -1; j++)
{
temperature[!alt][i][j] = temperature[alt][i][j] +
COEFF_DIFFUSION*(
temperature[alt][i][j+1]
+ temperature[alt][i+1][j]
+ temperature[alt][i][j-1]
+ temperature[alt][i-1][j]
- 4*temperature[alt][i][j]);
//On affiche la température du point considéré
setPixel(ecran, i, j, SDL_MapRGB(ecran->format, (temperature[!alt][i][j]), 0, 255-(temperature[!alt][i][j])));
}
}
SDL_UnlockSurface(ecran); /* On a fini de travailler sur l'écran, on débloque la surface */
SDL_Flip(ecran); //On actualise
// On recrée les points chauds pour que leur température soit constante
RecreerPtsChauds();
alt = !alt; //On change de matrice pour le calcul
}
}
/* La fonction setPixel permet de dessiner pixel par pixel dans une surface */
void setPixel(SDL_Surface *surface, int x, int y, Uint32 pixel)
{
int bpp = surface->format->BytesPerPixel;
Uint8 *p = (Uint8 *)surface->pixels + y * surface->pitch + x * bpp;
switch (bpp)
{
case 1:
*p = pixel;
break;
case 2:
*(Uint16 *)p = pixel;
break;
case 3:
if (SDL_BYTEORDER == SDL_BIG_ENDIAN)
{
p[0] = (pixel >> 16) & 0xff;
p[1] = (pixel >> 8) & 0xff;
p[2] = pixel & 0xff;
}
else
{
p[0] = pixel & 0xff;
p[1] = (pixel >> 8) & 0xff;
p[2] = (pixel >> 16) & 0xff;
}
break;
case 4:
*(Uint32 *)p = pixel;
break;
}
}
void ActualiserFenetre (void)
{
int i, j;
// Dessin de la surface
SDL_LockSurface(ecran); /* On bloque la surface ecran car on va directement modifier ses pixels */
for (i = 1 ; i < LARGEUR - 1; i++)
{
for (j = 1 ; j < HAUTEUR - 1; j++)
{
setPixel(ecran, i, j, SDL_MapRGB(ecran->format, (temperature[0][i][j]), 0, 255-(temperature[0][i][j])));
}
}
SDL_UnlockSurface(ecran); /* On a fini de travailler sur l'écran, on débloque la surface */
SDL_Flip(ecran);
}
void RecreerPtsChauds (void)
{
int i;
// Ajout des pts chauds
for (i=0; i<nb_pts_chauds; i++)
{
temperature[0][pts_chauds[i].position.x][pts_chauds[i].position.y] = pts_chauds[i].temperature;
temperature[0][pts_chauds[i].position.x-1][pts_chauds[i].position.y] = pts_chauds[i].temperature;
temperature[0][pts_chauds[i].position.x+1][pts_chauds[i].position.y] = pts_chauds[i].temperature;
temperature[0][pts_chauds[i].position.x][pts_chauds[i].position.y-1] = pts_chauds[i].temperature;
temperature[0][pts_chauds[i].position.x][pts_chauds[i].position.y+1] = pts_chauds[i].temperature;
temperature[0][pts_chauds[i].position.x+1][pts_chauds[i].position.y+1] = pts_chauds[i].temperature;
temperature[0][pts_chauds[i].position.x-1][pts_chauds[i].position.y-1] = pts_chauds[i].temperature;
temperature[0][pts_chauds[i].position.x-1][pts_chauds[i].position.y+1] = pts_chauds[i].temperature;
temperature[0][pts_chauds[i].position.x+1][pts_chauds[i].position.y-1] = pts_chauds[i].temperature;
temperature[1][pts_chauds[i].position.x][pts_chauds[i].position.y] = pts_chauds[i].temperature;
temperature[1][pts_chauds[i].position.x-1][pts_chauds[i].position.y] = pts_chauds[i].temperature;
temperature[1][pts_chauds[i].position.x+1][pts_chauds[i].position.y] = pts_chauds[i].temperature;
temperature[1][pts_chauds[i].position.x][pts_chauds[i].position.y-1] = pts_chauds[i].temperature;
temperature[1][pts_chauds[i].position.x][pts_chauds[i].position.y+1] = pts_chauds[i].temperature;
temperature[1][pts_chauds[i].position.x+1][pts_chauds[i].position.y+1] = pts_chauds[i].temperature;
temperature[1][pts_chauds[i].position.x-1][pts_chauds[i].position.y-1] = pts_chauds[i].temperature;
temperature[1][pts_chauds[i].position.x-1][pts_chauds[i].position.y+1] = pts_chauds[i].temperature;
temperature[1][pts_chauds[i].position.x+1][pts_chauds[i].position.y-1] = pts_chauds[i].temperature;
}
}

#include <stdlib.h>
#include <stdio.h>
#include <windows.h>
#include <SDL/SDL.h>
#include <SDL/SDL_getenv.h>

//Constantes de la fenètre
#define LARGEUR 200
#define HAUTEUR 200

//Constantes de calcul de la diffusion
#define DELAI_RAFRAICHISSEMENT      25     // En millisecondes

#define COEFF_DIFFUSION     0.15
#define TYPE    float

//Température initiale de toute la plaque
#define TEMPERATURE_INIT            50

//Nombre maximum de points chauds autorisés
#define NB_PTS_CHAUDS_MAX           100


struct pts_chauds
{
    SDL_Rect position;
    TYPE temperature;
} pts_chauds[NB_PTS_CHAUDS_MAX];


TYPE temperature[2][LARGEUR][HAUTEUR]; //Deux matrices pour le calcul de la température
SDL_Surface *ecran = NULL, *point = NULL; // Les pointeurs qui vont stocker la surface de l'écran et le pixel

int total = 0;
int nb_pts_chauds = 0;


void Initialisation (void);                 //On y fixe la température de la plaque et des bords à t = 0
void Simulation (void);                     //Calcul de la diffusion
void setPixel(SDL_Surface *surface, int x, int y, Uint32 pixel); //Pour afficher un pixel
void ActualiserFenetre (void);              //Pour afficher les points chauds après création
void AttentePtsChauds (void);               //Pour créer des points chauds temporaires avant de lancer la simulation
void RecreerPtsChauds (void);               //Pour permettre aux points chauds permanents de conserver leur température





int main(int argc, char *argv[])
{
    SDL_Init(SDL_INIT_VIDEO);
    putenv("SDL_VIDEO_WINDOW_POS=center"); //Pour centrer la fenêtre
    ecran = SDL_SetVideoMode(LARGEUR, HAUTEUR, 8, SDL_HWSURFACE|SDL_DOUBLEBUF); // On tente d'ouvrir une fenêtre
    if (ecran == NULL) // Si l'ouverture a échoué, on écrit l'erreur et on arrête
    {
        fprintf(stderr, "Impossible de charger le mode vidéo : %s\n", SDL_GetError());
        exit(EXIT_FAILURE);
    }

    SDL_WM_SetCaption("Diffusion", NULL);   //On change le nom de la fenêtre
    SDL_FillRect(ecran, NULL, SDL_MapRGB(ecran->format, TEMPERATURE_INIT, 0, 255-TEMPERATURE_INIT)); //On remplit l'écran
    SDL_Flip(ecran); /* Mise à jour de l'écran */


    Initialisation();       //On crée la plaque

    AttentePtsChauds();     //Ici l'utilisateur crée les points chauds

    Simulation();           //La simulation commence !

    SDL_Quit();             //On quitte

    return EXIT_SUCCESS;
}


// On crée la plaque
void Initialisation (void)
{
    int i,j;
    for (i=0; i<LARGEUR; i++)
    {
        for (j=0;j<HAUTEUR; j++)
        {
            temperature[0][i][j] = TEMPERATURE_INIT;
            temperature[1][i][j] = TEMPERATURE_INIT;
        }
    }

        for(j=0;j<HAUTEUR; j++){
            temperature[0][0][j]= 255; temperature[1][0][j]= 255;                       //Ici le coté gauche de la plaque est fixé à 255°C
            //temperature[0][HAUTEUR-1][j]= 255; temperature[1][HAUTEUR-1][j]= 255;     //Ici le coté droit de la plaque est fixé à 255°C
        }

        for(i=0;i<LARGEUR; i++){
            temperature[0][i][0]= 255; temperature[1][i][0]= 255;                       //Ici le haut de la plaque est fixé à 255°C
            //temperature[0][i][LARGEUR-1]= 255; temperature[1][i][LARGEUR-1]= 255;     //Ici le bas de la plaque est fixé à 255°C
        }
}



// Si l'utilisateur veut rajouter des points chauds temporaires à la main
void AttentePtsChauds (void)
{
    SDL_Event event;
    int continuer = 1;

    while (continuer)
    {
        SDL_WaitEvent(&event);
        switch(event.type)
        {

            case SDL_MOUSEBUTTONUP: //Souris
                if (event.button.button == SDL_BUTTON_RIGHT && nb_pts_chauds<NB_PTS_CHAUDS_MAX) /* Si on a fait un clic droit, on rajoute un point temporaire à 0°C */
                {
                    //Ajout du point froid temporaire (0°C)
                    temperature[0][event.button.x][event.button.y] = 0;
                    temperature[0][event.button.x-1][event.button.y] = 0;
                    temperature[0][event.button.x+1][event.button.y] = 0;
                    temperature[0][event.button.x][event.button.y-1] = 0;
                    temperature[0][event.button.x][event.button.y+1] = 0;

                    ActualiserFenetre();            //On affiche les points chauds
                }
                if (event.button.button == SDL_BUTTON_LEFT && nb_pts_chauds<NB_PTS_CHAUDS_MAX) /* Si on a fait un clic gauche, on rajoute un point temporaire à 255°C */
                {
                    //Ajout du point chaud temporaire (255°C)
                    temperature[0][event.button.x][event.button.y] = 255;
                    temperature[0][event.button.x-1][event.button.y] = 255;
                    temperature[0][event.button.x+1][event.button.y] = 255;
                    temperature[0][event.button.x][event.button.y-1] = 255;
                    temperature[0][event.button.x][event.button.y+1] = 255;

                    ActualiserFenetre();            //On affiche les points chauds
                }
                break;

            case SDL_QUIT:
                continuer = 0;
                break;

            case SDL_KEYDOWN: //Clavier
                switch(event.key.keysym.sym)
                {
                    case SDLK_RETURN:           //Si on appuie sur entrée, la simulation commence
                        continuer = 0;
                }
                break;
        }

    }

}




void Simulation (void)
{
    int i = 0, j = 0, alt = 0, tempsActuel = 0, tempsPrecedent = 0;
    int continuer = 1;
    SDL_Rect position;
    SDL_Event event;

// Ajout des pts chauds
    RecreerPtsChauds();


// Boucle
    while (continuer)
    {
        SDL_PollEvent(&event);
        switch (event.type)
        {
            case SDL_QUIT:
                continuer = 0;
                break;

            case SDL_MOUSEBUTTONUP: //Souris
                if (event.button.button == SDL_BUTTON_RIGHT && nb_pts_chauds<NB_PTS_CHAUDS_MAX) /* Si on a fait un clic droit, on rajoute un point à 0°C */
                {
                    //Ajout du point froid permanent ((0°C)
                    pts_chauds[nb_pts_chauds].temperature = 0;
                    pts_chauds[nb_pts_chauds].position.x = event.button.x;
                    pts_chauds[nb_pts_chauds].position.y = event.button.y;

                    temperature[0][event.button.x][event.button.y] = 0;
                    temperature[0][event.button.x-1][event.button.y] = 0;
                    temperature[0][event.button.x+1][event.button.y] = 0;
                    temperature[0][event.button.x][event.button.y-1] = 0;
                    temperature[0][event.button.x][event.button.y+1] = 0;

                    nb_pts_chauds++;

                }
                if (event.button.button == SDL_BUTTON_LEFT && nb_pts_chauds<NB_PTS_CHAUDS_MAX) /* Si on a fait un clic gauche, on rajoute un point à 255°C */
                {
                    //Ajout du point chaud permanent (255°C)
                    pts_chauds[nb_pts_chauds].temperature = 255;
                    pts_chauds[nb_pts_chauds].position.x = event.button.x;
                    pts_chauds[nb_pts_chauds].position.y = event.button.y;

                    temperature[0][event.button.x][event.button.y] = 255;
                    temperature[0][event.button.x-1][event.button.y] = 255;
                    temperature[0][event.button.x+1][event.button.y] = 255;
                    temperature[0][event.button.x][event.button.y-1] = 255;
                    temperature[0][event.button.x][event.button.y+1] = 255;

                    nb_pts_chauds++;
                }
                break;

            case SDL_KEYDOWN: //Clavier
                switch(event.key.keysym.sym)
                {
                    case SDLK_ESCAPE: //Echap
                        continuer = 0;
                        break;
                }
                break;
        }


// Gestion du temps
        tempsActuel = SDL_GetTicks();
        if (tempsActuel - tempsPrecedent < DELAI_RAFRAICHISSEMENT)
        {
            SDL_Delay(DELAI_RAFRAICHISSEMENT - (tempsActuel - tempsPrecedent));
        }
        tempsPrecedent = SDL_GetTicks();



// Dessin de la surface
        SDL_LockSurface(ecran); /* On bloque la surface ecran car on va directement modifier ses pixels */


//Calcul des nouvelles températures
        for (i = 1 ; i < LARGEUR -1 ; i++)
        {
            for (j = 1 ; j < HAUTEUR -1; j++)
            {
                temperature[!alt][i][j] = temperature[alt][i][j] +
                                          COEFF_DIFFUSION*(
                                              temperature[alt][i][j+1]
                                              + temperature[alt][i+1][j]
                                              + temperature[alt][i][j-1]
                                              + temperature[alt][i-1][j]
                                              - 4*temperature[alt][i][j]);

                //On affiche la température du point considéré
                setPixel(ecran, i, j, SDL_MapRGB(ecran->format, (temperature[!alt][i][j]), 0, 255-(temperature[!alt][i][j])));
            }
        }

        SDL_UnlockSurface(ecran); /* On a fini de travailler sur l'écran, on débloque la surface */
        SDL_Flip(ecran);        //On actualise

        // On recrée les points chauds pour que leur température soit constante
        RecreerPtsChauds();

        alt = !alt;     //On change de matrice pour le calcul
    }


}









/* La fonction setPixel permet de dessiner pixel par pixel dans une surface */
void setPixel(SDL_Surface *surface, int x, int y, Uint32 pixel)
{
    int bpp = surface->format->BytesPerPixel;

    Uint8 *p = (Uint8 *)surface->pixels + y * surface->pitch + x * bpp;

    switch (bpp)
    {
    case 1:
        *p = pixel;
        break;

    case 2:
        *(Uint16 *)p = pixel;
        break;

    case 3:
        if (SDL_BYTEORDER == SDL_BIG_ENDIAN)
        {
            p[0] = (pixel >> 16) & 0xff;
            p[1] = (pixel >> 8) & 0xff;
            p[2] = pixel & 0xff;
        }
        else
        {
            p[0] = pixel & 0xff;
            p[1] = (pixel >> 8) & 0xff;
            p[2] = (pixel >> 16) & 0xff;
        }
        break;

    case 4:
        *(Uint32 *)p = pixel;
        break;
    }
}





void ActualiserFenetre (void)
{
    int i, j;
    // Dessin de la surface

        SDL_LockSurface(ecran); /* On bloque la surface ecran car on va directement modifier ses pixels */

        for (i = 1 ; i < LARGEUR - 1; i++)
        {
            for (j = 1 ; j < HAUTEUR - 1; j++)
            {
                setPixel(ecran, i, j, SDL_MapRGB(ecran->format, (temperature[0][i][j]), 0, 255-(temperature[0][i][j])));
            }
        }
        SDL_UnlockSurface(ecran); /* On a fini de travailler sur l'écran, on débloque la surface */
        SDL_Flip(ecran);
}


void RecreerPtsChauds (void)
{
    int i;
    // Ajout des pts chauds

    for (i=0; i<nb_pts_chauds; i++)
    {
        temperature[0][pts_chauds[i].position.x][pts_chauds[i].position.y] = pts_chauds[i].temperature;
        temperature[0][pts_chauds[i].position.x-1][pts_chauds[i].position.y] = pts_chauds[i].temperature;
        temperature[0][pts_chauds[i].position.x+1][pts_chauds[i].position.y] = pts_chauds[i].temperature;
        temperature[0][pts_chauds[i].position.x][pts_chauds[i].position.y-1] = pts_chauds[i].temperature;
        temperature[0][pts_chauds[i].position.x][pts_chauds[i].position.y+1] = pts_chauds[i].temperature;

        temperature[0][pts_chauds[i].position.x+1][pts_chauds[i].position.y+1] = pts_chauds[i].temperature;
        temperature[0][pts_chauds[i].position.x-1][pts_chauds[i].position.y-1] = pts_chauds[i].temperature;
        temperature[0][pts_chauds[i].position.x-1][pts_chauds[i].position.y+1] = pts_chauds[i].temperature;
        temperature[0][pts_chauds[i].position.x+1][pts_chauds[i].position.y-1] = pts_chauds[i].temperature;

        temperature[1][pts_chauds[i].position.x][pts_chauds[i].position.y] = pts_chauds[i].temperature;
        temperature[1][pts_chauds[i].position.x-1][pts_chauds[i].position.y] = pts_chauds[i].temperature;
        temperature[1][pts_chauds[i].position.x+1][pts_chauds[i].position.y] = pts_chauds[i].temperature;
        temperature[1][pts_chauds[i].position.x][pts_chauds[i].position.y-1] = pts_chauds[i].temperature;
        temperature[1][pts_chauds[i].position.x][pts_chauds[i].position.y+1] = pts_chauds[i].temperature;

        temperature[1][pts_chauds[i].position.x+1][pts_chauds[i].position.y+1] = pts_chauds[i].temperature;
        temperature[1][pts_chauds[i].position.x-1][pts_chauds[i].position.y-1] = pts_chauds[i].temperature;
        temperature[1][pts_chauds[i].position.x-1][pts_chauds[i].position.y+1] = pts_chauds[i].temperature;
        temperature[1][pts_chauds[i].position.x+1][pts_chauds[i].position.y-1] = pts_chauds[i].temperature;
    }
}

Conclusion

Vous pouvez modifier toutes les températures des points, leur position, la température des bords, la taille de la plaque, etc.. à l'intérieur du code.

Pour ce qui est de la création de points chauds, je laisse les deux possibilités :
à t = 0, ce sont des points chauds "temporaires", (non renouvelés)
à t > 0, ce sont des points chauds "permanents"

Le coefficient de diffusion est arbitraire, et constant.

S'il y a le moindre problème, dites le.

Réalisé sous Code::Blocks, avec la SDL.

Fichier Zip

Les Membres Club peuvent télécharger directement un fichier contenu dans le zip sans télécharger le zip en entier !

Diffusion.cbp1 446 octets
Diffusion.depend185 octets
Diffusion.layout243 octets
Diffusion_exe10 752 octets
main.c 13 374 octets
SDL.dll258 048 octets

Télécharger le zip

Historique

07 août 2009 11:03:05 :: Possibilité de créer des points chauds pendant le calcul

Sources du même auteur

[C/API] SUPPRESSION DE FICHIERS EN DOUBLE (NOMS/CONTENUS)

SAVESCREEN

DECOUPEUR DE MP3

INTERFACE CREATOR

RECHERCHE DE ZÉRO D'UNE FONCTION PAR DICHOTOMIE

Toutes les sources de ncoder

Sources de la même categorie

PLANNING D'EQUIPE par grephit

APPLICATION DE DESSIN DE QUELQUES FIGURES par laguchori

HDR EXPOSURE FUSION par mecrosoft

IRC CLIENT MULTISERVEUR EN MFC (TXIRC) par TeniX

ENTETE DU FICHIER BMP (BIPMAP) par k.Lutchi

Toutes les sources de la catégorie Graphique

Sources en rapport avec celle ci

JEUX SERPENT par antho974

PENDU EN SDL par Damsou91

SIMULATION BASIQUE D'EAU par djine

SIMULATION DE VÉHICULE À 4 ROUES AVEC LA SDL par djine

CONSULTATION (SUITE DE SIMULATION DE LA DIFFRACTION) par titip03

Commentaires et avis

Commentaire de spipod le 07/08/2009 08:49:35

Ton code va surement être utile a beaucoup de débutant.

J'aurais personnelement des remarques sur la construction du code, mais cela n'est que de l'ordre du détail, et ce code est destiné aux débutants, donc très bien comme il est.

Par contre, sur la propagation de la chaleur, as-tu remarqué qu'un point chaud à droite ne disperse sa chaleur que vers le centre. les bords ne sont pas affectés ( i = 1 -> LARGEUR - 1).
Je pense que la procedure de calcul de dispersion est un peu trop simplifiée.

Commentaire de ncoder le 07/08/2009 09:05:05

Merci pour la remarque.

Je ne suis pas un spécialiste de la SDL, c'est vrai, je l'utilise comme "support" simple pour mes programmes.

Un point chaud à droite ne disperse sa chaleur que vers le centre car la température des bords est constante, même de ceux qui n'ont pas été considérés initialement comme "chauds" (on suppose un flux infini en bordure de plaque...).
Suivant le bord auquel tu veux enlever cette caractéristique, il faut corriger i=1 en i=0 ou LARGEUR-1 en LARGEUR.

Pour le calcul de diffusion, je suis d'accord c'est qqchose de simplifié.
Que changerais-tu ?

Commentaire de spipod le 07/08/2009 09:31:12

Je pense qu'il ne devrait pas y avoir de constante de température, juste des points chauds ou froids.
Pour moi, je peux me tromper, la dispersion est circulaire à partir d'un point chaud.
Plutôt que parcourir une matrice (dans un sens unique gauche->droite, bas->haut), il faut partir de chaque point chaud, et propager par cercle jusqu'a ce que la propagation n'affecte plus les points suivants.
Facile à dire, je sais.

Pour que le programme prenne un peu plus d'intérêt, recopie le "case SDL_MOUSEBUTTONUP:" de AttentePtsChauds dans Simulation en enlevant la création de points chaud, mais en laissant l'affectation de la température.

Commentaire de ncoder le 07/08/2009 10:55:06

Bien vu pour le "case SDL_MOUSEBUTTONUP:", je n'avais pas pensé à le mettre dans la simulation.

Pour ce qui est de la création de points chauds, je laisse les deux possibilités :
à t = 0, ce sont des points chauds "temporaires", (non renouvelés)
à t > 0, ce sont des points chauds "permanents"

(Car les pts chauds "temporaires" se résorbent vite du fait de la précision du programme : l'½il ne voit pas la différence entre un bleu (67, 0, 188) et un bleu (68, 0, 187) )

Commentaire de arout le 08/06/2010 21:50:19

salut
je trouve que ton code est intéressant et j'aimerai le comprendre .
mais étant débutant dans la matière, j'aime avoir des explications si vous le permettais bien sur

Ajouter un commentaire

Discussions en rapport avec ce code source dans le forum

problème SDL [ par Synhok ] J'ai essayer le tut sur la librairie SDL et quand je compile avec DevC++, le linker me marque :C:\...\BIN\ld.exe: cannot open -lSDLmain: No such file Texte & SDL en mode Video [ par Gaelle ] Bonjour à tous,Je suis en train d'interfacer une application C à l'aide de SDL. Pour ce qui est de la gestion vidéo et évènements souris aucun problèm librairie SDL et coonio.h [ par gloom ] salut peuple svp qq1 opurrait t'il m'aider voila je cherche a afficher une photo sous console DOS pour cela je sait qu'il ¸faute une librairie SDL mai multithreading et opengl [ par hdi ] je doit realiser une simulation avec un affichage 3D (opengl ou open inventor), pour cela je vais realiser deux threads, un qui s'occupe des calcul po multithreading et opengl [ par hdi ] je doit realiser une simulation avec un affichage 3D (opengl ou open inventor), pour cela je vais realiser deux threads, un qui s'occupe des calcul po SDL [ par sabinus ] SabinusBonjour tout le monde.Je recherche la description en français de toutes les instructions et paramétres de la SDL (Windows xp). Lib SDL -> probleme d'images [ par zmc ] bonjour, je debute plus ou moins en SDL, et je n'ai pas réussi a afficher un sprites (animé ou non) sur une image de fond. Selon la méthode que j'empl la librairie SDL pour dev-c++? [ par mioche ] Bonjour,G trouver des sources interessantes concernant la librairie SDL, moins puissante mais beaucoup plus simple que directx.Cependant, je cherche t Lire un wav avec SDL [ par jayce2703 ] Salut tout le monde!Je n'en peux plus : Je n'arrive pas à trouver un exemple SIMPLE de lecture d'un fichier wav sous SDL. Que ce soit sur le site offi Windows Media Player Simulation Unicast [ par ordy ] BonjourJ'aurais besoin de savoir si quelqu'un c'est comment stocker un flux unicast sans lancé le lecteur Windows Media ?Je dois reussir a stocker un

Nos sponsors

Derniers Blogs

IMAGINE CUP 2012, MAKE A SIGN EN FINALE par junarnoalg

KINECT 1.5 IS OUT ! par Vko

LES ACTUALITéS DE LA SEMAINE SUR C2I.FR (14 MAI - 20 MAI) par richardc

REACTIVE EXTENSIONS : CONSOMMER DES SERVICES AVEC RX PARTIE 3, LES PIèGES à éVITER par Groc

SHAREPOINT BLOG SITE, PROBLèME D'ARCHIVES par junarnoalg

Forum

SUPERVISION DE LA BATTERIE par romanticgirl

RE : OPENCV ET INTERFACE GRAPHIQUE VISUAL STUDIO par buno

RE : SAC A DOS par hadjkaddour

OPENCV ET INTERFACE GRAPHIQUE VISUAL STUDIO par ssaif

RE : EFFACER UNE IMAGE AVEC SDL EN C par mouni12imane

Côté IT

Offres d'emploi

Formations Video

Appels d'offres

Realisation d-une version 2 d-un site e-commerce existant sous prestashop-REPUBLICATION
Budget : 5 000€
Recherche Infogerance Serveur Debian Linux
Budget : 1 800€
Creation d-un site vitrine avec MagentO
Budget : 2 000€

Logiciels

sDEVIS-FACTURES vlPRO (8.1.0.3)
SDEVIS-FACTURES VLPRO (8.1.0.3)

sDEVIS-FACTURES vlPRO a été mis au point pour les particuliers, créateurs, entrepreneurs, artisa...

Cliquez pour télécharger sDEVIS-FACTURES vlPRO
974 Application Server (12.2.4.6)
974 APPLICATION SERVER (12.2.4.6)

Développez de puissantes applications dans un environnement de 'cloud computing', clusterisé, séc...

Cliquez pour télécharger 974 Application Server
vPicture (1.4.2.1)
VPICTURE (1.4.2.1)

Avec vPicture, hébergez vos images facilement et rapidement. vPicture est un utilitaire simple, ...

Cliquez pour télécharger vPicture
Easy-Planning (2.2.1.6)
EASY-PLANNING (2.2.1.6)

Easy-Planning permet de créer des plannings sous la représentation de diagrammes et est adapté au...

Cliquez pour télécharger Easy-Planning
COM-BACKUP (2.0)
COM-BACKUP (2.0)

COM-BACKUP est un logiciel de sauvegarde qui permet de planifier les sauvegardes de vos dossiers ...

Cliquez pour télécharger COM-BACKUP