module_oofem.cpp

Go to the documentation of this file.

#include "module_oofem.h"
#include "problem.h"

#include "librw.h"

#ifdef __OOFEM_MODULE
#include "mainlib.h"
#endif

#include <stdio.h>
#include <stdlib.h>
#include <unistd.h>
#include <time.h>


namespace midaspace {

void FEM_solve_OOFEM (Problem *pd)
{
  char *auxch;
  
  double sec = clock();
  pd->sodriver()->print_info_message (SODE_start_blue_arrow, "OOFEM analyse");
  
  
  fflush(stdout);
  //fpos_t pos;  fgetpos(stdout, &pos);
  int fd_out = dup(fileno(stdout));
  freopen(auxch = strdupl(pd->give_OUT_moFILE()->give_name(), ".oofem.log.stdout"), "w", stdout);
  delete [] auxch;
  
  fflush(stderr);
  //fpos_t pos;  fgetpos(stderr, &pos);
  int fd_err = dup(fileno(stderr));
  freopen(auxch = strdupl(pd->give_OUT_moFILE()->give_name(), ".oofem.log.stderr"), "w", stderr);
  delete [] auxch;
  
  const char *solver = pd->give_P_solverbinary()->give_name_or_null();
  
  int status;
  if (solver) {
    char cmd[1000];
    sprintf(cmd, "%s", solver);
    
    //argc = 1;
    //argv[argc++] = strdupl ("-rn");
    //argv[argc++] = strdupl ("-f");
    //argv[argc++] = strdupl (this->OUT_moFILE->give_name(), ".oofem.in");
    //argv[argc++] = strdupl ("-qo");
    //argv[argc++] = strdupl (this->OUT_moFILE->give_name(), ".oofem.log");
    //for (int i=1; i<argc; i++)  { strcat(cmd, " "); strcat(cmd, args[i]); }
    
    strcat(cmd, " -rn");
    strcat(cmd, " -f ");   strcat(cmd, pd->give_OUT_moFILE()->give_name());  strcat(cmd, ".oofem.in");
    strcat(cmd, " -qo ");  strcat(cmd, pd->give_OUT_moFILE()->give_name());  strcat(cmd, ".oofem.log");
    
    fflush(stdout);  fflush(stderr);
    status = system(cmd);
    fprintf(stdout,"\n");
  }
  else {
#ifdef __OOFEM_MODULE
    char fname[256];
    sprintf(fname, "%s%s", pd->give_OUT_moFILE()->give_name(), ".oofem.in");
    
    char oname[256];
    sprintf(oname, "%s%s", pd->give_OUT_moFILE()->give_name(), ".oofem.log");
    
    try {
      //oofem_main (, (char **)args, &pd->solver_file_ptr, &pd->solver_file_size);
      oofem_main (true, fname, oname);
      status = 0;
    }
    catch (...){
      status = 1;
    }
    
    // int status = oofem_main (argc, (char **)args, &pd->solver_file_ptr, &pd->solver_file_size);
    //if (status)  _errorr("OOFEM ma problemy Houstone");
#else
    _errorr("OOFEM binary is given as \"-\", however the module __OOFEM_MODULE is not linked. Activate the module in Makefile please.");
#endif
  }
  
  fflush(stdout);
  dup2(fd_out, fileno(stdout));
  close(fd_out);
  clearerr(stdout);
  //fsetpos(stdout, &pos);
  
  fflush(stderr);
  dup2(fd_err, fileno(stderr));
  close(fd_err);
  clearerr(stderr);
  //fsetpos(stderr, &pos);
  
  if (status) _errorr("OOFEM ma problemy Houstone");
  
  pd->sodriver()->print_info_time_green_ok (SODE_end_green_ok, (clock() - sec) / (double)CLOCKS_PER_SEC);
}

} // namespace midaspace