Discussion :

[salseropom]

Bonjour, après deux ans d'abandon du C++ (pour faire du C) je me remet au C++. J'avoue que j'ai dû oublié de nombreuses choses. Je m'entraine donc à refaire une classe Matrix templatée et c'est un peu plus dur que je ne le pensais. J'ai un pb de point virgule que je n'arrive pas à résoudre. Voici mon code :

main.cpp

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 <iostream>
//#include "matrix.h"
 
using namespace std;
 
int main(void)
{
  std::cout<<"Hello World"<<std::endl;
  /*
  Matrix<double> A;
  Matrix<double> B(2,2,2);
  //cout<<B<<endl;
  
  Matrix<Matrix<double> > C(4,3,B);
  //cout<<C<<endl;
*/
  return 0;//EXIT_SUCCESS;
}

matrix.cpp

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
 
#include<iostream>
#include "matrix.h"
 
/* default constructor */
template<typename T>
Matrix<T>::Matrix()
{
  std::cout<<"appel du constructeur par default"<<std::endl;
} 
 
  //Matrix(const Matrix &); /* copy constructor */
 
template<typename T>
Matrix<T>::Matrix(unsigned row,unsigned col,const T & d)
{
  std::cout<<"appel du constructeur"<<std::endl;
  nr=row;
  nc=col;
  data=new T * [ nr ];
 
  unsigned i,j;
 
  for(i=0;i<nr;i++)
  {
     data[i]=new T [ nc ];
     if(data[i]==NULL) std::cerr<<"Error"<<std::endl;
     for(j=0;j<nc;++j) data[i][j]=d;
  }
 
}

matrix.h

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
 
#ifndef MATRIX_H_PLP_20070925
#define MATRIX_H_PLP_20070925
 
#include<iostream>
 
template<typename T> class Matrix
{
private:
  unsigned nr; /* number of rows    */
  unsigned nc; /* number of columns */
  T ** data;
 
public:
  Matrix(); /* default constructors */
 
  Matrix(unsigned,unsigned, const T &); /* constructor */
 
 
}
 
#endif

et voici mon erreur :

../Sources/matrix.cpp:5: error: semicolon missing after declaration of `
Matrix<T>'
../Sources/matrix.cpp:5: error: erreur de syntaxe before `template'
make: *** [../Sources/matrix.o] Erreur 1

je suis un peu séché.

Autre question : j'ai dans mon privaite : T ** data. Ne vaudrait-il pas mieux faire un

Code :

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

1
2
 
vector<vector<T> > data;

ainsi, pour une augmentation ou diminution de ma matrice, cela sera plus facile non ?

Merci bcp

[rikau2]

[/CODE]

matrix.h

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
 
#ifndef MATRIX_H_PLP_20070925
#define MATRIX_H_PLP_20070925
 
#include<iostream>
 
template<typename T> class Matrix
{
private:
  unsigned nr; /* number of rows    */
  unsigned nc; /* number of columns */
  T ** data;
 
public:
  Matrix(); /* default constructors */
 
  Matrix(unsigned,unsigned, const T &); /* constructor */
 
 
}; // Elle est ici l'erreur
 
#endif

Il manque en effet un ; a la fin de la déclaration de ta classe.

Autre question : j'ai dans mon privaie : T ** data. Ne vaudrait-il pas mieux faire un

Code :

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

1
2
 
vector<vector<T> > data;

ainsi, pour une augmentation ou diminution de ma matrice, cela sera plus facile non ?

Merci bcp

Pour cela je te conseille de regarder les differents containers possibles.... Dans la FAQ est bien décrit les différents avantages et inconvénients de tel ou tel conteneur (set, vector, map, ...)

Pour bien gérer ce type de conteneur, je ne peux que te conseiller ce site:
CplusPlus Reference

Bonne reprise.

[Médinoc]

Pour la première question : Chaque déclaration de classe ou de structure doit être suivie d'un point-virgule:

Code :

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

1
2
3
4
5
6
7
8
#ifndef H_TRUC
#define H_TRUC

struct truc
{
	int a;
};
#endif

[salseropom]

Pour la première question : Chaque déclaration de classe ou de structure doit être suivie d'un point-virgule:

Code :

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

