Skip to content
Snippets Groups Projects
Commit 3270b21b authored by Elias Pipping's avatar Elias Pipping Committed by Elias Pipping
Browse files

Cleanup

parent 4e61ec10
No related branches found
No related tags found
No related merge requests found
...@@ -188,19 +188,13 @@ template <class MyConvexProblemTypeTEMPLATE> class MyBlockProblem { ...@@ -188,19 +188,13 @@ template <class MyConvexProblemTypeTEMPLATE> class MyBlockProblem {
linearization.b *= -1.0; linearization.b *= -1.0;
// apply truncation to system // apply truncation to system
typename Linearization::MatrixType::row_type::Iterator it;
typename Linearization::MatrixType::row_type::Iterator end;
for (size_t row = 0; row < linearization.A.N(); ++row) { for (size_t row = 0; row < linearization.A.N(); ++row) {
it = linearization.A[row].begin(); auto const end = linearization.A[row].end();
end = linearization.A[row].end(); for (auto it = linearization.A[row].begin(); it != end; ++it) {
for (; it != end; ++it) {
int const col = it.index(); int const col = it.index();
for (size_t i = 0; i < it->N(); ++i) { for (size_t i = 0; i < it->N(); ++i) {
typename Linearization::MatrixType::block_type::row_type::Iterator auto const blockEnd = (*it)[i].end();
blockIt = (*it)[i].begin(); for (auto blockIt = (*it)[i].begin(); blockIt != blockEnd; ++blockIt)
typename Linearization::MatrixType::block_type::row_type::
Iterator const blockEnd = (*it)[i].end();
for (; blockIt != blockEnd; ++blockIt)
if (linearization.truncation[row][i] or linearization if (linearization.truncation[row][i] or linearization
.truncation[col][blockIt.index()]) .truncation[col][blockIt.index()])
*blockIt = 0.0; *blockIt = 0.0;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment