GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-08-14 04:35:57 +08:00
Code Issues Packages Projects Releases Wiki Activity

Fixing bug #578. Thanks to Angelos <filiatra@gmail.com>

Browse Source
This commit is contained in:
Christoph Hertzberg 2013-04-03 16:29:16 +02:00
parent c3a6fa03a2
commit 9b33ab62da
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
Eigen/src/IterativeLinearSolvers/IncompleteLUT.h
View File

@ -29,6 +29,7 @@ int QuickSplit(VectorV &row, VectorI &ind, int ncut)
{ {
typedef typename VectorV::RealScalar RealScalar; typedef typename VectorV::RealScalar RealScalar;
using std::swap; using std::swap;
using std::abs;
int mid; int mid;
int n = row.size(); /* length of the vector */ int n = row.size(); /* length of the vector */
int first, last ; int first, last ;
@ -40,9 +41,9 @@ int QuickSplit(VectorV &row, VectorI &ind, int ncut)
do { do {
mid = first; mid = first;
RealScalar abskey = std::abs(row(mid)); RealScalar abskey = abs(row(mid));
for (int j = first + 1; j <= last; j++) { for (int j = first + 1; j <= last; j++) {
if ( std::abs(row(j)) > abskey) { if ( abs(row(j)) > abskey) {
++mid; ++mid;
swap(row(mid), row(j)); swap(row(mid), row(j));
swap(ind(mid), ind(j)); swap(ind(mid), ind(j));