#include <iostream>
#include "moab/Core.hpp"
#include <moab/ScdInterface.hpp>
#include <moab/ParallelComm.hpp>

int main(int argc, char **argv) {
  MPI_Init(&argc, &argv);

  moab::Core* mymbcore = new moab::Core();
  moab::Interface* mymb = dynamic_cast<moab::Interface*>(mymbcore);

  moab::ScdInterface* myscdiface;
  auto rval = mymb->query_interface(myscdiface);

  moab::EntityHandle rootset,partnset;
  rval = mymb->create_meshset(moab::MESHSET_SET, rootset);
  rval = mymb->create_meshset(moab::MESHSET_SET, partnset);

  MPI_Comm myComm = MPI_COMM_WORLD;
  const MPI_Comm* comm = const_cast<const MPI_Comm*> (&myComm);
  moab::ParallelComm* mypc = moab::ParallelComm::get_pcomm(mymb,partnset,comm);

  moab::ScdParData* mypar_data = new moab::ScdParData();
  mypar_data->pComm = mypc;

  mypar_data->gDims[0] = 0;
  mypar_data->gDims[1] = 0;
  mypar_data->gDims[2] = 0;
  mypar_data->gDims[3] = 64;
  mypar_data->gDims[4] = 16;
  mypar_data->gDims[5] = 0;

  mypar_data->partMethod = moab::ScdParData::SQIJ;

  moab::ScdBox* mybox;
  rval = myscdiface->construct_box(
    moab::HomCoord(), moab::HomCoord(), NULL, 0,
    mybox, NULL, mypar_data);
   
  MPI_Finalize();
  return 0;
}