Skip to content
Snippets Groups Projects
Commit f85f762c authored by Oliver Sander's avatar Oliver Sander Committed by sander@FU-BERLIN.DE
Browse files

fixed index bug 

[[Imported from SVN: r8257]]
parent 1f4aa4e8
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
dune/gfe/pktop1mgtransfer.hh
+ 3
− 3
View file @ f85f762c
...@@ -97,7 +97,7 @@ public: ...@@ -97,7 +97,7 @@ public:
for (size_t j=0; j<values.size(); j++) { for (size_t j=0; j<values.size(); j++) {
if (values[i] > 0.001) if (values[j] > 0.001)
{ {
size_t globalFine = fineBasis.index(*it, j); size_t globalFine = fineBasis.index(*it, j);
size_t globalCoarse = p1Basis.index(*it, i); size_t globalCoarse = p1Basis.index(*it, i);
...@@ -135,12 +135,12 @@ public: ...@@ -135,12 +135,12 @@ public:
for (size_t j=0; j<values.size(); j++) { for (size_t j=0; j<values.size(); j++) {
if (values[i] > 0.001) if (values[j] > 0.001)
{ {
size_t globalFine = fineBasis.index(*it, j); size_t globalFine = fineBasis.index(*it, j);
size_t globalCoarse = p1Basis.index(*it, i); size_t globalCoarse = p1Basis.index(*it, i);
(*this->matrix_)[globalFine][globalCoarse] = Dune::ScaledIdentityMatrix<double,TransferMatrixBlock::rows>(values[i]); (*this->matrix_)[globalFine][globalCoarse] = Dune::ScaledIdentityMatrix<double,TransferMatrixBlock::rows>(values[j]);
} }
} }
......
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