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

Ensure particle direction is normalized for plotting / volume calculations (#2816)

Browse source
This commit is contained in:
Paul Romano 2023-12-28 14:37:21 -06:00 • committed by GitHub
parent 411656668f
commit 2e4811b26b
No known key found for this signature in database
GPG key ID: 4AEE18F83AFDEB23
2 changed files with 3 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

3
include/openmc/plot.h
View file

@ -1,6 +1,7 @@
#ifndef OPENMC_PLOT_H
#define OPENMC_PLOT_H
#include <cmath>
#include <sstream>
#include <unordered_map>
@ -200,7 +201,7 @@ T SlicePlotBase::get_map() const
xyz[out_i] = origin_[out_i] + width_[1] / 2. - out_pixel / 2.;
// arbitrary direction
Direction dir = {0.7071, 0.7071, 0.0};
Direction dir = {1. / std::sqrt(2.), 1. / std::sqrt(2.), 0.0};
#pragma omp parallel
{

2
src/volume_calc.cpp
View file

@ -159,7 +159,7 @@ vector<VolumeCalculation::Result> VolumeCalculation::execute() const
p.n_coord() = 1;
Position xi {prn(&seed), prn(&seed), prn(&seed)};
p.r() = lower_left_ + xi * (upper_right_ - lower_left_);
p.u() = {0.5, 0.5, 0.5};
p.u() = {1. / std::sqrt(3.), 1. / std::sqrt(3.), 1. / std::sqrt(3.)};
// If this location is not in the geometry at all, move on to next block
if (!exhaustive_find_cell(p))