IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

 C Discussion :

Problème de passage de tableau 2D


Sujet :

C

  1. #1
    Membre averti
    Inscrit en
    Avril 2005
    Messages
    13
    Détails du profil
    Informations forums :
    Inscription : Avril 2005
    Messages : 13
    Par défaut Problème de passage de tableau 2D
    Bonjour à tous.

    Je rencontre un petit souci dans le passage d'un tableau 2D en paramètres.

    Je pars de ça dans mon main :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    int ** R = malloc(m*sizeof(double*));
    for ( l=0 ; l<m ; l++ ) { R[l]=malloc(m*sizeof(double)); }
    Juste avant d'appeler une fonction qui travaille sur R, un printf("%d\n",R[0][0]); me renvoie 3.

    Ensuite, je fais appel à une de mes fonctions (voici l'appel et le code de la dite fonction) :

    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
    29
    30
    31
    32
    33
    34
    35
    36
    37
    38
    39
    40
    41
    42
    43
    44
    45
    46
    double ** L = Choleski(R,m);
     
    /* Pour alléger, tout le contenu du main n'est pas affiché.
    Ici, je travaille un peu sur le tableau R, j'assigne des valeurs
    à tous ses éléments avant de l'expédier à la fonction Choleski,
    fonction source de mon problème vraisemblablement. */
     
    double ** Choleski(int ** A, int m)
    {
    	int i, j, k;
    	double sum;
     
    	double ** L = malloc(m*sizeof(double*));
    	for ( k=0 ; k<m ; k++ ) { L[k]=malloc(m*sizeof(double)); }
     
    	// Initialization of L (0 above the diagonal of the matrix)
    	for ( j=1 ; j<m ; j++ )
    	{
    		for ( i=0 ; i<j ; i++ ) { L[i][j]=0.0; }
    	}
     
    	// First colum of L
    	L[0][0] = sqrt((double)A[0][0]);
            printf("%e\n",A[0][0]); /* CE PRINTF AFFICHE UN NOMBRE EN 10^-314 */
    	for ( j=1 ; j<m ; j++ )
    	{
    		L[j][0] = ((double)A[0][j])/L[0][0];
    	}
     
    	// Other columns
    	for ( i=1 ; i<m ; i++ )
    	{
    		printf("%d iterations remaining for Choleski-descomposition of RP.\n",m-1-i);
    		sum=0.0;
    		for ( k=0 ; k<i ; k++ ) { sum += pow(L[i][k],2.0); }
    		L[i][i] = sqrt(((double)A[i][i])-sum);
    		for ( j=i+1 ; j<m ; j++ )
    		{
    			sum=0.0;
    			for ( k=0 ; k<i ; k++ ) { sum += L[i][k]*L[j][k]; }
    			L[j][i] = (((double)A[i][j])-sum)/L[i][i];
    		}
    	}
     
    	return L;
    }
    Le printf mis en évidence à la ligne 19 dans le code de la fonction Choleski me donne un nombre de l'ordre de 10^(-314), au lieu de 3... Tous les autres coefficients de ma matrice R, une fois passés en paramètres ont également cette valeur.
    Il s'en suit pas mal de problèmes de NaN dans les divisions du reste de cette fonction.

    Quelqu'un saurait me dire où se trouve mon erreur ?

    Pas de lien vers la FAQ s'il vous plaît, j'avoue qu'elle ne m'a pas trop aidé et je souhaiterais savoir qu'est-ce qui, dans ce code, fait qu'il ne fonctionne pas...

    Merci d'avance.

  2. #2
    Invité(e)
    Invité(e)
    Par défaut
    Bonjour,
    Est-ce voulu de définir R comme pointeur de pointeur sur int et de faire les allocations avec le type double ?

    Pour éviter ce genre de confusion, on peut écrire
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    int ** R = malloc(m * sizeof *R);
    for ( l=0 ; l<m ; l++ ) { 
        R[l] = malloc( m * sizeof *R[l]); 
    }

  3. #3
    Membre averti
    Inscrit en
    Avril 2005
    Messages
    13
    Détails du profil
    Informations forums :
    Inscription : Avril 2005
    Messages : 13
    Par défaut
    Ah, non, ce n'était pas voulu. En manipulant les deux, je me suis un peu embrouillé, c'est corrigé.

    Dans ton petit algo mabu, sizeof *R est équivalent à sizeof(int) ?

    La grosse erreur était dans mon printf aussi : j'ai utilisé un %e mais, comme le nombre étant déclaré en int, le chiffre affiché était aberrant.

  4. #4
    Invité(e)
    Invité(e)
    Par défaut
    Citation Envoyé par jozog47 Voir le message
    Dans ton petit algo mabu, sizeof *R est équivalent à sizeof(int) ?
    sizeof *R va renvoyer la taille nécessaire au stockage d'un élément stocké dans R.
    Si tu changes ce qui est pointé par R, tu n'as pas à modifier ton malloc
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
     
    double *D;
    int *I;
     
    /* sizeof *D <=> sizeof(double) */
    /* sizeof *I <=> sizeof(int) */;
     
    /*allocation d'un tableau de 10 éléments double et 
      d'un tableau de 10 éléments int */
    D = malloc(10*sizeof *D);
    I = malloc(10*sizeof *I);

  5. #5
    Membre averti
    Inscrit en
    Avril 2005
    Messages
    13
    Détails du profil
    Informations forums :
    Inscription : Avril 2005
    Messages : 13
    Par défaut
    Très bonne astuce, il faudra que j'y pense à l'avenir.

Discussions similaires

  1. Problème de passage de tableau de Python vers Fortran 90
    Par P.Math dans le forum Interfaçage autre langage
    Réponses: 0
    Dernier message: 31/03/2010, 11h01
  2. Problème de passage de tableau entre amfphp et php
    Par stefvat dans le forum Flash
    Réponses: 0
    Dernier message: 04/11/2009, 11h02
  3. Problème de passage de tableau Shell --> AWK
    Par LestoK dans le forum Shell et commandes GNU
    Réponses: 5
    Dernier message: 22/07/2009, 16h07
  4. Problème de passage de tableau en paramètre
    Par Tueur_a_gage dans le forum ASP
    Réponses: 9
    Dernier message: 16/11/2006, 20h09
  5. Réponses: 5
    Dernier message: 02/08/2006, 17h05

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo