[ET Trac] [Einstein Toolkit] #625: CapretIOHDF5 sliced output does not output symmetry points unless it also outputs buffer points

[Einstein Toolkit]

#625: CapretIOHDF5 sliced output does not output symmetry points unless it also
outputs buffer points
--------------------------+-------------------------------------------------
 Reporter:  hinder        |       Owner:  eschnett
     Type:  defect        |      Status:  new     
 Priority:  major         |   Milestone:          
Component:  Carpet        |     Version:          
 Keywords:  CarpetIOHDF5  |  
--------------------------+-------------------------------------------------
 The parameter CarpetIOHDF5::output_symmetry_points = yes does not work
 when CarpetIOHDF5::output_buffer_points = no (the default).  This is
 probably due to the line

     if (not output_buffer_points) {
       exts &= allactive;
     }

 near line 1071 in OutputSlice.cc.  I suspect that allactive does not
 include the symmetry points.  The attached parameter file should be run on
 one processor.  You can then use

   h5dump -d '/GRID::x it=0 tl=0 rl=0' output_symmetry/x.x.h5

 Near the top you will see

   DATASET "/GRID::x it=0 tl=0 rl=0" {
      DATATYPE  H5T_IEEE_F64LE
      DATASPACE  SIMPLE { ( 5 ) / ( 5 ) }
      DATA {
      (0): 0, 0.1, 0.2, 0.3, 0.4
      }

 i.e. the first point output is x = 0.  For this parameter file, there
 should be a symmetry point at x = -0.1.  If you modify the parameter file
 to use CarpetIOHDF5::output_buffer_points = yes, then the symmetry point
 appears in the output.

 The workaround is to always use

   CarpetIOHDF5::output_buffer_points = yes

 if you want to output symmetry points.

-- 
Ticket URL: <https://trac.einsteintoolkit.org/ticket/625>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit