#446: 1D HDF5 output nonfunctional with multipatch - std::out_of_range ---------------------+------------------------------------------------------ Reporter: hinder | Owner: eschnett Type: defect | Status: new Priority: major | Milestone: Component: Carpet | Version: Resolution: | Keywords: ---------------------+------------------------------------------------------
Comment (by rhaas):
Hello Ian,
The OutputSlice.cc:564 line is the GetAllActive call here:
// re-compute the active (non-buffered) region ibset allactive; GetAllActive (dd, hh, m, rl, allactive);
Hmm, so this happens even with my initial patch it seems. GetAllActive's result is really only needed when you set output_buffer_points (or whatever it is actually called) to "no". This happens only for MultiPatch runs (and only for the "extra" non-refined patches)? If you need a quick fix so that you can actually run something (and don't need output_buffer_zones) then the attached patch disables GetAllActive when output_buffer_points="no" (which is the default). I put in the GetAllActive function (it's more or less a copy of a piece of dh::regrid). But I cannot reproduce your error with grid functions (ie. all of out[0123]D_var = "grid::x") works fine for me). Grid array output also seems to work at least for out[012]D_var = "sphericalsurface::sf_radius". This really seems to be a multipatch only issue.