[C] Besoin d'aide pour débugger mon code

[Polo]

Bonjour à tous,
J'essaie actuellement d'adapter une formule mathématique (permettant de calculer la latitude à partir de la latitude isométrique) en langage C, dont le schéma séquentiel est le suivant :


Donc pour l'instant j'ai réussit à faire ce code, mais j'ai manifestement fait une erreur quelque part car il me retourne toujours le même nombre, à savoir 4.940656e-324 :

Code C :

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <errno.h>


int main()
{
    double L = 1.00552653648, e = 0.08199188998;
    float pi = 3.141592654;
    int i = 0;
    double epsilon = -1 * pow(10, 11);

    double *Phi = malloc(i * sizeof(double));
    if (Phi == NULL)
        exit(2);

    Phi[0] = 2 * atan( exp(L) ) - pi/2;

    do
    {
        i++;

        Phi = realloc(Phi, i * sizeof(double));

        if (Phi == NULL)
            exit(2);

        Phi[i] = log( tan((pi/4) + (L/2)) * pow( ((1-e*sin(L)) / (1+e*sin(L))), e/2 ) );

        printf("%le", &Phi[i]);

    } while(abs(Phi[i]-Phi[i-1]) < epsilon);

    //double Phi_f = Phi[i];

    free(Phi);
}
 

PS : j'ai oublié de vous prévenir que mon niveau était plutôt classé dans "débutant" en C

Toute aide est la bienvenue ... :/

EDIT : J'ai oublié de préciser que j'ai des valeurs de test dont je connais les futurs résultats :

---

Y'a du nouveau ... je viens de réaliser que je m'étai trompé dans la valeur de epsilon (j'avais mis -1 * pow(10, 11) au lieu de pow(10,-11) ) et l'erreur qu'il me renvoie est désormais :

Code :

realloc(): inval id next size: 0x00000000013fc010

:/

[supersnail]

Bonjour,

Déjà

Code :

&Phi[i]

ça te renvoie l'adresse en mémoire de la valeur Phi[i]: donc là tu affiches à l'écran l'adresse de ton pointeur interprété comme un double par printf :')

Un

Code C :

 printf("%le", Phi[i]);

suffirait ici

Enfin (remarque purement "esthétique"), vaut mieux allouer avec malloc un tableau d'une certaine taille, et le réallouer que lorsque ça "déborde" plûtot que de realloc à chaque fois qu'on veut rajouter un truc du temps (realloc ça peut pomper au niveau perfs, donc autant en réduire l'usage )

Edit: erreur du pow corrigée (ou presque) :')

[Polo]

Merci pour toutes ces remarques je vais commencer par corriger tout ça

Concernant mon epsilon j'ai corrigé l'erreur (cf en bas de mon premier post) c'était une erreur d'étourderie, comme bien souvent ^^

Par contre concernant ta remarque "estéthique", qui serait apparemment une solution à l'erreur que j'ai avec realloc (cf en bas du premier post ... ), j'en ai tenu compte et j'ai modifié mon code comme ceci :

Code C :


int main()
{
    double L = 1.00552653648, e = 0.08199188998;
    float pi = 3.141592654;
    int i = 10, n = 0; // n sert de compteur dans la boucle et i définit la taille de base du tableau de Phi
    double x = pow(10, -11);

    double *Phi = malloc(i * sizeof(double));
    if (Phi == NULL)
        exit(2);

    Phi[0] = 2 * atan( exp(L) ) - pi/2;
    printf("%le \n", Phi[0]);
    do
    {
        n++;
        if(sizeof(Phi) >= n) // si ça commence à déborder...
        {
            i += 10; // on en rajoute 10 d'un coup ...
            Phi = realloc(Phi, i * sizeof(double));
            if (Phi == NULL)
                exit(3);
        }


       // on continue ...

        Phi[i] = log( tan((pi/4) + (L/2)) * pow( ((1-e*sin(L)) / (1+e*sin(L))), e/2 ) );

        printf("%le", Phi[i]);

    } while(abs(Phi[i]-Phi[i-1]) < x); // tant que notre marge d'erreur n'est pas atteinte, on continue...

    //double Phi_f = Phi[i];

    free(Phi); // on libère la place en mémoire

    return 0;
}
 

Donc ça marche déjà mieux (façon de parler), mais je n'ai toujours pas les bons résultats ; j'obtiens ceci :

Code :

8.69343e-01
1.230778e+00

Donc là je comprends pas car il est apparement arrivé à la fin de la boucle et ... j'ai pas ce que je devrais avoir (cad : 0.872 664 626 00) oO
Correction : il est bien arrivé à la fin de la boucle mais n'itère qu'une seule fois... Je vais chercher pourquoi (vous pouvez vous aussi hein )...

[notfound]

Je passe vite fait ! J'ai regardé 2 secondes ton truc, t'as FAIL au niveau de :

Code C :


 printf("%le", Phi[i]);
/* J'ai remplacé par %f par exemple */
 

Le résultat m'a l'air cohérent au détail près que :

Code BASH :


┌ NotFound@Debian_squeeze:/tmp$
└→  ./a.out
Résultat : 0.040996
 

Je n'obtiens que la moitié du résultat (0.040996 * 2 = 0,081991)
En espérant t'avoir aidé un minimum

Hackmicalement

[Polo]

Merci de ces conseils

Hier MacYavel m'a apporté beaucoup d'aide sur IRC, et on a finalement réussit à aboutir à ce code (ça peut peut-être intéresser qqun ) :

Code C :


int main ()
{
    double L = 1.00552653648, e = 0.08199188998;
    double pi = 3.141592654;
    double x = 1e-11;
    double a, b = 0;

    b = 2 * atan( exp(L) ) - pi/2;

    do {

    a = b;
    // formule secondaire de calcul de b :
    // b = log( tan((pi/4) + (L/2)) * pow( ((1+e*sin(a)) / (1-e*sin(a))), e/2 ) );
    b = 2*atan( pow( ((1+e*sin(a)) / (1-e*sin(a))), e/2 ) * exp(L)) - pi/2;


    } while(abs(a-b) > x);

    printf("%le", b);
    return b;
}
 

Ce code est normalement fonctionnel.