Skip to content
Snippets Groups Projects
Commit 7b6291f5 authored by Samuel Kemmler's avatar Samuel Kemmler
Browse files

Add fix for NANs after starting from checkpoint

parent 2b2543d4
No related merge requests found
......@@ -332,11 +332,11 @@ private:
{
if( singleCast(particleIdx, ac, containsPointFctr, ac, Vector3<real_t>(cx,cy,cz)) )
{
if(boundaryHandling->isDomain(x,y,z))
{
//if(boundaryHandling->isDomain(x,y,z))
//{
boundaryHandling->forceBoundary(obstacleFlag, x, y, z);
(*particleField)(x,y,z) = ac.getUid(particleIdx);
}
//}
// no else -> will not overwrite pre-existing boundary cells, e.g. of previously applied applied mapping (e.g. when first mapping global particles)
}
cx += dx;
......
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