Method extendOffsetDataset

SRC/recorder/VTKHDF_Recorder.cpp:3564–3651

}

int VTKHDF_Recorder::extendOffsetDataset(hid_t group_id, 
                  const char* dataset_name, 
                  const void* new_value, 
                  hid_t datatype,
                  size_t extra_rows) {
    // Open dataset
    hid_t dset_id = H5Dopen(group_id, dataset_name, H5P_DEFAULT);
    if (dset_id < 0) {
        opserr << "Error opening dataset '" << dataset_name << "'\n";
        return -1;
    }

    // Get current dataspace and dimensions
    hid_t space_id = H5Dget_space(dset_id);
    if (space_id < 0) {
        opserr << "Error getting dataspace for '" << dataset_name << "'\n";
        H5Dclose(dset_id);
        return -1;
    }

    // Get current dimensions
    hsize_t current_dims[1];
    if (H5Sget_simple_extent_dims(space_id, current_dims, NULL) < 0) {
        opserr << "Error getting current dimensions\n";
        H5Sclose(space_id);
        H5Dclose(dset_id);
        return -1;
    }

    // Calculate new dimensions
    hsize_t new_dims[1] = {current_dims[0] + static_cast<hsize_t>(extra_rows)};

    // Extend dataset
    if (H5Dset_extent(dset_id, new_dims) < 0) {
        opserr << "Error setting new extent for '" << dataset_name << "'\n";
        H5Sclose(space_id);
        H5Dclose(dset_id);
        return -1;
    }

    // Get new dataspace
    hid_t new_space_id = H5Dget_space(dset_id);
    if (new_space_id < 0) {
        opserr << "Error getting new dataspace for '" << dataset_name << "'\n";
        H5Sclose(space_id);
        H5Dclose(dset_id);
        return -1;
    }

    // Select hyperslab for new data
    hsize_t offset[1] = {current_dims[0]};
    hsize_t count[1] = {1};
    if (H5Sselect_hyperslab(new_space_id, H5S_SELECT_SET, offset, NULL, count, NULL) < 0) {
        opserr << "Error selecting hyperslab\n";
        H5Sclose(new_space_id);
        H5Sclose(space_id);
        H5Dclose(dset_id);
        return -1;