Discussion :

[Arry_C]

Bon l'histoire c'est que je dois remplire des matrices carrée d'ordre 2 avec des nombres complexes et multiplier ces matrices entre elles.
Mais je sais pas du tout, mais alors pas du tout comment utiliser les complexes en C.
J'ai chercher dans les cours/FAQ j'ai rien trouver ..


Merci

[Thierry Chappuis]

Bon l'histoire c'est que je dois remplire des matrices carrée d'ordre 2 avec des nombres complexes et multiplier ces matrices entre elles.
Mais je sais pas du tout, mais alors pas du tout comment utiliser les complexes en C.
J'ai chercher dans les cours/FAQ j'ai rien trouver ..


Merci

Si ton compilateur support la dernière norme du langage (C99), le fichier d'en-tête définit des macros et fonctions viant à supporter l'arithmétique des nombres complexes. Exemple simple:

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
#include <stdio.h>
#include <stdlib.h>
#include <complex.h>
 
int main(void)
{
    double complex a = 1 + 2*I;
    double complex b = 2 + 2*I;
 
    /* Calcul la somme des nombres complexes a et b */
    double complex z = a + b;
    printf("a + b vaut %.2f%+.2fi\n", creal(z), cimag(z));
 
    /* Calcul le conjugé de a + b */
    z = conj(z);
    printf("le conjugé de a + b vaut %.2f%+.2fi\n", creal(z), cimag(z));
 
    return EXIT_SUCCESS;
}

ATTENTION: Le support des nombres complexes via <complex.h> a le status "broken" dans gcc: http://gcc.gnu.org/c99status.html

Thierry

[Thierry Chappuis]

Une autre solution est de créer ton propre type de nombres complexes:

Code :

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

1
2
3
4
typedef struct MyComplex_ {
    double real;
    double imag;
} MyComplex_s

ansi que les fonctions nécessaires pour manipuler cette structure:

Code :

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

1
2
3
4
5
6
 
void MyComplex_add(MyComplex_s *result, 
                              MyComplex_s const *z1, 
                              MyComplex_s const *z2);
void MyComplex_conj(MyComplex_s *result, 
                              MyComplex_s const *z1);

etc.

Thierry

[Arry_C]

Merci bien

J'ai fais la methode 2 qui à l'avantage de passer avec n'importe quel compileur