1
2
3
4
5
6
7
8
#ifndef H_TRUC
#define H_TRUC

struct truc
{
	int a;
};
#endif

OK, merci. Honte à moi !

[Médinoc]

Pour la seconde question, je dirais que non. À bas niveau comme pour ta classe, je pense qu'utiliser vector n'est pas approprié, et que le changement de dimensions passe de toute façon par la copie des éléments de l'ancienne matrice à la nouvelle.

PS: Pour la création, je serais plus du genre à faire une allocation par dimension, plutôt que nr+1 allocations...

[salseropom]

Pour la seconde question, je dirais que non. À bas niveau comme pour ta classe, je pense qu'utiliser vector n'est pas approprié, et que le changement de dimensions passe de toute façon par la copie des éléments de l'ancienne matrice à la nouvelle.

PS: Pour la création, je serais plus du genre à faire une allocation par dimension, plutôt que nr+1 allocations...

Oui, je comprends bien ce que tu veux dire mais si je fais qqch du genre

Code :

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

1
2
3
4
5
6
7
8
9
 
class Matrix
{
private :
  unsigned nr,nc;
 T * data;
public :
 blabla;
}

avec un

Code :

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

1
2
 
data=new T[nc*nr];

pour accéder à l'élément A(i,j), je vais devoir faire

Code :

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

1
2
 
data[i*nc+j]

cad 2 opérations, alors qu'avec nc+1 allocations, je ne ferais que

Code :

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

1
2
 
data[i][j]

soit 0 opération. Et pour les produits matricels, il y a N^3 opérations environs (matrice carrée de taille N)

[rikau2]

Code :

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

1
2
3
4
5
6
7
8
9
 
class Matrix
{
private :
  unsigned nr,nc;
 T * data;
public :
 blabla;
}

Attention il manque encore un ";"...!
Pour ce qui est de l'accès aux éléments de ta matrice, il faut que tu réfléchisses comment se feront ces dit accès... Je m'explique:
Si ta matrice est énorme et que tu accède les éléments de manière successive tu risque d'avoir des page fault en fonction de l'ordre dàccés. Donc si le tableau est a deux dimensions, réfléchir à l'ordre de stockage pour utiliser le "buffer" de lecture du HDD.

[screetch]

early optimization is the root of all evil. comemnce par faire un truc qui marche, ensuite tu pourras commencer a comparer les methodes et voir ce qui est le mieux.

[Médinoc]

Oui, je comprends bien ce que tu veux dire mais si je fais qqch du genre

Code :

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

1
2
3
4
5
6
7
8
class Matrix
{
private :
  unsigned nr,nc;
 T * data;
public :
 blabla;
};

avec un

Code :

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

1
2
 
data=new T[nc*nr];

<snip les conséquences d'un tel design>

Non.
Je parlais de ceci:

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
class Matrix
{
private:
	unsigned nr,nc;
	T** pRowPointers;
public:
	blabla;
};
 
...
 
Matrix::Matrix(nr,nc etc.)
{
	pRowPointers = new T* [nr];
	T * pCells = new T[nr * nc]
	for(i=0 ; i<nr ; i++)
		pRowPointers[i] = &(pCells[i*nc]);
}
 
Matrix::~Matrix()
{
	T * pCells = pRowPointers[0];
	delete[] pCells;
	delete[] pRowPointers;
}

C'est ça, une allocation par dimension...
On réduit drastiquement le nombre d'allocations sans augmenter les temps d'accès...

[salseropom]

Non.
Je parlais de ceci:

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
class Matrix
{
private:
	unsigned nr,nc;
	T** pRowPointers;
public:
	blabla;
};
 
...
 
Matrix::Matrix(nr,nc etc.)
{
	pRowPointers = new T* [nr];
	T * pCells = new T[nr * nc]
	for(i=0 ; i<nr ; i++)
		pRowPointers[i] = &(pCells[i*nc]);
}
 
Matrix::~Matrix()
{
	T * pCells = pRowPointers[0];
	delete[] pCells;
	delete[] pRowPointers;
}

C'est ça, une allocation par dimension...
On réduit drastiquement le nombre d'allocations sans augmenter les temps d'accès...

Ah oui, c'est très malin ce que tu viens de me dire !!
Je vais bien y réfléchir.

Merci encore.