diff options
| author |  Gael Guennebaud <g.gael@free.fr> | 2016-12-05 15:11:57 +0100 |
|---|---|---|
| committer | Gael Guennebaud <g.gael@free.fr> | 2016-12-05 15:11:57 +0100 |
| commit | 62acd679037400b55f81aac709edea75dc5ae876 (patch) | |
| tree | c42443e9e0061b91bf9dd43b7de560f4b0339b8a /Eigen/src/SparseLU | |
| parent | 0db6d5b3f434ae2e0c2e8b78402e062e67e86339 (diff) | |
remove temporary in SparseLU::solve
Diffstat (limited to 'Eigen/src/SparseLU')
| -rw-r--r-- | Eigen/src/SparseLU/SparseLU_SupernodalMatrix.h | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/Eigen/src/SparseLU/SparseLU_SupernodalMatrix.h b/Eigen/src/SparseLU/SparseLU_SupernodalMatrix.h index f0856db85..2aa746f84 100644 --- a/Eigen/src/SparseLU/SparseLU_SupernodalMatrix.h +++ b/Eigen/src/SparseLU/SparseLU_SupernodalMatrix.h @@ -276,7 +276,7 @@ void MappedSuperNodalMatrix<Scalar,Index_>::solveInPlace( MatrixBase<Dest>&X) co // Matrix-vector product new (&A) Map<const Matrix<Scalar,Dynamic,Dynamic, ColMajor>, 0, OuterStride<> > ( &(Lval[luptr+nsupc]), nrow, nsupc, OuterStride<>(lda) ); - work.block(0, 0, nrow, nrhs) = A * U; + work.block(0, 0, nrow, nrhs).noalias() = A * U; //Begin Scatter for (Index j = 0; j < nrhs; j++) |