PlasCom2.H

Go to the documentation of this file.

#ifndef __PC2_H__
#define __PC2_H__

#include "Stencil.H"
#include "Simulation.H"

#ifdef USE_OVERKIT
#include "overkit.h"
#endif

namespace plascom2 {
  
  typedef simulation::grid::parallel_blockstructured                grid_t;
  typedef simulation::state::base                                   state_t;
  typedef simulation::grid::halo                                    halo_t;
  typedef simulation::domain::base<grid_t,state_t>                  domain_t;
  typedef simulation::application<domain_t>                         application_t;
  
  // struct domain_info {
  //   int numDomains;
  //   std::vector<int> globalIndex;
  // };
  
  struct process_info {
    int worldRank;
    int ioRank;
    int computeRank;
  };

#ifdef USE_OVERKIT
  struct overkit_data {
    ovk_context *context;
    std::vector<ovk_domain *> domains;
  };
#endif

  class application : public application_t
  {
    

  public:

    application(int numArgs,char **argList,MPI_Comm mpiCommunicator) :
      application_t(numArgs,argList,mpiCommunicator), isRestart(0),
      numSteps(0), numGrids(0)
    {};
    virtual int ConfigureApplication();
    virtual int ConfigureSimulation();
    virtual int ConfigureDomains();
    virtual int InitializeApplication();
    virtual int RunApplication();
    virtual int FinalizeApplication();
    virtual int WriteDomains();
    virtual int ReadDomains();
    virtual int Restart() { return(isRestart); };
    void SetRestart(int trufals = 1){ isRestart = trufals; };
    int RunEulerPBSPeriodic();
    int RunPlasCom2();
    int RunMaxwellSolver();
    int ConfigureGeometries();
    void SetNumGrids(int nGrids)
    {
      numGrids = nGrids;
      numPointsGrid.resize(numGrids);
    };
    int SetNumPointsGrid(int iGrid,size_t numPoints)
    {
      if(iGrid > (numGrids - 1))
        return(1);
      numPointsGrid[iGrid] = numPoints;
      return(0);
    };
    void SetNumSteps(size_t nSteps) {numSteps = nSteps;};
#ifdef USE_OVERKIT
    int InitializeOverkit();
    int ConfigureOverkit();
    int FinalizeOverkit();
    overkit_data *GetOverkitData() { return &overkitData; }
#endif

  protected:
    pcpp::PlatformInfoType platformInfo;
    pcpp::VirtualNodeType  virtualNode;
    int isRestart;
#ifdef USE_OVERKIT
    overkit_data overkitData;
#endif
    size_t numSteps;
    int numGrids;
    std::vector<size_t> numPointsGrid;
  };
}
  
#endif