Merge pull request #5381 from ye-luo/fix-spline-dump-read

Spline dump should only be read by rank 0.

Author: prckent

src/QMCWaveFunctions/BsplineFactory/HybridRepSetReader.cpp

@@ -148,7 +148,7 @@ std::unique_ptr<SPOSet> HybridRepSetReader<SA>::create_spline_set(const std::str
   bool foundspline = spline_reader_.createSplineDataSpaceLookforDumpFile(bandgroup, *bspline);
   typename SA::HYBRIDBASE& hybrid_center_orbs = *bspline;
   hybrid_center_orbs.resizeStorage(bspline->myV.size());
-  if (foundspline)
+  if (foundspline && myComm->rank() == 0)
   {
     Timer now;
     hdf_archive h5f(myComm);

src/QMCWaveFunctions/BsplineFactory/SplineSetReader.cpp

@@ -42,7 +42,7 @@ std::unique_ptr<SPOSet> SplineSetReader<SA>::create_spline_set(const std::string
   auto bspline = std::make_unique<SA>(my_name, use_offload);
   app_log() << "  ClassName = " << bspline->getClassName() << std::endl;
   bool foundspline = createSplineDataSpaceLookforDumpFile(bandgroup, *bspline);
-  if (foundspline)
+  if (foundspline && myComm->rank() == 0)
   {
     Timer now;
     hdf_archive h5f(myComm);
@@ -80,6 +80,7 @@ std::unique_ptr<SPOSet> SplineSetReader<SA>::create_spline_set(const std::string
   }
 
   {
+    myComm->barrier();
     Timer now;
     bspline->bcast_tables(myComm);
     app_log() << "  Time to bcast the table = " << now.elapsed() << std::endl;