PC2Initialize.C

Go to the documentation of this file.

#include "PlasCom2.H"
#include "PCPPUtil.H"
void pcpp::commandline::Initialize()
{
  InitializeDefaults();
  AddOption('c',"config",2,"configfile");
  AddHelp('c',"Set the configuration file.");
  AddArgument("restartFile");
  AddArgHelp("restartFile","HDF5 restart file.");
  AddArgument("targetFile"); 
  AddArgHelp("targetFile","HDF5 file with target data.");
  SetDescription(std::string("PlasCom2 alpha development."));
  SetNotes("Notes for PlasCom2.\n");
}

namespace plascom2 {

  std::string Splash()
  {
    std::ostringstream outStream;
    outStream << "*******************************" << std::endl
              << "*        XPACC PlasCom2       *" << std::endl
              << "*       Alpha development     *" << std::endl
              << "*******************************" << std::endl;
    return(outStream.str());
  };
  

  int application::InitializeApplication()
  {
    int error = pcpp::parallelprogram::Initialize();
    if(error < 0){
      ErrOut("Parallel program failed to initialize.\n");
      return(1);
    } else if (error > 0){
      return(1);
    }
    error = 0;
    StdOut(Splash(),1);
    // Configuration file processing
    globalCommunicator.Barrier();
    return(error);
  }

  //  int application::ConfigureApplication()
  //  {
    // if(InitializeNodeCommunicator())
    //   return(1);
    // if(InitializeIOCommunicators())
    //   return(1);
    // if(SetupIOTopology())
    //   return(1);
    // if(PartitionGridsLogical())
    //   return(1);
    // // if(PrepareHalos())
    // //   return(1);
    // std::ostringstream Ostr;
    // Ostr << "Plascom domain configuration" << std::endl
    //      << "----------------------------" << std::endl;
    // plascomDomain.Report(Ostr);
    // Ostr << "----------------------------" << std::endl;
    // StdOut(Ostr.str(),2);
    // if(InitializeGrids())
    //   return(1);
    // if(DistributeGridsLogical())
    //   return(1);
    // if(GenerateGrids())
    //   return(1);
    // // Can be moved around (fixme: production)
    // if(DistributeGrids())
    //   return(1);
    // if(ExchangeIONeighborhood())
    //   return(1);
    // if(ExchangeComputeNeighborhood())
    //   return(1);
    // if(CollideComputeNeighbors())
    //   return(1);
    // if(ExchangeCommunicationMaps())
    //   return(1);
    // if(DistributeCommunicationMaps())
    //   return(1);
    // if(CreateFieldDataBuffers())
    //   return(1);
    // if(CreateCommunicationBuffers())
    //   return(1);
    // // Not needed - included here as a test (fixme: production)
    // //  if(GatherGrids())
    // //    return(1);
    // globalCommunicator.Barrier();
  //    FunctionExit("Setup");
  //    return(0);
  //  }
  
  int application::FinalizeApplication()
  {
#ifdef USE_OVERKIT
    FinalizeOverkit();
#endif
    int numProcs = NProc();
//     int numSteps = 0;
//     size_t numPoints = 0;
//     int numGrids = 0;
    //    maxSteps = plascomConfig.GetValue<int>("Stepping:NumSteps");
    // if(!plascomDomain.empty()){
    //   numGrids = 1;
    //   pcpp::geometry::grid &a1Grid(plascomDomain[0]);
    //   pcpp::geometry::gridinfo &gridInfo(a1Grid.Info());
    //   std::vector<int> &globalSize(gridInfo.globalSize);
    //   std::vector<int>::iterator sizeIt = globalSize.begin();
    //   numPoints = 1;
    //   while(sizeIt != globalSize.end())
    //         numPoints *= *sizeIt++;
    // }
    std::ostringstream timerReportStream;
    timerReportStream << "#Number of procs:     " 
                      << numProcs << std::endl
                      << "#Number of threads:   " 
                      << numThreads << std::endl
                      << "#Number of steps:     " 
                      << numSteps << std::endl
                      << "#Number of grids:     " 
                      << numGrids << std::endl;
    if(numGrids > 0){
      timerReportStream << "#Number of points:    ";
      for(int iGrid = 0;iGrid < numGrids;iGrid++)
        timerReportStream << numPointsGrid[iGrid] << " ";
      timerReportStream << std::endl;
    }

    Report(timerReportStream);
    std::ostringstream timingFileNameStream;
    timingFileNameStream << "PlasCom2Timing_";
    int nZeros = (numProcs > 99999 ? 0  :
                  numProcs > 9999  ? 1  :
                  numProcs > 999   ? 2  :
                  numProcs > 99    ? 3  :
                  numProcs > 9     ? 4  : 5);
    int zeroCount = 0;
    while(zeroCount++ < nZeros)
      timingFileNameStream << "0";
    timingFileNameStream << numProcs << ".txt";
    std::ofstream timingFile;
    if(Rank() == 0){
      timingFile.open(timingFileNameStream.str().c_str());
      timingFile << timerReportStream.str() << std::endl;
      timingFile.close();
    } else if(OutStreamReady()) {
      OutStream() << timerReportStream.str() << std::endl;
    }
    if(pcpp::ParallelProgramType::Finalize()){
      ErrOut("Failed to finalize parallel program.\n");
      if(ErrStreamReady())
        DumpErrors(ErrStream());
      return(1);
    }
    StdOut("Program finished - goodbye.\n",1);
    return(0);
  }
}