Function main

Tools/Postprocessing/C_Src/MultiFabToMatLab.cpp:72–132

int main(int argc, char* argv[])
{
  amrex::Initialize(argc, argv);
  {
    if(argc == 1)
      PrintUsage(argv[0]);

    if (ParallelDescriptor::NProcs() > 1){
      amrex::Error("This is an inherently serial program!");
    }
    ParmParse pp;

    std::string name;
    pp.get("infile", name);

    //
    // MatLab expects native floating-point format.
    //
    FArrayBox::setFormat(FABio::FAB_NATIVE);

    VisMF::IO_Buffer io_buffer(VisMF::IO_Buffer_Size);
    char buf[128];
    std::string file = name;
    file += ".mat";
    std::ofstream os;

    os.rdbuf()->pubsetbuf(io_buffer.dataPtr(), io_buffer.size());

    os.open(file.c_str(), std::ios::out|std::ios::binary);

    if(os.fail()) {
      amrex::FileOpenFailed(file);
    }


    MultiFab in;
    VisMF::Read(in, name);
    const BoxArray ba = in.boxArray();

    //Fake a xlo to xhi
    for(int i = 0; i < ba.size(); ++i)
    {
      const Box& b = ba[i];
      Real xlo[BL_SPACEDIM], xhi[BL_SPACEDIM];
      for(int d = 0; d <BL_SPACEDIM; ++d)
      {
        xlo[d] = b.loVect()[d];
        xhi[d] = b.hiVect()[d];
        os.write((char*)&xlo,sizeof(Real));
        os.write((char*)&xhi,sizeof(Real));
      }
    }
    //Write the Fab Data
    for(int i = 0; i < ba.size(); ++i)
    {
      WriteFab(os, in[i], buf);
    }
    os.close();