boost::numeric::ublas matrix.assign( identity_matrix )

Code:

---

1 2 3 4 5 6

m_Matrice = matrix<float>(4,4); m_Identite = identity_matrix<float>(4); m_Zero = zero_matrix<float>(4,4);   m_Matrice.assign(m_Identite);

---

Problème si j'affiche m_Matrice apres ça , on vois qu'elle est composée que de zéros.
J'ai alors affecté à la main

Code:

---

1 2 3 4 5

m_Matrice(0,0) = 1.0; m_Matrice(1,1) = 1.0; m_Matrice(2,2) = 1.0; m_Matrice(3,3) = 1.0;

---

pour voir si cela venait de l'affichage . Ce qui n'est pas le cas. J'ai alors affiché m_Identite et on voit bien la diagonale de 1.
J'ai essayé avec

Code:

---

1 2	m_Matrice = m_Identite

---

Cela non plus ne fonctionne pas. Et avec l'opérateur += l'appliaction plante.

cela fonctionne bien si je fais

Code:

---

1 2	m_Matrice = identity_matrix<float>(4);

---

mais bon j'aimerais bien assigner des matrices existantes entre elles.

Documention ublas
Documention du type matrix

http://osdir.com/ml/lib.boost.ublas/.../msg00106.html
Comme on peut le voir sur ce lien pour créer une matrice d'identitée, on lui passe en paramètres deux tailles

Code:

---

1 2	identity_matrix<value_type>( a.size1(), a.size2() )

---

Hors je croyais en regardant dans la doc de ublas que la création se faisait avec une seule taille puisqu'une matrice d'identitée est carrée.

Code:

---

1 2	identity_matrix (size_type size) Constructs an identity_matrix that holds size rows of size elements.

---

Ce soir ou ce midi je vais donc essayer avec deux paramètres.

Faites vos matrices avec des pointeurs c'est moche pour les opérateurs mais quand on créer un objet avec des templates vos mieux le faire nikel du premier cou que d'affecter un objet par copie lors de la construction. C'est la conclusion que j'en tire .

EDIT:

Code:

---

1 2	m_Identite = new identity_matrix<float>(4,4); //ca marche tres bien aussi

---