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

avoid unneeded const_cast

Browse Source
This commit is contained in:
Desire NUENTSA 2013-07-16 15:56:05 +02:00
parent 3e094af410
commit cfd7f9b84a
1 changed files with 10 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
Eigen/src/SparseLU/SparseLU.h
View File

@ -377,12 +377,13 @@ void SparseLU<MatrixType, OrderingType>::analyzePattern(const MatrixType& mat)
if (m_perm_c.size()) { if (m_perm_c.size()) {
m_mat.uncompress(); //NOTE: The effect of this command is only to create the InnerNonzeros pointers. FIXME : This vector is filled but not subsequently used. m_mat.uncompress(); //NOTE: The effect of this command is only to create the InnerNonzeros pointers. FIXME : This vector is filled but not subsequently used.
//Then, permute only the column pointers //Then, permute only the column pointers
Index * outerIndexPtr; const Index * outerIndexPtr;
if (mat.isCompressed()) outerIndexPtr = const_cast<Index *>(mat.outerIndexPtr()); if (mat.isCompressed()) outerIndexPtr = mat.outerIndexPtr();
else else
{ {
outerIndexPtr = new Index[mat.cols()+1]; Index *outerIndexPtr_t = new Index[mat.cols()+1];
for(Index i = 0; i <= mat.cols(); i++) outerIndexPtr[i] = m_mat.outerIndexPtr()[i]; for(Index i = 0; i <= mat.cols(); i++) outerIndexPtr_t[i] = m_mat.outerIndexPtr()[i];
outerIndexPtr = outerIndexPtr_t;
} }
for (Index i = 0; i < mat.cols(); i++) for (Index i = 0; i < mat.cols(); i++)
{ {
@ -461,12 +462,13 @@ void SparseLU<MatrixType, OrderingType>::factorize(const MatrixType& matrix)
{ {
m_mat.uncompress(); //NOTE: The effect of this command is only to create the InnerNonzeros pointers. m_mat.uncompress(); //NOTE: The effect of this command is only to create the InnerNonzeros pointers.
//Then, permute only the column pointers //Then, permute only the column pointers
Index * outerIndexPtr; const Index * outerIndexPtr;
if (matrix.isCompressed()) outerIndexPtr = const_cast<Index *>(matrix.outerIndexPtr()); if (matrix.isCompressed()) outerIndexPtr = matrix.outerIndexPtr();
else else
{ {
outerIndexPtr = new Index[matrix.cols()+1]; Index* outerIndexPtr_t = new Index[matrix.cols()+1];
for(Index i = 0; i <= matrix.cols(); i++) outerIndexPtr[i] = m_mat.outerIndexPtr()[i]; for(Index i = 0; i <= matrix.cols(); i++) outerIndexPtr_t[i] = m_mat.outerIndexPtr()[i];
outerIndexPtr = outerIndexPtr_t;
} }
for (Index i = 0; i < matrix.cols(); i++) for (Index i = 0; i < matrix.cols(); i++)
{ {