OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

HDF5File: add recursive group creation in write

Browse Source
This commit is contained in:
Arne Morten Kvarving 2023-02-08 12:21:52 +01:00
parent e5c2fcfaf2
commit a90189c78d
1 changed files with 20 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
opm/simulators/utils/HDF5File.cpp
View File

@ -22,6 +22,8 @@
#include <config.h> #include <config.h>
#include <opm/simulators/utils/HDF5File.hpp> #include <opm/simulators/utils/HDF5File.hpp>
#include <opm/common/utility/String.hpp>
#include <filesystem> #include <filesystem>
#include <stdexcept> #include <stdexcept>
@ -80,9 +82,26 @@ void HDF5File::write(const std::string& group,
if (groupExists(m_file, group)) { if (groupExists(m_file, group)) {
grp = H5Gopen2(m_file, group.c_str(), H5P_DEFAULT); grp = H5Gopen2(m_file, group.c_str(), H5P_DEFAULT);
} else { } else {
auto grps = split_string(group, '/');
std::string curr;
for (size_t i = 0; i < grps.size()-1; ++i) {
curr += '/';
curr += grps[i];
if (!groupExists(m_file, curr)) {
hid_t subgrp = H5Gcreate2(m_file, curr.c_str(), 0, H5P_DEFAULT, H5P_DEFAULT);
if (subgrp == H5I_INVALID_HID) {
throw std::runtime_error("HDF5File: Failed to create group '" + curr + "'");
}
H5Gclose(subgrp);
}
}
grp = H5Gcreate2(m_file, group.c_str(), 0, H5P_DEFAULT, H5P_DEFAULT); grp = H5Gcreate2(m_file, group.c_str(), 0, H5P_DEFAULT, H5P_DEFAULT);
} }
if (grp == H5I_INVALID_HID) {
throw std::runtime_error("HDF5File: Failed to create group '" + group + "'");
}
hsize_t size = buffer.size(); hsize_t size = buffer.size();
hsize_t start = 0; hsize_t start = 0;
@ -90,7 +109,7 @@ void HDF5File::write(const std::string& group,
hid_t dataset_id = H5Dcreate2(grp, dset.c_str(), H5T_NATIVE_CHAR, space, hid_t dataset_id = H5Dcreate2(grp, dset.c_str(), H5T_NATIVE_CHAR, space,
H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT); H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
if (dataset_id == H5I_INVALID_HID) { if (dataset_id == H5I_INVALID_HID) {
throw std::runtime_error("HDF5File: Trying to write already existing dataset " + group + '/' + dset); throw std::runtime_error("HDF5File: Trying to write already existing dataset '" + group + '/' + dset + "'");
} }
if (size > 0) { if (size > 0) {