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

Populating c-side structures. Updating Fortran universe population in read_geometry_xml.

Browse source
This commit is contained in:
shriwise 2018-05-16 17:24:52 -05:00 • committed by pshriwise
parent 8e5f424ae1
commit df342211cc
2 changed files with 22 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

27
src/cad.cpp
View file

@ -15,17 +15,32 @@ void load_cad_geometry_c()
rval = DAGMC->init_OBBTree();
MB_CHK_ERR_CONT(rval);
int32_t cad_univ_id = 0; // universe is always 0 for CAD
// initialize cell objects
openmc::n_cells = DAGMC->num_entities(3);
for(int i = 0; i < openmc::n_cells; i++)
{
// set cell ids using global IDs
openmc::CADCell c = openmc::CADCell();
c.id = DAGMC->id_by_index(3, i);
c.dagmc_ptr = DAGMC;
c.universe = 0; // set to zero for now
openmc::cells_c.push_back(&c);
openmc::cell_dict[c.id] = i;
openmc::CADCell* c = new openmc::CADCell();
c->id = DAGMC->id_by_index(3, i);
c->dagmc_ptr = DAGMC;
c->universe = cad_univ_id; // set to zero for now
c->material.push_back(40); // TEMPORARY
openmc::cells_c.push_back(c);
openmc::cell_dict[c->id] = i;
// Populate the Universe vector and dict
auto it = openmc::universe_dict.find(cad_univ_id);
if (it == openmc::universe_dict.end()) {
openmc::universes_c.push_back(new openmc::Universe());
openmc::universes_c.back()-> id = cad_univ_id;
openmc::universes_c.back()->cells.push_back(i);
openmc::universe_dict[cad_univ_id] = openmc::universes_c.size() - 1;
}
else {
openmc::universes_c[it->second]->cells.push_back(i);
}
}
// initialize surface objects

2
src/input_xml.F90
View file

@ -435,7 +435,7 @@ contains
u % cells(:) = 0
end associate
end do
root_universe = find_root_universe() + 1
root_universe = 0 + 1
do i = 1, n_cells
! Get index in universes array