Discussion :

[Ouranis]

Bonjour,
Je réalise un programme pour calculer la trajectoire d'une bille. Jusqu'ici le programme marche parfaitement, mais quand j'intègre les cosinus et les tangentes au calcul de l'ordonnée de la bille, la fenetre ce referme immédiatement et renvoie 3... Y a-t-il quelque chose de spécial a faire pour la trigo?

Voici mes déclarations des variables de calcul et mon code d'affichage du point.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
int continuer = 1, continuerD = 1, sourisX, sourisY, tempsActuel = 0, tempsPrecedent = 0, tempsDepart = 0, X = 12, Y = 472;
 
double distance, gravite = 9.81, vitesse = sqrt(3*9.81), angle = 0, angleRad = 0, Sin2Teta = 0, Teta2 = 0 ,Z = 472, A = 12;
 
               if(event.button.x >= 745 && event.button.y >= 510 && event.button.x <= 793 && event.button.y <= 558)
                {
                    tempsDepart = SDL_GetTicks();
                    while(tempsActuel - tempsDepart < 10000)
                    {
                        tempsActuel = SDL_GetTicks();
                        if (tempsActuel - tempsPrecedent > 100)
                        {
                            Sin2Teta=(distance*gravite)/(vitesse*vitesse);
                            Teta2=asin(Sin2Teta);
                            angleRad=Teta2/2;
                            angle = (180*angleRad)/M_PI;
                            X++;
                            Z =  Z - ( - gravite / (2 * pow(vitesse, 2) * pow(cos(angleRad), 2)) * X * X + tan(angleRad) * X);
                            Y = (int)Z;
 
                            positionBille.x = X;
                            positionBille.y = Y;
                            tempsPrecedent = tempsActuel;
                            SDL_BlitSurface(bille, NULL, ecran, &positionBille);
                            SDL_Flip(ecran);
                        }
                    }
                }

[Sve@r]

Salut
Tu devrais rajouter du printf() pour être sûr que ça vient bien de ces fonctions. Et aussi pour vérifier si le dénominateur de ta division ne vaut pas 0. Et j'espère aussi que tu as inclus <math.h>...

[Neckara]

Pour utiliser les fonctions cosinus, tangente et sinus, on nous conseille de ne pas utiliser ces fonctions.

Mais plutôt de créer un tableau puis de lire dedans cosinus[x] pour cos(x), ce qui est d'après mes sources plus rapide.

Sinon, return 3, c'est un SEGFAULT, tu touches un espace mémoire qui ne t'appartient pas.

Dans le cas d'une division par 0, c'est un autre type d'erreur.

Essaye surtout de regarder au niveau de tes pointeurs

[LittleWhite]

Bonjour,

Mais plutôt de créer un tableau puis de lire dedans cosinus[x] pour cos(x), ce qui est d'après mes sources plus rapide.

, oui, sur Amiga ou les consoles portables qui n'ont pas de gestion de virgule flottante. Mais en général, maintenant, les PCs sont capable de répondre en un temps appréciable aux fonctions cos() / sin() (même s'il est vrai que le précalcul sera toujours plus rapide).

[Sve@r]

Mais plutôt de créer un tableau puis de lire dedans cosinus[x] pour cos(x), ce qui est d'après mes sources plus rapide.

Pas évident car comment affuter ce tableau ? De 1 en 1 ? De 0.1 en 0.1 ?? De 0.01 en 0.01 ???

[Neckara]

C'est à lui de le déterminer en fonction de son utilisation.

Après comme cela je ne peux évidemment pas donner de chiffres.

Vu que pour la gravité, il prend 9,81 je dirais qu'il devrait affuter de 0,001 en 0,001 (comme ça il est sûr d'avoir une précision qui lui convienne mais de 0,01 en 0,01 pourrait marcher aussi)

[anacharsis]

salut !
le cosinus ne peut pas engendrer d'erreur ! par contre la tangente, l'arc sinus oui.

A+

[souviron34]

dans ce débat la question est surtout

Pour utiliser les fonctions cosinus, tangente et sinus, on nous conseille de ne pas utiliser ces fonctions.

Qui te conseille ça et pourquoi ?

Dans 99.999999999% des cas aujourdhui cela n'a pas d'influence...

Franchement j'aimerais bien connaître les raisons et le contexte d'un tel conseil !!!

[anacharsis]

on nous conseille de ne pas utiliser ces fonctions.

les légendes ont la vie dure ...

A+

[dba01]

J'avais fait un truc comme ça en Basic il y a plus de 20 ans, car je n'avais pas les fonctions trigo dans ma machine.
J'ai retrouvé un vieux texte dans mes archives. Mais je ne sais plus si c'était un brouillon ou un truc en état de marche. D'ailleurs, j'ai deux définitions du COS, et il me semble me rappeler que l'ATN fonctionnait une fois sur deux.

Donc c'est à revérifier et réécrire.


Pi.......................................3.1415926535897932384626433
COS(radians)........................SIN(radians + (Pi / 2) )
TAN(radians)........................SIN(radians) / COS(radians)
COT(radians)........................1 / TAN(radians)


____



function FACTOR(INIFACT) // Factorielle(nombre)

RESUFACT := INIFACT
for BASFACT = (INIFACT-1) to 1 step -1
RESUFACT *= BASFACT
next

return RESUFACT

function SIN(INISINUS) // Sinus(radians)

IN2SINUS := INISINUS % (2 * Pi)
RESUSINUS := IN2SINUS
RESUSINUS -= ( IN2SINUS **3 ) / FACTOR(3 )
RESUSINUS += ( IN2SINUS **5 ) / FACTOR(5 )
RESUSINUS -= ( IN2SINUS **7 ) / FACTOR(7 )
RESUSINUS += ( IN2SINUS **9 ) / FACTOR(9 )
RESUSINUS -= ( IN2SINUS **11) / FACTOR(11)
RESUSINUS += ( IN2SINUS **13) / FACTOR(13)
RESUSINUS -= ( IN2SINUS **15) / FACTOR(15)
RESUSINUS += ( IN2SINUS **17) / FACTOR(17)
RESUSINUS -= ( IN2SINUS **19) / FACTOR(19)

return RESUSINUS

function COS(INICOSIN)

IN2SINUS := INICOSIN % (2 * Pi)
RESUSINUS := 1
RESUSINUS -= ( IN2SINUS **2 ) / FACTOR(2 )
RESUSINUS += ( IN2SINUS **4 ) / FACTOR(4 )
RESUSINUS -= ( IN2SINUS **6 ) / FACTOR(6 )
RESUSINUS += ( IN2SINUS **8 ) / FACTOR(8 )
RESUSINUS -= ( IN2SINUS **10) / FACTOR(10)
RESUSINUS += ( IN2SINUS **12) / FACTOR(12)
RESUSINUS -= ( IN2SINUS **14) / FACTOR(14)
RESUSINUS += ( IN2SINUS **16) / FACTOR(16)
RESUSINUS -= ( IN2SINUS **18) / FACTOR(18)

return RESUSINUS

function ATN(INIARCTN) // ArcTAN(tang) => rad

if INIARCTN < 10 .and. INIARCTN > -10
RESUARCTN := Pi * 1.5
TPNARCTN := RESUARCTN / 20
do while abs(TAN(RESUARCTN) - INIARCTN) > 0.00001
do while TAN(RESUARCTN) > INIARCTN
RESUARCTN -= TPNARCTN
enddo
do while TAN(RESUARCTN) < INIARCTN
RESUARCTN += TPNARCTN
enddo
TPNARCTN *= 0.1
enddo
elseif INIARCTN < -10000
RESUARCTN := (2 * Pi * 90.02 / 360)
elseif INIARCTN < -1000
RESUARCTN := (2 * Pi * 90.3 / 360)
elseif INIARCTN < -100
RESUARCTN := (2 * Pi * 92 / 360)
elseif INIARCTN < 100
RESUARCTN := (2 * Pi * 88 / 360)
elseif INIARCTN < 1000
RESUARCTN := (2 * Pi * 89.7 / 360)
else
RESUARCTN := (2 * Pi * 89.98 / 360)
endif
// if RESUARCTN > Pi ; RESUARCTN -= Pi ; endif

return (RESUARCTN - Pi)


PS : et en plus je n'ai jamais été un as en maths !