Skip to content
Snippets Groups Projects
Commit dcccc18b authored by Oliver Sander's avatar Oliver Sander Committed by sander
Browse files

Loop over algebraic entities instead of grid entities after all

[[Imported from SVN: r9764]]
parent b72d8940
No related branches found
No related tags found
No related merge requests found
...@@ -28,9 +28,8 @@ private: ...@@ -28,9 +28,8 @@ private:
std::vector<TransferVectorTuple> localVectorEntries; std::vector<TransferVectorTuple> localVectorEntries;
// Translate vector entries // Translate vector entries
const auto& gridView = guIndex.getGridView(); for (int k=0; k<localVector.size(); k++)
for (auto it = gridView.template begin<2>(); it != gridView.template end<2>(); ++ it) localVectorEntries.push_back(TransferVectorTuple(guIndex.globalIndex(k), localVector[k]));
localVectorEntries.push_back(TransferVectorTuple(guIndex.globalIndex(*it), localVector[guIndex.localIndex(*it)]));
// Get number of vector entries on each process // Get number of vector entries on each process
localVectorEntriesSizes = MPIFunctions::shareSizes(guIndex.getGridView(), localVectorEntries.size()); localVectorEntriesSizes = MPIFunctions::shareSizes(guIndex.getGridView(), localVectorEntries.size());
...@@ -89,7 +88,7 @@ public: ...@@ -89,7 +88,7 @@ public:
// And translate solution again // And translate solution again
for (size_t k = 0; k < localVectorEntries.size(); ++k) for (size_t k = 0; k < localVectorEntries.size(); ++k)
x[guIndex.localIndex(localVectorEntries[k].globalIndex_)] = localVectorEntries[k].value_; x[k] = localVectorEntries[k].value_;
return x; return x;
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment