Affichage d'un double

Hello
je réalise l'affichage d'une valeur de type double dans une boucle a
interval régulier, calculée a partir de nombreuses autres variables de type float/double et int.
Au début l'affichage se réalise correctement, mais
au bout d'un moment, la valeur affichée devient:

82.5
83.58
-1.#IND
-1.#IND

Quelqu'un connait-il la signification de cet affichage et à la suite de quelle type d'erreur il apparait?

Citation:

---

Envoyé par vinzzzz

je réalise l'affichage d'une valeur de type double dans une boucle a
82.5
83.58
-1.#IND
-1.#IND

Quelqu'un connait-il la signification de cet affichage et à la suite de quelle type d'erreur il apparait?

---

La valeur ne représente pas un flottant valide (Not a Number ou NaN). IND est peut être une francisation (INDéfini ou INDéterminé...)

Ca sent le débordement de tableau ou la valeur non initialisé, ou le pointeiur fou...

Bref, sans voir le code, ça va être dur d'en dire plus...

La fonction qui calcule la valeur:

Code:

---

1 2 3 4 5 6 7

double calcul_maj(s_param *params, double mass, float subst_c, double uptk_rate) { double n = mass* (params->growth_rt / params->max_rdt) * subst_c ; double d = (params->sat_cst) + subst_c ; double result = uptk_rate*n/d ; return result ; }

---

Les paramètres de la structure params sont tous de type float. Pas de tableau ni de pointeur sauf pour la structure params qui contient des valeurs correctement initialisées. J'ai vraiment aucune idée de ce qui peut se passer...

J'ai affiché l'intégralité des valeurs utilisées pour les calculs et voici ce que j'ai obtenu (affichage console):

Code:

---

1 2 3 4

In calcul_maj: (((1176.96*(0.0205/0.444)*1)) / (0.00234+1))*1 = 54.2147 (...) In calcul_maj: (((1178.13*(0.0205/0.444)*1)) / (0.00234+1))*1 = -1.#IND

---

Ca m'inspire pas beaucoup... (au passage, mon printf utilise %g pour afficher les valeurs)

Avec d'autres valeurs...:

Code:

---

1 2 3 4

In calcul_maj: (((131.761*(0.0205/0.444)*0.001)) / (0.00234+0.001))*1 = 1.82143 In calcul_maj: (((131.761*(0.0205/0.444)*0.001)) / (0.00234+0.001))*1 = 1.82143 In calcul_maj: (((131.798*(0.0205/0.444)*0.001)) / (0.00234+0.001))*1 = -1.#IND

---

...

Citation:

---

Envoyé par vinzzzz

Code:

---

1 2 3 4

<...> double result = uptk_rate*n/d ; <...>

---

---

Ah oui, j'avais oublié çà. Je vois des divisions, pas de divisions par 0 ?

d n'est jamais egal à 0 ? C'est testé ? Validé ?

d est toujours positif. Dans:

Code:

---

1 2	d = (params->sat_cst) + subst_c

---

subst_c est toujours strictement positif c'est vérifié.

Code:

---

1 2 3 4 5 6 7 8

Console:   In calcul_maj: d = 1.00234 (((1175.79*(0.0205/0.444)*1)) / (0.00234+1))*1 = 54.1609   In calcul_maj: d = 1.00234 (((1178.13*(0.0205/0.444)*1)) / (0.00234+1))*1 = -1.#IND

---

8O

J'ai encore mieux, en utilisant pour mass une valeur constante, au bout d'un certain nombre d'itération j'obtient:

Code:

---

1 2 3 4 5 6 7 8 9 10

Console:   In calcul_maj: d = 1.00234 (((1568.07*(0.0205/0.444)*1)) / (0.00234+1))*1 = 72.2304   (...)   In calcul_maj: d = 1.00234 (((1568.07*(0.0205/0.444)*1)) / (0.00234+1))*1 = -1.#IND

---

Tu as une autre division dans ton code:

Code:

---

1 2	double n = mass* (params->growth_rt / params->max_rdt) * subst_c ;

---

Si tu es sûr de ton d, alors la division par 0 viendra sûrement de là...

Jc

En fait toutes mes valeures sont strictement positives (on peut d'ailleur le voir dans les affichages du dessu), et dans ma fonction seul le paramètre "mass" est différent a chaque appel, c'est pourquoi j'ai testé l'appel de la fonction en envoyant une constante numérique dirèctement, sans succès (voir au dessu...).
C'est vraiment bizard.

Citation:

---

Envoyé par vinzzzz

En fait toutes mes valeures sont strictement positives (on peut d'ailleur le voir dans les affichages du dessu), et dans ma fonction seul le paramètre "mass" est différent a chaque appel, c'est pourquoi j'ai testé l'appel de la fonction en envoyant une constante numérique dirèctement, sans succès (voir au dessu...).
C'est vraiment bizard.

---

Tu peux poster un programme complet et minimal permettant de reproduire le problème?

Je vais essayer mais ca va pas être facile jpense que j'aurai pas le temps avant demain.
RDV demain donc.