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

centre to centre_

Browse source
This commit is contained in:
RemDelaporteMathurin 2022-10-07 08:40:04 +00:00
parent 21200d9b6e
commit 896ae2792e
1 changed files with 8 additions and 8 deletions
Download patch file Download diff file Expand all files Collapse all files

16
src/mesh.cpp
View file

@ -955,7 +955,7 @@ CylindricalMesh::CylindricalMesh(pugi::xml_node node) : StructuredMesh {node}
grid_[0] = get_node_array<double>(node, "r_grid");
grid_[1] = get_node_array<double>(node, "phi_grid");
grid_[2] = get_node_array<double>(node, "z_grid");
centre = get_node_array<double>(node, "centre");
centre_ = get_node_array<double>(node, "centre");
if (int err = set_grid()) {
fatal_error(openmc_err_msg);
@ -985,9 +985,9 @@ StructuredMesh::MeshIndex CylindricalMesh::get_indices(
}
// TODO: pass centre as argument
mapped_r[0] += centre[0]
mapped_r[1] += centre[1]
mapped_r[2] += centre[2]
mapped_r[0] += centre_[0]
mapped_r[1] += centre_[1]
mapped_r[2] += centre_[2]
MeshIndex idx = StructuredMesh::get_indices(mapped_r, in_mesh);
@ -1193,7 +1193,7 @@ SphericalMesh::SphericalMesh(pugi::xml_node node) : StructuredMesh {node}
grid_[0] = get_node_array<double>(node, "r_grid");
grid_[1] = get_node_array<double>(node, "theta_grid");
grid_[2] = get_node_array<double>(node, "phi_grid");
centre = get_node_array<double>(node, "centre");
centre_ = get_node_array<double>(node, "centre");
if (int err = set_grid()) {
fatal_error(openmc_err_msg);
@ -1224,9 +1224,9 @@ StructuredMesh::MeshIndex SphericalMesh::get_indices(
}
// TODO: pass centre as argument
mapped_r[0] += centre[0]
mapped_r[1] += centre[1]
mapped_r[2] += centre[2]
mapped_r[0] += centre_[0]
mapped_r[1] += centre_[1]
mapped_r[2] += centre_[2]
MeshIndex idx = StructuredMesh::get_indices(mapped_r, in_mesh);