crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 14:15:42 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Adding coincidence check when comparing d_surf, d_lat.

Browse source
This commit is contained in:
Patrick Shriwise 2019-03-30 15:42:11 -05:00
parent c69d53ee86
commit e62dd935a6
1 changed files with 5 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

6
src/geometry.cpp
View file

@ -32,6 +32,10 @@ std::vector<int64_t> overlap_check_count;
// Non-member functions
//==============================================================================
bool coincident(const double& d1, const double& d2) {
return std::abs(d1 - d2) < FP_COINCIDENT;
}
bool check_cell_overlap(Particle* p)
{
int n_coord = p->n_coord_;
@ -396,7 +400,7 @@ BoundaryInfo distance_to_boundary(Particle* p)
// a higher level then we need to make sure that the higher level boundary
// is selected. This logic must consider floating point precision.
double& d = info.distance;
if (d_surf < d_lat) {
if (d_surf < d_lat && !coincident(d_surf, d_lat)) {
if (d == INFINITY || (d - d_surf)/d >= FP_REL_PRECISION) {
d = d_surf;