PlasCom2IO.C

Go to the documentation of this file.

#include "PlasCom2.H"
#include "PCPPCommUtil.H"
#include "PCPPIO.H"
#include "PCPPReport.H"
#include "PCPPIntervalUtils.H"

namespace plascom2 {
  namespace io {
    namespace hdf5 {

      template<typename GridType,typename StateType,typename ConfigType>
      int OutputSingle(const std::string &fileName,const GridType &inGrid,const StateType &inState,
                       const ConfigType &simConfig,const ConfigType &gridConfig,const ConfigType &stateConfig) {

        pcpp::CommunicatorType &gridComm(inGrid.Communicator());


        // Begin PlasCom2 output
        if(ix::sys::FILEEXISTS(fileName))
          ix::sys::Remove(fileName);
        pcpp::io::hdf5::base pc2File(fileName,gridComm);
        pc2File.CreateGroup("/PlasCom2");
        pc2File.CreateGroup("/PlasCom2/Configuration");
        std::ostringstream configStream;
        configStream << simConfig;
        std::string configString(configStream.str());
        if(!configString.empty())
          pc2File.CreateAttribute("/PlasCom2/Configuration/Parameters",configString);

        pc2File.CreateGroup("/PlasCom2/Geometry");
        std::vector<size_t> attSize(1,1);
        pc2File.CreateAttribute("/PlasCom2/Geometry/numGeometries",attSize,&numGeometries);
        pc2File.CreateAttribute("/PlasCom2/Geometry/geometryNames",geometryName);
        std::string geomGroupName(std::string("/PlasCom2/Geometry/")+geometryName);
        pc2File.CreateGroup(geomGroupName);
        int numGrids = 1;
        pc2File.CreateAttribute(geomGroupName+"/numGrids",attSize,&numGrids);
        pc2File.CreateAttribute(geomGroupName+"/gridNames",gridName);
        pcpp::io::hdf5::WriteGrid(testGrid,gridName,geomGroupName,pc2File);
        
        
        pc2File.CreateGroup("/PlasCom2/Simulation");
        pc2File.CreateAttribute("/PlasCom2/Simulation/numDomains",attSize,&numDomains);
        pc2File.CreateAttribute("/PlasCom2/Simulation/domainNames",domainName);
        pcpp::io::hdf5::WriteStateAttributes(simState,'s',"/PlasCom2/Simulation",pc2File);
        pcpp::io::hdf5::WriteStateAttributes(paramState,'s',"/PlasCom2/Simulation",pc2File);
        std::string domainGroupName(std::string("/PlasCom2/Simulation/")+domainName);
        pc2File.CreateGroup(domainGroupName);
        pc2File.CreateGroup(domainGroupName+"/Configuration");
        std::ostringstream domainConfigStream;
        domainConfigStream << domainConfig;
        pc2File.CreateAttribute(domainGroupName+"/Configuration/Parameters",domainConfigStream.str());
        pc2File.CreateAttribute(domainGroupName+"/numGrids",attSize,&numDomainGrids);
        pc2File.CreateAttribute(domainGroupName+"/gridNames",domainGridName);
        pcpp::io::hdf5::WriteStateAttributes(simState,'d',domainGroupName,pc2File);
        pcpp::io::hdf5::WriteStateAttributes(paramState,'d',domainGroupName,pc2File);
        std::string gridGroupName(domainGroupName+"/"+domainGridName);
        pc2File.CreateGroup(gridGroupName);
        std::vector<size_t> partitionSizes(partInterval.Sizes());
        std::vector<size_t> partitionStarts(partInterval.Starts());
        std::vector<size_t> partitionBufferStarts(partitionBufferInterval.Starts());
        pcpp::io::hdf5::WriteStateNodeData(simState,gridSizes,partitionSizes,partitionStarts,
                                           bufferSizes,partitionBufferStarts,gridGroupName,pc2File);  
        pc2File.Close();     
      }
    }
    
  }
}