Skip to content

Random Ray Explicit Void Treatment #3299

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
paulromano merged 15 commits into from
Feb 18, 2025
Merged

Random Ray Explicit Void Treatment #3299

Changes from 1 commit
Commits
Show all changes
15 commits
Select commit Hold shift + click to select a range
5d734d9
refactored and I think its all correct. Now to update the source regi…
jtramm Feb 12, 2025
9875f1b
added fields to SourceRegion
jtramm Feb 12, 2025
3e77ad1
calling flux attenuation function now
jtramm Feb 12, 2025
3084798
pass through complete
jtramm Feb 12, 2025
130d4db
ran git clang format
jtramm Feb 12, 2025
d29562b
removed change to adjoint source computation as this is unrelated to …
jtramm Feb 12, 2025
64d84e8
added void estimator test
jtramm Feb 12, 2025
bd6c53d
added check for sigma_t = zero at program init
jtramm Feb 12, 2025
6fe236e
simplified testing script
jtramm Feb 12, 2025
c96e74f
simplified testing script
jtramm Feb 12, 2025
25bd717
ran clang-format v15
jtramm Feb 12, 2025
b8c3d80
removed unneeded MPI allreduce
jtramm Feb 12, 2025
a9e1ec0
merging with develop
jtramm Feb 17, 2025
f8e8fbd
added some comments to explain SourceRegion fields
jtramm Feb 18, 2025
2a4b55b
Use Doxygen format for data member comments
paulromano Feb 18, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
removed change to adjoint source computation as this is unrelated to … 
…current PR
  • Loading branch information
@jtramm
jtramm committed Feb 12, 2025
commit d29562b113d0ec7db3fec87eb7ed4d049cf13869
9 changes: 2 additions & 7 deletions src/random_ray/flat_source_domain.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1074,13 +1074,8 @@ void FlatSourceDomain::set_adjoint_sources(const vector<double>& forward_flux)
#pragma omp parallel for
for (int64_t sr = 0; sr < n_source_regions_; sr++) {
for (int g = 0; g < negroups_; g++) {
// Protect against a case where the forward flux is zero or negative.
if (forward_flux[sr * negroups_ + g] > 0.0) {
source_regions_.external_source(sr, g) =
1.0 / forward_flux[sr * negroups_ + g];
} else {
source_regions_.external_source(sr, g) = 0.0;
}
source_regions_.external_source(sr, g) =
1.0 / forward_flux[sr * negroups_ + g];
}
}

Expand Down