Function WriteFab

Tools/Postprocessing/C_Src/MultiFabToMatLab.cpp:37–69

}

static
void
WriteFab (std::ostream&         os,
          const FArrayBox& fab,
          const char*      name)

{
    int nx = fab.box().length(0);
    int ny = fab.box().length(1);
    int dim = BL_SPACEDIM;


#if (BL_SPACEDIM == 2)
    for (int i = 0; i < nx; i++)
     for (int j = 0; j < ny; j++)
     {
       int index = j*nx + i;
       const Real * ptr = fab.dataPtr();
       os.write((char*)(ptr+index),sizeof(Real));
     }

#elif (BL_SPACEDIM == 3)
    int nz = fab.box().length(2);
    for (int i = 0; i < nx; i++)
     for (int j = 0; j < ny; j++)
      for (int k = 0; k < nz; k++)
      {
        int index = k*(nx*ny) + j*nx + i;
        const Real * ptr = fab.dataPtr();
        os.write((char*)(ptr+index),sizeof(Real));
      }
#endif
}


int main(int argc, char* argv[])