diff options
author | Desire NUENTSA <desire.nuentsa_wakam@inria.fr> | 2013-01-11 17:16:14 +0100 |
---|---|---|
committer | Desire NUENTSA <desire.nuentsa_wakam@inria.fr> | 2013-01-11 17:16:14 +0100 |
commit | 91b3b3aaab19fc11db18d95a28c1a0be9ae9d9cd (patch) | |
tree | d5858697f227ce67812944a181627a593df76552 /Eigen/src/MetisSupport | |
parent | 1ccd90a927e7386574ff845ff0d326733352e9d1 (diff) |
Add a sparse QR factorization and update the elimination tree in SparseLU
Diffstat (limited to 'Eigen/src/MetisSupport')
-rw-r--r-- | Eigen/src/MetisSupport/MetisSupport.h | 3 |
1 files changed, 1 insertions, 2 deletions
diff --git a/Eigen/src/MetisSupport/MetisSupport.h b/Eigen/src/MetisSupport/MetisSupport.h index a762d96f6..3a723b384 100644 --- a/Eigen/src/MetisSupport/MetisSupport.h +++ b/Eigen/src/MetisSupport/MetisSupport.h @@ -122,10 +122,9 @@ public: //NOTE: If Ap is the permuted matrix then perm and iperm vectors are defined as follows // Row (column) i of Ap is the perm(i) row(column) of A, and row (column) i of A is the iperm(i) row(column) of Ap - // To be consistent with the use of the permutation in SparseLU module, we thus keep the iperm vector matperm.resize(m); for (int j = 0; j < m; j++) - matperm.indices()(j) = iperm(j); + matperm.indices()(iperm(j)) = j; } |