diff options
author | Desire NUENTSA <desire.nuentsa_wakam@inria.fr> | 2012-09-03 14:49:03 +0200 |
---|---|---|
committer | Desire NUENTSA <desire.nuentsa_wakam@inria.fr> | 2012-09-03 14:49:03 +0200 |
commit | 3a22c47fb5ab1d34498fdaa5b1263c1426d97a37 (patch) | |
tree | 59551e9411eb2fef3db3ed922f31ec2d1fb3a938 /Eigen | |
parent | 63d2dcfb7045ce4346b67e09c5c836e49fb75c6e (diff) |
Bug in blas 3 2D block update
Diffstat (limited to 'Eigen')
-rw-r--r-- | Eigen/src/SparseLU/SparseLU_panel_bmod.h | 7 |
1 files changed, 3 insertions, 4 deletions
diff --git a/Eigen/src/SparseLU/SparseLU_panel_bmod.h b/Eigen/src/SparseLU/SparseLU_panel_bmod.h index bfe13b38a..90d3297a1 100644 --- a/Eigen/src/SparseLU/SparseLU_panel_bmod.h +++ b/Eigen/src/SparseLU/SparseLU_panel_bmod.h @@ -95,7 +95,7 @@ void LU_panel_bmod(const int m, const int w, const int jcol, const int nseg, Sca // if the blocks are large enough, use level 3 // TODO find better heuristics! - if(false && nsupc >= 50 && nrow > 50 && u_cols>6) + if(true /*nsupc >= 50 && nrow > 50 && u_cols>6*/) { Map<Matrix<Scalar,Dynamic,Dynamic> > U(tempv.data(), u_rows, u_cols); @@ -126,7 +126,6 @@ void LU_panel_bmod(const int m, const int w, const int jcol, const int nseg, Sca u_col++; } - // solve U = A^-1 U luptr = glu.xlusup(fsupc); no_zeros = (krep - u_rows + 1) - fsupc; @@ -162,14 +161,14 @@ void LU_panel_bmod(const int m, const int w, const int jcol, const int nseg, Sca { int irow = glu.lsub(isub++); dense_col(irow) = U.coeff(i+off,u_col); - U.coeffRef(i,u_col) = 0; + U.coeffRef(i+off,u_col) = 0; } // Scatter l into SPA dense[] for (int i = 0; i < nrow; i++) { int irow = glu.lsub(isub++); - dense_col(irow) -= L.coeff(i+off,u_col); + dense_col(irow) -= L.coeff(i,u_col); L.coeffRef(i,u_col) = 0; } u_col++; |