pc2compare.C

Go to the documentation of this file.

#include "PC2Util.H"

int 
main(int argc,char *argv[])
{
  double errorTolerance = 1e-15;
  if(argc == 1){
    std::cout << "This utility compares the state data in two specified"          << std::endl
              << "PlasCom2 or PlasComCM HDF5 files. Cross-comparisons between"    << std::endl
              << "PlasCom2 and PlasComCM are supported. A plottable summary"      << std::endl
              << "of the comparison is generated on stdout. Return values are:"   << std::endl
              << "  0  : All tests passed."         << std::endl
              << "  1  : One or more tests failed." << std::endl
              << " -1  : Tests not run due to some error or incompatibility."     << std::endl
              << std::endl << "Some limitations:" << std::endl
              << "* Currently, only rho, rhoV, and rhoE fields are compared."       << std::endl
              << "* Comparisons on subsets of domains are not yet supported."       << std::endl
              << "* Only single grid datasets are supported (default to 1st grid)." << std::endl
              << std::endl
              << " Usage: " << std::endl
              << argv[0] << " <red file> <blue file> [tolerance]" << std::endl
              << std::endl
              << "<red/blue file> := Required PlasCom2 or PlasComCM HDF5 files." << std::endl
              << "tolerance       := Optional error tolerance defaults to 1e-15." << std::endl;
    return(-1);
  }

  if(argc < 3){
    std::cerr << argv[0] 
              << ":Error: Missing some required arguments. Invoke without arguments for usage summary."
              << std::endl;
    return(-1);
  }
  
  if(argc > 3){
    std::istringstream Istr(argv[3]);
    Istr >> errorTolerance;
  }

  return(plascom2::util::PC2Compare(argv[1],argv[2],errorTolerance,std::cout));

}