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

Updating how the libmesh solution set is handled. Also applying data for a single score, single nuclide tally.

Browse source
This commit is contained in:
Patrick Shriwise 2019-12-18 17:16:20 -06:00
parent 2226dc9f6e
commit 2410a4dfba
2 changed files with 4 additions and 9 deletions
Download patch file Download diff file Expand all files Collapse all files

1
include/openmc/mesh.h
View file

@ -521,7 +521,6 @@ private:
std::unique_ptr<libMesh::Mesh> m_;
std::unique_ptr<libMesh::PointLocatorBase> point_locator_;
std::unique_ptr<libMesh::EquationSystems> equation_systems_;
std::unique_ptr<libMesh::NumericVector<libMesh::Real>> sol_;
std::map<std::string, unsigned int> variable_map_;
std::string eq_system_name_;
libMesh::Elem* first_element_;

12
src/mesh.cpp
View file

@ -2064,8 +2064,6 @@ LibMesh::LibMesh(pugi::xml_node node) : UnstructuredMeshBase(node) {
libMesh::ExplicitSystem& eq_sys =
equation_systems_->add_system<libMesh::ExplicitSystem>(eq_system_name_);
sol_ = eq_sys.current_local_solution->clone();
point_locator_ = m_->sub_point_locator();
point_locator_->enable_out_of_mesh_mode();
@ -2154,17 +2152,15 @@ LibMesh::set_variable(const std::string& var_name, int bin, double value) {
std::vector<libMesh::dof_id_type> dof_indices;
dof_map.dof_indices(e, dof_indices, var_num);
Ensures(dof_indices.size() == 1);
(*sol_).set(dof_indices[0], value);
(*eqn_sys.solution).set(dof_indices[0], value);
}
void LibMesh::write() const {
// update solution
auto& eqn_sys = equation_systems_->get_system(eq_system_name_);
*eqn_sys.solution = *sol_;
libMesh::ExodusII_IO exo(*m_);
std::set<std::string> systems_out = {eq_system_name_};
exo.write_discontinuous_exodusII("test.exo", *equation_systems_, &systems_out);
std::stringstream output_filename;
output_filename << "umesh_" << id_ << ".exo";
exo.write_discontinuous_exodusII(output_filename.str(), *equation_systems_, &systems_out);
}