GAUSS-JORDAN
econs
Messages postés
4030
Date d'inscription
mardi 13 mai 2003
Statut
Membre
Dernière intervention
23 décembre 2008
-
20 avril 2006 à 13:58
cs_leFeu Messages postés 22 Date d'inscription samedi 25 mars 2006 Statut Membre Dernière intervention 18 février 2009 - 17 févr. 2009 à 14:28
cs_leFeu Messages postés 22 Date d'inscription samedi 25 mars 2006 Statut Membre Dernière intervention 18 février 2009 - 17 févr. 2009 à 14:28
Cette discussion concerne un article du site. Pour la consulter dans son contexte d'origine, cliquez sur le lien ci-dessous.
https://codes-sources.commentcamarche.net/source/37098-gauss-jordan
https://codes-sources.commentcamarche.net/source/37098-gauss-jordan
17 févr. 2009 à 14:28
y a il une solution????
merci
30 août 2006 à 16:44
Karim Hasnaoui
16 juil. 2006 à 20:34
comme je l'ai écrit je ne sais plus où, je n'ai pas codé cela, c'est une autre personne qui m'a demandé de le mettre en ligne car il avait un pb avec la notation du code qui ne marchait pas je crois.
Donc vu que je n'ai pas utilisé ce code, je ne peux pas répondre à vos questions.
Désolée et bon courage.
13 juil. 2006 à 18:04
k+c+q*dim >= 32*32
Que veux-tu dire par tenseurs d'ordre 2 ?
Ta fonction est bien sensée marcher quelque soit la taille (carrée) de la matrice non?
6 mai 2006 à 15:04
j'ai essayé d'utiliser ce programme pour inverser une matrice 7*7 , le programme me donne un message d'erreur.
le probleme reside dans la valeur de "c", car à un moment donné le [k+c+q*dim] se trouvant dans cette boucle s'accroi et depasse la dimension des vecteurs.
while(abs(a)<0.000000001)
{
c++;
for(int q=0;q<dim;q++)
{
temp[k+q*dim]=temp[k+q*dim]+temp[k+c+q*dim];
mat2[k+q*dim]=mat2[k+q*dim]+mat2[k+c+q*dim];
}
a=temp[k+k*dim];
je vous remercie d'avance pour les eventuels eclaircissements.
20 avril 2006 à 13:58
Petite suggestion pour cette procédure identité, qui peut, à mon avis être améliorée du point de vue de la vitesse d'exécution :
Dans ton exemple, pour une matrice 10x10, tu fais deux boucles imbriquées (çà c'est bien), mais aussi 100 instructions IF.
Avec le code ci-dessous, tu fais 0 instruction IF, remplacées par une boucle de 10 itérations.
J'imagine que pour de grandes matrices, il doit y avoir un certain gain de temps (pas testé).
# void identite(double mat[],int dim)
# {
# for(int i=0;i<dim;i++)
# {
# for(int j=0;j<dim;j++)
# {
# real(mat[i+j*dim])=0.;
# imag(mat[i+j*dim])=0.;
# }
# }
# for(int k=0;k<dim;k++)
# {
# real(mat[k*(dim+1)])=1.;
# imag(mat[k*(dim+1)])=1.;
# }
# }