remove temporary in SparseLU::solve

(grafted from 62acd679037400b55f81aac709edea75dc5ae876
)
This commit is contained in:
Gael Guennebaud 2016-12-05 15:11:57 +01:00
parent a9bb9796e0
commit 28bf2bf070

View File

@ -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++)