From 1fcb18ac54ed0aef7d75a8ca2c316d03c30f059e Mon Sep 17 00:00:00 2001 From: Oliver Hahn Date: Mon, 25 Nov 2013 16:18:33 +0100 Subject: [PATCH] fixed a bug caused by non-existing refinement masks not being handled properly in unigrid case --- mesh.hh | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/mesh.hh b/mesh.hh index ae7fa05..58a6378 100644 --- a/mesh.hh +++ b/mesh.hh @@ -939,8 +939,7 @@ public: // 1 = in mask and not refined (i.e. cell exists only on this level) - //if( ilevel == levelmax() ) - // return false; + if( bhave_refmask ){ short v = (*m_ref_masks[ilevel])(i,j,k); @@ -950,6 +949,9 @@ public: //if( ilevel == levelmax()-1 && bhave_refmask ) // return (*m_ref_masks[ilevel])(i,j,k); + + if( !bhave_refmask && ilevel == levelmax() ) + return false; if( i < offset(ilevel+1,0) || i >= offset(ilevel+1, 0)+(int)size(ilevel+1,0)/2 || j < offset(ilevel+1,1) || j >= offset(ilevel+1, 1)+(int)size(ilevel+1,1)/2 ||