Modifier Eigen matrice diagonale

Question

J'ai un Eigen::MatrixXd et j'ai besoin de modifier la valeur des éléments dans sa diagonale. En particulier j'ai un autre Eigen::MatrixXd avec une seule colonne et le même nombre de lignes de la première matrice.

J'ai besoin de soustraire à la diagonale de la première matrice la valeur des éléments de la seconde matrice.

Exemple:

A 
1 2 3 
4 5 6 
7 8 9 

B 
1 
1 
1 


A' 
0 2 3 
4 4 6 
7 8 8 

Comment puis-je faire?

Answer 1

Cela fonctionne pour moi:

A_2=A-B.asDiagonal(); 

Answer 2

for(int i = 0; i < matrix1.rows(); ++i) 
    matrix1(i, i) -= matrix2(i, 0); 

Ce code effectue une itération sur chaque rangée de la matrice (matrix1.rows()) et soustrait la valeur correspondante de matrix2 (matrix2(i, 0)) à partir des diagonales dans matrix1 (matrix1(i, i)).

Answer 3

manipulation de la matrice dans Eigen fonctionne très similaire à celle des tableaux. L'indexation commence à partir de zéro et c'est la ligne principale. La documentation (Eigen: The Matrix Class est bien écrit et peut vous aider à résoudre les problèmes futurs.

Pour votre problème posé la solution est donnée ci-dessous.

#include <iostream> 
#include <eigen3/Eigen/Dense> 
using Eigen::MatrixXd; 
int main() 
{ 
    MatrixXd matA(3,3), matB(3,1); 
    matA<<1,2,3, 
    4,5,6, 
    7,8,9; 
    matB<<1,1,1; 
    for(int i=0; i<3;i++) 
    matA(i,i) -= matB(i); 
    std::cout<<matA<<std::endl; 
    return 0; 
} 

Cependant, je voudrais utiliser matrix3D et Vector3D pour la première et la deuxième . matrice, respectivement

Answer 4

le plus simple et plus rapide pour ne parvenir est:

Eigen::MatrixXd A1(3,3), B(3,1), A2; 
... 
A2 = A1; 
A2.diagonal() -= B; 

bien sûr, mieux utiliser le VectorXd type pour les vecteurs (ici pour B), et enfin si B est constante, vous pouvez utiliser les installations de tableau:

A2.diagonal().array() -= 1;