From cf09ef9022498fe6126ffc98305e9e7f1fc670a3 Mon Sep 17 00:00:00 2001
From: Richard Comblen <richard.comblen@uclouvain.be>
Date: Thu, 4 Mar 2010 13:18:15 +0000
Subject: [PATCH] Removed bug in splitMultirate

---
 Solver/dgGroupOfElements.cpp | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/Solver/dgGroupOfElements.cpp b/Solver/dgGroupOfElements.cpp
index 854e23e1aa..b4f8c979d1 100644
--- a/Solver/dgGroupOfElements.cpp
+++ b/Solver/dgGroupOfElements.cpp
@@ -1170,7 +1170,7 @@ double dgGroupCollection::splitGroupsForMultirate(int maxLevels,dgConservationLa
           int gId=interface.connectedElements[iConn].first;
           int eId=interface.connectedElements[iConn].second;
           int newGroupId=newGroupIds[gId][eId];
-          if(newGroupId >= _lowerLevelGroupIdStart && newGroupId<_lowerLevelGroupIdEnd){
+          if(newGroupId>=0 /*newGroupId >= _lowerLevelGroupIdStart && newGroupId<_lowerLevelGroupIdEnd*/){
             toAdd=true;
             continue;
           }
@@ -1246,6 +1246,9 @@ double dgGroupCollection::splitGroupsForMultirate(int maxLevels,dgConservationLa
             currentNewGroupId++;
           }
 
+          for(int i=0;i<forInnerBuffer.size();i++){
+            newGroupIds[it->first][forInnerBuffer[i]]=currentNewGroupId;
+          }
           if(!forInnerBuffer.empty()){
             std::vector<MElement*>forInnerBufferV;
             for(int i=0;i<forInnerBuffer.size();i++){
-- 
GitLab