main.cpp

Go to the documentation of this file.

#include "mainapp.h"
#include <stdlib.h>
#include <signal.h>
#include <fenv.h>
#include <limits.h>
#include "hdf5orthowriter.h"
#include "flashco2brine.h"
#ifdef USE_MPI
#include"netmpi.h"
#include "mainappmpi.h"
#endif

#include "timer.h"


MainApp *global_app=NULL;
bool abort_once=false;




void closeFiles(int sig_num)
{
  if (!abort_once)
  {
    abort_once=true;
    printf("An abort evend was signaled......\n");
    global_app->abortEvent();
    printf("Trying to output the data......Ok\n");
  }
  else
    printf("Trying to output the data......Failed\n");
  printf("Closing Output file\ngoodby\n");
  HDF5OrthoWriter::getHDF5OrthoWriter().close();
 
}




void fpehandler(int sig_num)
{
        signal(SIGFPE, fpehandler);
        printf("SIGFPE %d: FLOATING POINT EXCEPTION OCCURRED\n",sig_num);
        printf("Segmentation Fault signal will be sent to cause\n");
        printf("production of a core file\n");
        closeFiles(0);
        kill(getpid(),SIGSEGV);
        abort();
}






int main(int argc,char **argv)
{
  

  //Build the message to print if bad arguments were passed
  char errorMessage[1000];
  
  //FlashCO2Brine::printFlashPropertiesTable(5,800,100,353,0.02);return 0;
  //FlashCO2Brine::printCO2RichMolarVolume(5,800,300,353);return 0;
  

  //enable arithmetic exception
  feenableexcept(FE_DIVBYZERO|FE_INVALID);
  signal(SIGFPE,fpehandler);
  signal(SIGKILL,closeFiles);
  signal(SIGABRT,closeFiles);
  signal(SIGTERM,closeFiles);

  sprintf(errorMessage,"Invalid arguments passed\n\
Usage: %s [Options] <CONFILE FILE NAME>\n\
Options:\n\
-g\tActivate unit tests mode\n",argv[0]);


  if (argc != 3 && argc!= 2)
  {
      printf("%s",errorMessage);
      return 0;
  }
  //read parameter -g and set the
  //variables bunit_test, configFile
  //depending on it
  std::string configFile;
  bool bunit_test=false;

  
  
  //check if -g option was passed
  if (argc == 3)
  {
    std::string option = argv[1];
    if (option == "-g")
      bunit_test=true;
    else
    {
      printf("%s",errorMessage);
      return 0;
    }
    configFile=argv[2];
  }
  else if (argc == 2)
    configFile=argv[1];

  
  

  // if the program arrived here, then the option requested
  // is to run the simulator with the configuration file
  // passed as argument. 
#ifndef USE_MPI
  Timer t1(PROCESSOR_TIME);
   t1.start();

  //Execute the MainApp class
  MainApp app;
  global_app = &app;
  app.execute(configFile,bunit_test);

  t1.end();
  std::cout << "Total Exec Time: " <<  t1.printTime() << std::endl;
  return 0;
  
#else

  
  MPI_Init(&argc,&argv);
  if (argc != 2)
   {
     printf("%s",errorMessage);
     return 0;
   }
  system("echo $HOSTNAME $SSH_HOST $SSH_CONNECTION");
  MainAppMPI app;
  app.execute(argv[1]);
  
  MPI_Finalize();
  return 0;

#endif


  

}









  /*
  //If first argument is the keyword "--ref" then
  //begin the process to coarse the mesh"
  if (argc == 7 && strcmp(argv[1],"--ref")==0)
  {
    unsigned nX,nY,nZ;
    char fileIn[1000],fileOut[1000];
    if (argc != 7)//check if the numbers of arguments are correct
    {
      printf("%s",errorMessage);
      return 0;
    }
    //Read the arguments
    if (sscanf(argv[2],"%d",&nX) != 1 ||
        sscanf(argv[3],"%d",&nY) != 1 ||
        sscanf(argv[4],"%d",&nZ) != 1 ||
        sscanf(argv[5],"%s",fileIn) != 1 ||
        sscanf(argv[6],"%s",fileOut) != 1)
    {
      printf("%s",errorMessage);
      return 0;
    }

    //Coarse the mesh
    MainApp app;
    app.coarseMesh(nX,nY,nZ,fileIn,fileOut);
    return 0;
  }
  */

---