Commit 39c4d613 authored by Alvarez, Gonzalo's avatar Alvarez, Gonzalo
Browse files

LanczosVectors: swap, don't copy

parent a12147c6
......@@ -268,13 +268,13 @@ public:
SizeType steps() const {return steps_; }
const DenseMatrixType& lanczosVectors() const
void lanczosVectorsSwap(DenseMatrixType& V)
{
const DenseMatrixType* ptr = lanczosVectors_.data();
DenseMatrixType* ptr = lanczosVectors_.data();
if (!ptr)
throw RuntimeError("LanczosCore::lanczosVectors() called but no data stored\n");
return *(ptr);
ptr->swap(V);
}
const MatrixType& matrix() const { return mat_; }
......
......@@ -91,9 +91,9 @@ public:
return ls_.decomposition(initVector, ab, 1);
}
const typename LanczosCoreType::DenseMatrixType& lanczosVectors() const
void lanczosVectorsSwap(typename LanczosCoreType::DenseMatrixType& V)
{
return ls_.lanczosVectors();
ls_.lanczosVectorsSwap(V);
}
SizeType steps() const {return ls_.steps(); }
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment