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

Pass particle history into DAGMC for more efficient normal calls.

Browse source
This commit is contained in:
Patrick Shriwise 2020-07-14 23:02:52 -05:00
parent aba8af458c
commit 6d41c669b0
2 changed files with 16 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

8
include/openmc/position.h
View file

@ -63,6 +63,14 @@ struct Position {
return std::sqrt(x*x + y*y + z*z);
}
inline Position reflect(Position n) {
const double projection = n.dot(*this);
const double magnitude = n.dot(n);
n *= (2.0 * projection / magnitude);
return *this -= n;
}
//! Rotate the position based on a rotation matrix
Position rotate(const std::vector<double>& rotation) const;

12
src/surface.cpp
View file

@ -172,11 +172,9 @@ Surface::reflect(Position r, Direction u, Particle* p) const
// Determine projection of direction onto normal and squared magnitude of
// normal.
Direction n = normal(r);
const double projection = n.dot(u);
const double magnitude = n.dot(n);
// Reflect direction according to normal.
return u -= (2.0 * projection / magnitude) * n;
return u.reflect(n);
}
Direction
@ -279,7 +277,13 @@ Direction DAGSurface::reflect(Position r, Direction u, Particle* p) const
{
Expects(p);
p->history_.reset_to_last_intersection();
p->last_dir_ = Surface::reflect(r, u, p);
moab::ErrorCode rval;
moab::EntityHandle surf = dagmc_ptr_->entity_by_index(2, dag_index_);
double pnt[3] = {r.x, r.y, r.z};
double dir[3];
rval = dagmc_ptr_->get_angle(surf, pnt, dir, &p->history_);
MB_CHK_ERR_CONT(rval);
p->last_dir_ = u.reflect(dir);
return p->last_dir_;
}