diff --git a/src/mesh.cpp b/src/mesh.cpp index 30553c8c74..333a5dc331 100644 --- a/src/mesh.cpp +++ b/src/mesh.cpp @@ -60,7 +60,7 @@ std::vector> meshes; #ifdef LIBMESH namespace settings { std::unique_ptr libmesh_init; -const libMesh::Parallel::Communicator* libmesh_comm; +const libMesh::Parallel::Communicator* libmesh_comm {nullptr}; } #endif @@ -2155,6 +2155,9 @@ void LibMesh::initialize() fatal_error("Attempting to use an unstructured mesh without a libMesh communicator."); } + // assuming that unstructured meshes used in OpenMC are 3D + n_dimension_ = 3; + m_ = std::make_unique(*settings::libmesh_comm, n_dimension_); m_->read(filename_); m_->prepare_for_use();