crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-29 06:35:48 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Don't initially reset surface when doing exhaustive find cell search

Browse source
This commit is contained in:
Paul Romano 2022-05-26 10:01:15 -05:00
parent 67d5449127
commit 256078ef10
1 changed files with 2 additions and 6 deletions
Download patch file Download diff file Expand all files Collapse all files

8
src/particle.cpp
View file

@ -454,9 +454,7 @@ void Particle::cross_surface()
// ==========================================================================
// COULDN'T FIND PARTICLE IN NEIGHBORING CELLS, SEARCH ALL CELLS
// Remove lower coordinate levels and assignment of surface
auto surface_index = surface();
surface() = 0;
// Remove lower coordinate levels
n_coord() = 1;
bool found = exhaustive_find_cell(*this);
@ -466,6 +464,7 @@ void Particle::cross_surface()
// the particle is really traveling tangent to a surface, if we move it
// forward a tiny bit it should fix the problem.
surface() = 0;
n_coord() = 1;
r() += TINY_BIT * u();
@ -479,9 +478,6 @@ void Particle::cross_surface()
return;
}
}
// Reassign surface to avoid tracking errors
surface() = surface_index;
}
void Particle::cross_vacuum_bc(const Surface& surf)