Discussion :

[vincho]

Bonjour,
Je m'intéresse à la meta-prog pour optimiser des calculs entre des vecteurs(class Point). Cependant je remarque que le calcul fait sur une recursion et specialisation pour boucler va deux moins vite que d'écrire directement le code.
Est-ce normal? Ai-je mal codé?
Merci de vos réponse.

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
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
 
template <int Begin, int End> struct Boucle
{
  template<class T>
  static void Add( T & p1, const T &p2)
    {
      p1[Begin]+=p2[Begin];
      Boucle<Begin+1,End>::Add(p1,p2);
    }
};
 
template <int N> struct Boucle<N, N>
{
  template<class T>
  static void Add( T & p1, const T &p2) {}
};
 
template <int dim, class T>
class Point
{
public:
 
  T dat[dim];
  inline Point<dim,T> & Add(const Point<dim,T>& p)
    {
      Boucle<0, dim>::Add(*this, p);
      return *this;
    }
    inline Point<dim,T> & Addf(const Point<dim,T>& p)
    {
      (*this)[0]+=p[0];
      (*this)[1]+=p[1];
      (*this)[2]+=p[2];
      (*this)[3]+=p[3];
      return *this;
    }
  inline const T & operator []( int i)
    const
    {
      return   this->dat[i];
    }
  inline  T & operator []( int i)
 
    {
      return   this->dat[i];
    }
};
 
int main(void){
 
  Point<4,int > p1,p2;
 
  int t1=time(NULL);
  for(int i=0;i<100000;i++)
    for(int j=0;j<1000;j++)
  {
    p1.Add(p2);
  }
  int t2=time(NULL);
  cout<< t2-t1<<endl;
 
  t1=time(NULL);
  for(int i=0;i<100000;i++)
    for(int j=0;j<1000;j++)
  {
    p1.Addf(p2);
  }
 t2=time(NULL);
 cout<< t2-t1<<endl;
}

Affichage console:

Code :

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

1
2
3
 
8
4

[PierroElMito]

tu compiles en release?
théoriquement, c'est censé donner le même resultat. essaie peut-etre de rajouter des __forceinline dans ta struct Boucle. sinon, il faut regarder le code assembleur généré pour comprendre la différence.

[vincho]

tu compiles en release?

je compile avec g++-4.0

Code :

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

1
2
 
g++-4.0 _rm.cpp -o exe

Je n'ai que des vagues souvenirs d'assembleur....

Dans l'article de Gomila, il ne parle pas de _Forceinline.
Je pense que le compilateur est assez fort pour générer un bon code.

[PierroElMito]

je connais pas trop gcc, mais je dirais qu'il faut rajouter -O3 a ta ligne de compilation pour activer les optimisations.
pour __forceinline, c'est spécifique à visual c++ me semble t'il, donc ça marchera pas dans ton cas. il faut quand même peut-etre mettre juste 'inline'.

[vincho]

je connais pas trop gcc, mais je dirais qu'il faut rajouter -O3 a ta ligne de compilation pour activer les optimisations.

J'ai multiplié par 50 le nombres de tours de boucles:
J'obtiens:

Code :

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

1
2
3
 
17
16

Donc à peu de chose près la même chose.

Il faudrait dire à Gomila de rajouter les spécifications de compilation dans son article

Merci PierroElMito

[Laurent Gomila]

Il faudrait dire à Gomila de rajouter les spécifications de compilation dans son article

Ca n'a pas vraiment de rapport avec la méta-programmation, ce serait un algo de tri ou un jeu d'échecs ce serait pareil : pas de comparaison valable sans option d'optimisation.
Pour ça il faut plutôt regarder les tutoriels du genre "bien utiliser son environnement de développement"

Et pour ceux qui ont un petit background ASM, regarder le code généré sera toujours plus parlant et fiable qu'un benchmark.

[vincho]

Intérêt par rapport à un code classique

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
  t1=time(NULL);
  for(int i=0;i<100000;i++)
    for(int j=0;j<50000;j++)
  {
    for(int k=0;k<4;k++)
    {
      p1[k]+=p2[k];
    }
  }
 t2=time(NULL);
 cout<< t2-t1<<endl;

Le benchmark est de 46s donc on a gagné plus de 3fois en vitesse.

Ce que je trouve étrange est pourquoi mon compilateur avec l'option -O3 ou moins -O4 ne déroule pas la boucle???

[PierroElMito]

il faut peut-etre rajouter le parametre:
-funroll-loops
a ta ligne de compilation.