Add failing test from bug report

Author: franzpoeschel

CMakeLists.txt

@@ -791,6 +791,7 @@ if(openPMD_BUILD_TESTING)
             list(APPEND ${out_list}
                 test/Files_ParallelIO/read_variablebased_randomaccess.cpp
                 test/Files_ParallelIO/iterate_nonstreaming_series.cpp
+                test/Files_ParallelIO/bug_1655_bp5_writer_hangup.cpp
             )
         elseif(${test_name} STREQUAL "Core")
             list(APPEND ${out_list}

test/Files_ParallelIO/ParallelIOTests.hpp

@@ -75,4 +75,9 @@ namespace iterate_nonstreaming_series
 auto iterate_nonstreaming_series() -> void;
 }
 
+namespace bug_1655_bp5_writer_hangup
+{
+auto bug_1655_bp5_writer_hangup() -> void;
+}
+
 #endif

test/Files_ParallelIO/bug_1655_bp5_writer_hangup.cpp

@@ -0,0 +1,46 @@
+#include "ParallelIOTests.hpp"
+
+#include <mpi.h>
+#include <openPMD/openPMD.hpp>
+
+namespace bug_1655_bp5_writer_hangup
+{
+auto worker(std::string const &ext) -> void
+{
+    int mpi_size;
+    int mpi_rank;
+
+    MPI_Comm_size(MPI_COMM_WORLD, &mpi_size);
+    MPI_Comm_rank(MPI_COMM_WORLD, &mpi_rank);
+
+    auto const value = float(mpi_size * 100 + mpi_rank);
+    std::vector<float> local_data(10 * 300, value);
+
+    std::string filename = "../samples/ptl_%T." + ext;
+
+    Series series = Series(filename, Access::CREATE, MPI_COMM_WORLD);
+
+    Datatype datatype = determineDatatype<float>();
+
+    auto myptl = series.writeIterations()[1].particles["ion"];
+    Extent global_ptl = {10ul * mpi_size * 300};
+    Dataset dataset_ptl = Dataset(datatype, global_ptl, "{}");
+    myptl["charge"].resetDataset(dataset_ptl);
+
+    series.flush();
+
+    if (mpi_rank == 0) // only rank 0 adds data
+        myptl["charge"].storeChunk(local_data, {0}, {3000});
+
+    series.flush(); // hangs here
+    series.close();
+}
+
+auto bug_1655_bp5_writer_hangup() -> void
+{
+    for (auto const &ext : testedFileExtensions())
+    {
+        worker(ext);
+    }
+}
+} // namespace bug_1655_bp5_writer_hangup

test/ParallelIOTest.cpp

@@ -2183,3 +2183,10 @@ TEST_CASE("iterate_nonstreaming_series", "[serial][adios2]")
 }
 
 #endif // openPMD_HAVE_ADIOS2 && openPMD_HAVE_MPI
+
+#if openPMD_HAVE_MPI
+TEST_CASE("bug_1655_bp5_writer_hangup", "[parallel]")
+{
+    bug_1655_bp5_writer_hangup::bug_1655_bp5_writer_hangup();
+}
+#endif