[Scilab] Test d'égalité

[invitee6353f29]

Bonjour,

Je travaille actuellement sous Scilab (auquel je ne suis pas habitué) sur un algorithme qui utilise, comme beaucoup d'autres, un test d'égalité dans un "if".
Je vous joint la partie qui pose problème:

dT = 0.001; 
T = 3; 
n = 0; 
k = 0; 

for t = 0:dT:T
        
        if(k == n*dT) then
            n = n + 1;              
        end;
 
     printf('=> k: %f; n: %d\n',k,n)     
   
     k = k + dT;
        
end;

Le test d'égalité se fait parfaitement jusqu'à ce que n prenne la valeur 10, pas au delà. Par contre si je remplace k par t cela marche. Or j'ai besoin de k puisque je dois le réinitialiser à 0 de temps à autre.
De puis si je fais le test "à la main":

        
        if(k = 10*dT) then
            n = n + 1;              
         end;

il ne se passe quand même rien, le test n'est pas validé alors qu'un simple printf me permet de vérifier visuellement l'égalité.

Auriez vous une idée de la source du problème ?

Merci d'avance !
-----

[Jack]

if(k = 10*dT) then

peut-être une erreur involontaire, mais il faut remplacer le "=" par "==".

if(k == n*dT) then

Un test d'égalité sur des flottants est toujours un peu risqué vu que toutes les valeurs ne sont pas codables. Tu as essayé avec un encadrement de la valeur à tester?

[invitee6353f29]

Merci pour ta réponse.

En effet c'est une erreur involontaire lors de la copie sur le forum, j'ai du l'effacer par mégarde.
Tu entends quoi par un encadrement de la valeur à tester ? Elle est déjà encadrée à cause de la boucle, ma variable k est comprise entre 0 et T, comme t d'ailleurs (et comme je l'ai dit cela marche avec t).
J'ai réussi à faire marcher le code avec le test "round(k/dT) == n" mais cela ne me dit quand même pas d'où vient le problème initial. :s

[Jack]

Essaie-voir avec
if((k >= (n*dT) - p) || (k <= (n*dT) + p)) then

|| étant le OU logique et p la précision attendue

[A voir en vidéo sur Futura]

[invitee6353f29]

Ah, j'avais mal compris ce que tu voulais me dire, j'avais déjà essayé et ça marche.
Cependant comme c'est du "bricolage" je ne trouvais pas ça approprié et je voulais faire au plus simple.
Comme apparemment c'est un problème de logiciel plus que de code je vais garder ça.
Merci pour ton aide.