module_t3d.cpp

Go to the documentation of this file.

#include "module_t3d.h"
#include "problem.h"

#include "librw.h"
#include "gelib.h"

#ifdef __T3D_MODULE
#include "t3d.h"
#endif

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


namespace midaspace {

void mesh_generate_T3d (const Problem *pd, const char *solver, const char *options_string)
{
  char *auxch;
  
  //double sec = clock();
  //Pd->sodriver()->print_info_message (SODE_start_blue_arrow, "Model[%ld] - mesh generation by T3d", this->give_ID());
  
  
  fflush(stdout);
  //fpos_t pos;  fgetpos(stdout, &pos); // kdyz to oddebuguju, tak se zde promenna pos nezmeni, tim padem mi to dole ...
  int fd_out = dup(fileno(stdout));
  freopen(auxch = strdupl(pd->give_OUT_moFILE()->give_name(), ".T3d.log.stdout"), "w", stdout);
  delete [] auxch;
  
  fflush(stderr);
  //fpos_t pos;  fgetpos(stderr, &pos); // kdyz to oddebuguju, tak se zde promenna pos nezmeni, tim padem mi to dole ...
  int fd_err = dup(fileno(stderr));
  freopen(auxch = strdupl(pd->give_OUT_moFILE()->give_name(), ".T3d.log.stderr"), "w", stderr);
  delete [] auxch;
  
  
  int status;
  if (solver) {
    char cmd[1000];
    sprintf(cmd, "%s %s", solver, options_string);
    //fprintf (stdout,"%s\n", cmd);
    
    fprintf (stdout,"\n");
    fflush(stdout);  fflush(stderr);
    status = system(cmd);
    fprintf(stdout,"\n");
  }
  else {
#ifdef __T3D_MODULE
    try {
      t3d_main (NULL, (char *)options_string);
      status = 0;
    }
    catch (...){
      status = 1;
    }
#else
    _errorr("T3d binary is given as \"-\", however the module __T3D_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);  // ... rve, ze pos neni inicializovana, proto tento radek zatim zakomentovavam
  
  fflush(stderr);
  dup2(fd_err, fileno(stderr));
  close(fd_err);
  clearerr(stderr);
  //fsetpos(stderr, &pos);  // ... rve, ze pos neni inicializovana, proto tento radek zatim zakomentovavam
  
  if (status) _errorr("T3d ma problemy Houstone");
  //print_info_message_timeok ((clock() - sec) / (double)CLOCKS_PER_SEC);
}

} // namespace midaspace

// these functions are have to be outside of midspace
// nejaky problemy s t3d funkcema, casem to oprav aby to tady nebylo
void print_message(FILE *stream, char *buffer)
{
  fprintf (stream, "%s", buffer);
}
void exit_termination(int exit_code)
{
  fflush (stdout);
  CHANGE_CONSOLE_COLOUR(stderr, TC_D_RED);     fprintf (stderr,"\n\n T3d  ERROR:  exit code %d", exit_code);
  CHANGE_CONSOLE_COLOUR(stderr, TC_DEFAULT);   fprintf (stderr," !!!\n\n");
}
void process_exit(int exit_code)
{
  throw 0;
}