diff --git a/apps/cpu/LaminarPlaneFlow/LaminarPlaneFlow.cfg b/apps/cpu/LaminarPlaneFlow/LaminarPlaneFlow.cfg index 161c21057c9e343b79399d7a9e043186d3b67381..a767d4d2a6cca4a697ea6e6fb13b03c39c43b0ed 100644 --- a/apps/cpu/LaminarPlaneFlow/LaminarPlaneFlow.cfg +++ b/apps/cpu/LaminarPlaneFlow/LaminarPlaneFlow.cfg @@ -1,15 +1,15 @@ pathname = ./output/LaminarPlaneFlow -numOfThreads = 18 +numOfThreads = 1 blocknx = 16 16 1 -boundingBox = 32 32 1 +boundingBox = 64 64 1 deltax = 1 -refineLevel = 0 +refineLevel = 1 newStart = true restartStep = 1000000 -cpStep = 1000000 -cpStepStart = 1000000 +cpStep = 1000000 +cpStart = 1000000 outTime = 10 endTime = 10 diff --git a/src/cpu/core/Visitors/SetInterpolationConnectorsBlockVisitor.cpp b/src/cpu/core/Visitors/SetInterpolationConnectorsBlockVisitor.cpp index 4f882353f1ce6d7675271e753bd1ebf9f9962c8e..691c65ba8c5006b6e2ec89968e545621692ee1ab 100644 --- a/src/cpu/core/Visitors/SetInterpolationConnectorsBlockVisitor.cpp +++ b/src/cpu/core/Visitors/SetInterpolationConnectorsBlockVisitor.cpp @@ -59,6 +59,10 @@ void SetInterpolationConnectorsBlockVisitor::visit(SPtr<Grid3D> grid, SPtr<Block { if(!block) return; + UbTupleInt3 blockNX = grid->getBlockNX(); + if (val<1>(blockNX) % 2 == 0 && val<2>(blockNX) % 2 == 0 && val<3>(blockNX) % 2 == 0) + UB_THROW(UbException(UB_EXARGS, "Odd number of nodes: The number of nodes in each dimension (x,y,z) has to be even.")); + UBLOG(logDEBUG5, "SetInterpolationConnectorsBlockVisitor::visit() - start"); UBLOG(logDEBUG5, block->toString());