magic_plot.cpp

Aller à la documentation de ce fichier.

//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//  MAGiC
//####//  Jean Christophe Cuilliere et Vincent FRANCOIS
//####//  Departement de Genie Mecanique - UQTR
//####//------------------------------------------------------------
//####//  MAGIC est un projet de recherche de l equipe ERICCA
//####//  du departement de genie mecanique de l Universite du Quebec a Trois Rivieres
//####//  http://www.uqtr.ca/ericca
//####//  http://www.uqtr.ca/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       magic_plot.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:57 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "magic_plot.h"
#include  <sys/types.h>
#include  <unistd.h>
 
 
 
MAGIC_PLOT::MAGIC_PLOT()
{
tenseur_couleur=VERT;
tenseur_titre="Titre";
tenseur_titrebarcouleur="Valeur";
tenseur_minbarcouleur=-1e301;
tenseur_maxbarcouleur=1e301;
tenseur_format1="%.4e";
tenseur_format2="%.4e";
tenseur_nom_fichier_eps="tenseur.eps";
 
histogramme_couleur=GRIS;
std::string histogramme_titre="Titre";
std::string histogramme_nom_fichier_eps="histogramme.eps";
std::string histogramme_axeX_nom="Nom axe x";
std::string histogramme_axeY_nom="Nom axe y";
}
 
MAGIC_PLOT::~MAGIC_PLOT()
{
}
 
void MAGIC_PLOT::plot_tenseur(int num,OT_TENSEUR tens1,OT_TENSEUR tens2,std::string output,bool aveccouleur)
{
plot_script_tenseur(num,tens1,tens2,output,aveccouleur);
char chaine[500];
sprintf(chaine,"gnuplot -p %s",output.c_str());
if (fork()==0) {system(chaine);exit(0);}
 
 
}
 
void MAGIC_PLOT::plot_script_tenseur(int num,OT_TENSEUR tens1,OT_TENSEUR tens2,std::string output,bool aveccouleur)
{
FILE *out=fopen(output.c_str(),"wt");
int n=tens1.Get_NbRows();
int m=tens1.Get_NbCols();
 
fprintf(out,"$data_moyenne << EOD\n");
for (int i=0;i<n;i++)
  {
  for (int j=0;j<m;j++)
    fprintf(out,"%e ",tens1(i,j).get_x());
  fprintf(out,"\n");
  }
fprintf(out,"EOD\n");
if (num==2)
{
fprintf(out,"$data_ecart_type << EOD\n");
for (int i=0;i<n;i++)
  {
  for (int j=0;j<m;j++)
    fprintf(out,"%e ",tens2(i,j).get_x());
  fprintf(out,"\n");
  }
fprintf(out,"EOD\n");
}
 
fprintf(out,"#set terminal postscript eps size 3.0,2.00 enhanced color\n");
fprintf(out,"#set output %s\n",tenseur_nom_fichier_eps.c_str());
 
fprintf(out,"set title \"%s\" font ',14'\n",tenseur_titre.c_str());
fprintf(out,"unset key\n");
fprintf(out,"set border linewidth 2\n");
if (!aveccouleur) fprintf(out,"set palette rgbformula 2,2,2\n");
else
  {
  if (tenseur_couleur==VERT) fprintf(out,"set palette rgbformula -3,2,-3\n");  
  if (tenseur_couleur==ROUGE) fprintf(out,"set palette rgbformula 2,-3,-3\n");  
  if (tenseur_couleur==BLEU) fprintf(out,"set palette rgbformula -3,-3,2\n");  
  if (tenseur_couleur==JAUNE) fprintf(out,"set palette rgbformula 2,2,-3\n");  
  if (tenseur_couleur==VIOLET) fprintf(out,"set palette rgbformula 2,-3,2\n");  
  if (tenseur_couleur==GRIS) fprintf(out,"set palette rgbformula -8,-8,-8\n");  
  
  }
 
 
fprintf(out,"unset xtics\n");
fprintf(out,"unset ytics\n");
if (aveccouleur) 
  {
  if ((tenseur_maxbarcouleur<1e300) && (tenseur_minbarcouleur>-1e300)) fprintf(out,"set cbrange [%lf:%lf]\n",tenseur_minbarcouleur,tenseur_maxbarcouleur);
  if ((tenseur_maxbarcouleur>1e300) && (tenseur_minbarcouleur<-1e300)) fprintf(out,"set cbrange []\n");
  if ((tenseur_maxbarcouleur<1e300) && (tenseur_minbarcouleur<-1e300)) fprintf(out,"set cbrange [:%lf]\n",tenseur_maxbarcouleur);
  if ((tenseur_maxbarcouleur>1e300) && (tenseur_minbarcouleur>-1e300)) fprintf(out,"set cbrange [%lf:]\n",tenseur_minbarcouleur);
  
  fprintf(out,"set cblabel \"%s\"\n",tenseur_titrebarcouleur.c_str());
  }
else fprintf(out,"unset colorbox\n");
  
fprintf(out,"set xrange []\n");
fprintf(out,"set yrange [] reverse\n");
 
for (int i=1;i<n;i++)
  fprintf(out,"set arrow from -0.5,%d.5 to %d.5,%d.5 nohead front lc rgb 'black' lw 1\n",i-1,m-1,i-1);
 
for (int j=1;j<m;j++)
  fprintf(out,"set arrow from %d.5,-0.5 to %d.5,%d.5 nohead front lc rgb 'black' lw 1\n",j-1,j-1,n-1);
 
fprintf(out,"plot '$data_moyenne' matrix with image,\\\n");
if (num==2)
    {
    fprintf(out,"'$data_moyenne' matrix using 1:2:(sprintf(\"%s\", $3)) with labels font ',14' offset 0,1 ,\\\n",tenseur_format1.c_str());
    fprintf(out,"'$data_ecart_type' matrix using 1:2:(sprintf(\"(%s)\", $3)) with labels font ',14' offset 0,-0.8\n",tenseur_format2.c_str());
    }
    else fprintf(out,"'$data_moyenne' matrix using 1:2:(sprintf(\"%s\", $3)) with labels font ',14'\n",tenseur_format1.c_str());
 
 
fclose(out);
}
 
void MAGIC_PLOT::plot_histogramme(std::string data_histogramme, std::string output, bool aveccouleur)
{
  plot_script_histogramme(data_histogramme,output,aveccouleur);
  char chaine[500];
  sprintf(chaine,"gnuplot -p %s",output.c_str());
  if (fork()==0) {system(chaine);exit(0);}
}
 
void MAGIC_PLOT::plot_script_histogramme(std::string data_histogramme, std::string output, bool aveccouleur)
{
  FILE *out=fopen(output.c_str(),"wt");
  fprintf(out,"#set terminal postscript eps size 3.0,2.00 enhanced color\n");
  fprintf(out,"#set output %s\n",histogramme_nom_fichier_eps.c_str());
  fprintf(out,"set title \"%s\" font ',14'\n",histogramme_titre.c_str());
  fprintf(out,"set xrange []\n");
  fprintf(out,"set yrange []\n");
  fprintf(out,"set xlabel '%s'\n",histogramme_axeX_nom.c_str());
  fprintf(out,"set ylabel '%s'\n",histogramme_axeY_nom.c_str());
  fprintf(out,"plot '%s' using 3:4 with impulses notitle lc rgb 'blue' lw 2",data_histogramme.c_str());
  fclose(out);
}