Discussion :

[le_voisin]

Salut,

j essaye d implementer la structure suivante :

Code :

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

1
2
3
4
5
6
7
 
while(for(int col=0; col<m.size2(); col=col+1)
			{
				for(int lig=col; lig<m.size1(); lig=lig+1)
				{m(lig,col)!=0;}
			})
{ ....}

Mais ca ne compile et ca me sort plein d erreur, alors je pense qu on ne peut pas implementer une telle condition dans un while, est ce que quelqu un peut me confirmer cela ?

merci, a+

[roulious]

Je confirme. La condition du while doit etre une valeur convertible en int, et une boucle for n'a pas de valeur. Que veux-tu faire exactement ?

[le_voisin]

En fait le but est de realiser en decomposition LU d une matrice. Et cette structure permet de gerer un cas ou l algorithme rencontre un probleme lorsqu il y a un 0 au maivais endroit dans la matrice. Dans cas cas l algorithme saute une etape et la decomposition n est alors pas termine.
Avec cette structure while, je verifie que la decomposition est bien terminee, et si elle n est pas termine (dû a ce 0 mal place), je relance l algorithme qui en partant du debut sur la matrice issue de la premeire decomposition avortee devrait aboutir.
Voilà, j espere que c est clair !

Pour ceux que ca interesse je met l algorithme en entier :

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
47
48
49
50
51
52
53
54
 
#include <iostream> 
#include <vector>
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/io.hpp>
using namespace std;
using namespace boost::numeric::ublas;
 
std::vector<matrix<float>> decomp_lu(matrix<float> m)
{
	std::vector<matrix<float>> V_results(m.size1());
 
	matrix<float> l(m.size1(), m.size1());
	cout << "dimension de la matrice soumise au LU :"<<m.size1()<<"\n";
 
	while(for(int col=0; col<m.size2(); col=col+1)
			{
				for(int lig=col; lig<m.size1(); lig=lig+1)
				{m(lig,col)!=0;}
			})
	{
		for(int i=0; i<=m.size1()-2; i=i+1)//calcul de L, l'operation est faite (n-1) fois
		{
			//coef à 0 de la matrice L
			for(int c=0; c<m.size1(); ++ c)
			for(int d=0; d<m.size2(); ++ d)
			l(c,d)=0;
			//coef à 1 sur la diagonale de la matrice L
			for(int c=0; c<m.size1(); ++ c)
			l(c,c)=1;
 
			if (m(i,i)!=0)
			{
			//calcul des coef non nuls de la matrice L
			for(int j=i; j<=m.size1()-2; ++ j)
			l(j+1,i)=-m(j+1,i)/m(i,i);
			}
			if (m(i,i)==0)
			{
				cout<< "SAUT \n";
				goto jump;
			}
 
			V_results[i]=l; //ecriture de la ieme matrice de L
 
			cout << "operation numero "<< i <<", L : "<< l <<"\n";
			m=prod(l, m);
			cout << "operation numero "<< i <<", U : "<< m <<"\n";
			jump:;
		}
	}//fin du while
	V_results[m.size1()-1]=m;
	return(V_results);
}

[roulious]

Donc tu veux quelque chose du genre :

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
 
bool decompositionTerminee = false;
while(!decompositionTerminee)
{
    // Decomposition ...
 
    decompositionTerminee = false;
 
    for(int col=0; col<m.size2(); col=col+1)
    {
        for(int lig=col; lig<m.size1(); lig=lig+1)
        {
            decompositionTerminee &= (m(lig,col)!=0);
        }
    }
}

Pendant que j'y suis :

Code :

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

1
2
 
std::vector<matrix<float>> decomp_lu(matrix<float> m);

ne compilera pas : la chaine >> est un operateur, donc ton compilateur ne comprendra pas la declaration. Il faut mettre un espace entre les deux :

Code :

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

1
2
 
std::vector<matrix<float> > decomp_lu(matrix<float> m);

[le_voisin]

euh ca a l air interessant ce que tu me dis, par contre comme je commence en c++ je connais pas tout les symboles :
while(!decompositionTerminee)

decompositionTerminee &= (m(lig,col)!=0);

si tu peux me dire en francais ce que cela signifie.

Sinon concernant : std::vector<matrix<float>> decomp_lu(matrix<float> m);
ca compile sans probleme.

Merci, a+

[roulious]

while(!decompositionTerminee)

Le ! est l'operateur booleen NON. Donc while(!condition) se lit en francais : tant que ma condition n'est pas remplie.

decompositionTerminee &= (m(lig,col)!=0);

Le & est l'operation booleenne ET. Les deux boucles for imbriquees verifient que tous les elements de m sont non nuls. Si l'un d'eux est nul, decompositionTerminee passe a false, et y reste.

Sinon concernant : std::vector<matrix<float>> decomp_lu(matrix<float> m);
ca compile sans probleme.

Etrange. Quel compilateur utilises-tu ?