diff --git a/Solver/dgSlopeLimiter.cpp b/Solver/dgSlopeLimiter.cpp index f4ce4cb54ba0c732525ebb2fc39c83f21c18a9ff..3490c0f196f194afc7a53db5a15980842b43b757 100644 --- a/Solver/dgSlopeLimiter.cpp +++ b/Solver/dgSlopeLimiter.cpp @@ -104,13 +104,8 @@ bool dgSlopeLimiter::apply ( dgDofContainer &solution, if (solutionEClipped){ for (int iElement=0 ; iElement<egroup->getNbElements() ;++iElement) { solutionE.setAsProxy(solGroup, iElement*nbFields, nbFields ); - fullMatrix<double> Temp; - Temp.setAsProxy(solGroup, nbFields*iElement, nbFields ); - cacheElement.set(egroup->getElement(iElement)); - for (int K=0;K<Temp.size1();K++) - for (int L=0;L<Temp.size2();L++) - Temp(K,L) = (*solutionEClipped)(K,L); - } + solutionE.set((*solutionEClipped)()); + } delete solutionEClipped; } }