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

avoid segfault with idiomatic thresholding, rays through pixel centers

Browse source
This commit is contained in:
Gavin Ridley 2023-04-05 11:50:45 -04:00
parent 5c4144d520
commit a92f2d1f9f
2 changed files with 8 additions and 6 deletions
Download patch file Download diff file Expand all files Collapse all files

12
src/plot.cpp
View file

@ -1273,8 +1273,10 @@ void ProjectionPlot::create_output() const
// Generate the starting position/direction of the ray
if (orthographic_width_ == 0.0) { // perspective projection
double this_phi = -horiz_fov_radians / 2.0 + dphi * horiz;
double this_mu = -vert_fov_radians / 2.0 + dmu * vert + M_PI / 2.0;
double this_phi =
-horiz_fov_radians / 2.0 + dphi * horiz + 0.5 * dphi;
double this_mu =
-vert_fov_radians / 2.0 + dmu * vert + M_PI / 2.0 + 0.5 * dmu;
Direction camera_local_vec;
camera_local_vec.x = std::cos(this_phi) * std::sin(this_mu);
camera_local_vec.y = std::sin(this_phi) * std::sin(this_mu);
@ -1436,9 +1438,9 @@ void ProjectionPlot::create_output() const
if (i * i + j * j < wireframe_thickness_ * wireframe_thickness_) {
// Check if wireframe pixel is out of bounds
int w_i = std::max(std::min(horiz + i, pixels_[0]), 0);
int w_j = std::max(std::min(vert + j, pixels_[1]), 0);
data(i, j) = wireframe_color_;
int w_i = std::max(std::min(horiz + i, pixels_[0] - 1), 0);
int w_j = std::max(std::min(vert + j, pixels_[1] - 1), 0);
data(w_i, w_j) = wireframe_color_;
}
}
}