mstruct_analyse_polycristaux.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       mstruct_analyse_polycristaux.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:57 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "gestionversion.h"
#include "mstruct_analyse_polycristaux.h"
#include <math.h>
#include <mutex>
 
 
MSTRUCT_POLY_ANALYSE_RESULTAT::MSTRUCT_POLY_ANALYSE_RESULTAT(char *nomparam)
{
param.lire(nomparam);
}
 
 
MSTRUCT_POLY_ANALYSE_RESULTAT::~MSTRUCT_POLY_ANALYSE_RESULTAT()
{
}
 
 
 
MSTRUCT_POLY_RESULTAT::~MSTRUCT_POLY_RESULTAT()
{
}
 
 
int MSTRUCT_POLY_ANALYSE_RESULTAT::ident(bool verif)
{
int nbessai=0;
for (int i=0;i<tab.size();i++)
    {
    bool essaivalide=false;
    if (verif==false) essaivalide=true;
    if (((tab[i].ECH-tab[i].Evoigt)*(tab[i].ECH-tab[i].Ereuss) < 0) && ((tab[i].EDH-tab[i].Evoigt)*(tab[i].EDH-tab[i].Ereuss) < 0)) essaivalide=true;
    if ((essaivalide==true) && (nbessai==0))
        {
        tab[0].essai=0;
        nbessai++;
        std::vector<int> lst;
        listeessai.push_back(lst);
        listeessai[nbessai-1].push_back(i);
        }
    else
        if (essaivalide==true)
        {
        bool trouve=false;
        for (int j=0;j<listeessai.size();j++)
            if (tab[i].nbcristaux==tab[listeessai[j][0]].nbcristaux)
            if (tab[i].degre==tab[listeessai[j][0]].degre)
            if (tab[i].nbphase==tab[listeessai[j][0]].nbphase)
                {
                tab[i].essai= tab[listeessai[j][0]].essai;
                listeessai[tab[i].essai].push_back(i);
                trouve=true;
                }
        if (trouve==false)
                {
                tab[i].essai=nbessai;
                nbessai++;
                std::vector<int> lst;
                listeessai.push_back(lst);
                listeessai[nbessai-1].push_back(i);
                }
        }
        else listeessairatee.push_back(i);
    }
return listeessai.size();
    
}
 
 
void MSTRUCT_POLY_ANALYSE_RESULTAT::execute(int nbessai)
{
OT_VECTEUR_3D zero(0.,0.,0.);
cdmm=zero;cdme=zero;
KDH.zero();GDH.zero();EDH.zero();
KCH.zero();GCH.zero();ECH.zero();
masse.zero();volume.zero();
cumcdm.resize(get_nbechantillon(nbessai));
int nbph=tab[listeessai[nbessai][0]].nbphase;
cdmphase.resize(nbph);
cdmephase.resize(nbph);
cumcdmphase.resize(nbph);
masseph.resize(nbph);
volumeph.resize(nbph);
tens1m.reinit(3,3);tens2m.reinit(3,3);tens3m.reinit(3,3);
tens1e.reinit(3,3);tens2e.reinit(3,3);tens3e.reinit(3,3);
for (int i=0;i<3;i++)
    for (int j=0;j<3;j++)
        {
        tens1m(i,j)=0.;tens2m(i,j)=0.;tens3m(i,j)=0.;    
        tens1e(i,j)=0.;tens2e(i,j)=0.;tens3e(i,j)=0.;    
        }
for (int i=0;i<nbph;i++)
    {
    cdmphase[i]=zero;
    cdmephase[i]=zero;
    masseph[i].zero();
    volumeph[i].zero();
    cumcdmphase[i].resize(get_nbechantillon(nbessai));
    cumEDH.resize(get_nbechantillon(nbessai));
    cumECH.resize(get_nbechantillon(nbessai));
    cumiso.resize(get_nbechantillon(nbessai));
    }
for (int i=0;i<get_nbechantillon(nbessai);i++)
    {
    KDH.moyenne=KDH.moyenne+tab[listeessai[nbessai][i]].KDH;
    GDH.moyenne=GDH.moyenne+tab[listeessai[nbessai][i]].GDH;
    EDH.moyenne=EDH.moyenne+tab[listeessai[nbessai][i]].EDH;
    KCH.moyenne=KCH.moyenne+tab[listeessai[nbessai][i]].KCH;
    GCH.moyenne=GCH.moyenne+tab[listeessai[nbessai][i]].GCH;
    ECH.moyenne=ECH.moyenne+tab[listeessai[nbessai][i]].ECH;
    masse.moyenne=masse.moyenne+tab[listeessai[nbessai][i]].masse;
    volume.moyenne=volume.moyenne+tab[listeessai[nbessai][i]].volume;
    cdmm=cdmm+tab[listeessai[nbessai][i]].cdm;
    cumcdm[i]=cdmm/((double)i+1);
    cumEDH[i]=EDH.moyenne/((double)i+1);
    cumECH[i]=ECH.moyenne/((double)i+1);
    for (int j=0;j<nbph;j++)
        {
         volumeph[j].moyenne=volumeph[j].moyenne+tab[listeessai[nbessai][i]].volumephase[j];
         masseph[j].moyenne=masseph[j].moyenne+tab[listeessai[nbessai][i]].massephase[j];
         cdmphase[j]=cdmphase[j]+tab[listeessai[nbessai][i]].cdmphase[j];
         cumcdmphase[j][i]=cdmphase[j]/((double)i+1);  
        }
    for (int ii=0;ii<3;ii++)
        for (int jj=0;jj<3;jj++)
            {
            tens1m(ii,jj)=tens1m(ii,jj)+tab[listeessai[nbessai][i]].tens1(ii,jj);    
            tens2m(ii,jj)=tens2m(ii,jj)+tab[listeessai[nbessai][i]].tens1(ii,jj);    
            tens3m(ii,jj)=tens3m(ii,jj)+tab[listeessai[nbessai][i]].tens1(ii,jj);    
            }
    }
KDH.moyenne=KDH.moyenne/(double)get_nbechantillon(nbessai);
GDH.moyenne=GDH.moyenne/(double)get_nbechantillon(nbessai);
EDH.moyenne=EDH.moyenne/(double)get_nbechantillon(nbessai);
KCH.moyenne=KCH.moyenne/(double)get_nbechantillon(nbessai);
GCH.moyenne=GCH.moyenne/(double)get_nbechantillon(nbessai);
ECH.moyenne=ECH.moyenne/(double)get_nbechantillon(nbessai);
masse.moyenne=masse.moyenne/(double)get_nbechantillon(nbessai);
volume.moyenne=volume.moyenne/(double)get_nbechantillon(nbessai);
cdmm=cdmm/(double)get_nbechantillon(nbessai);
for (int j=0;j<nbph;j++)
    {
    masseph[j].moyenne=masseph[j].moyenne/(double)get_nbechantillon(nbessai); 
    volumeph[j].moyenne=volumeph[j].moyenne/(double)get_nbechantillon(nbessai); 
    cdmphase[j]=cdmphase[j]/(double)get_nbechantillon(nbessai); 
    }
for (int ii=0;ii<3;ii++)
    for (int jj=0;jj<3;jj++)
            {
            tens1m(ii,jj)=tens1m(ii,jj)/(double)get_nbechantillon(nbessai);    
            tens2m(ii,jj)=tens2m(ii,jj)/(double)get_nbechantillon(nbessai);    
            tens3m(ii,jj)=tens3m(ii,jj)/(double)get_nbechantillon(nbessai);    
            }
for (int i=0;i<get_nbechantillon(nbessai);i++)
    {
    KDH.ecart=KDH.ecart+(tab[listeessai[nbessai][i]].KDH-KDH.moyenne)*(tab[listeessai[nbessai][i]].KDH-KDH.moyenne);
    GDH.ecart=GDH.ecart+(tab[listeessai[nbessai][i]].GDH-GDH.moyenne)*(tab[listeessai[nbessai][i]].GDH-GDH.moyenne);
    EDH.ecart=EDH.ecart+(tab[listeessai[nbessai][i]].EDH-EDH.moyenne)*(tab[listeessai[nbessai][i]].EDH-EDH.moyenne);
    KCH.ecart=KCH.ecart+(tab[listeessai[nbessai][i]].KDH-KCH.moyenne)*(tab[listeessai[nbessai][i]].KCH-KCH.moyenne);
    GCH.ecart=GCH.ecart+(tab[listeessai[nbessai][i]].GDH-GCH.moyenne)*(tab[listeessai[nbessai][i]].GCH-GCH.moyenne);
    ECH.ecart=ECH.ecart+(tab[listeessai[nbessai][i]].EDH-ECH.moyenne)*(tab[listeessai[nbessai][i]].ECH-ECH.moyenne);
    masse.ecart=masse.ecart+(tab[listeessai[nbessai][i]].masse-masse.moyenne)*(tab[listeessai[nbessai][i]].masse-masse.moyenne);
    volume.ecart=volume.ecart+(tab[listeessai[nbessai][i]].volume-volume.moyenne)*(tab[listeessai[nbessai][i]].volume-volume.moyenne);
    cdme.change_x(cdme.get_x()+(tab[listeessai[nbessai][i]].cdm.get_x()-cdmm.get_x())*(tab[listeessai[nbessai][i]].cdm.get_x()-cdmm.get_x()));
    cdme.change_y(cdme.get_y()+(tab[listeessai[nbessai][i]].cdm.get_y()-cdmm.get_y())*(tab[listeessai[nbessai][i]].cdm.get_y()-cdmm.get_y()));
    cdme.change_z(cdme.get_z()+(tab[listeessai[nbessai][i]].cdm.get_z()-cdmm.get_z())*(tab[listeessai[nbessai][i]].cdm.get_z()-cdmm.get_z()));
    for (int j=0;j<nbph;j++)
        {
        masseph[j].ecart=masseph[j].ecart+(tab[listeessai[nbessai][i]].massephase[j]-masseph[j].moyenne)*(tab[listeessai[nbessai][i]].massephase[j]-masseph[j].moyenne);
        volumeph[j].ecart=volumeph[j].ecart+(tab[listeessai[nbessai][i]].volumephase[j]-volumeph[j].moyenne)*(tab[listeessai[nbessai][i]].volumephase[j]-volumeph[j].moyenne);
        cdmephase[j].change_x(cdmephase[j].get_x()+(tab[listeessai[nbessai][i]].cdmphase[j].get_x()-cdmphase[j].get_x())*(tab[listeessai[nbessai][i]].cdmphase[j].get_x()-cdmphase[j].get_x()));
        cdmephase[j].change_y(cdmephase[j].get_y()+(tab[listeessai[nbessai][i]].cdmphase[j].get_y()-cdmphase[j].get_x())*(tab[listeessai[nbessai][i]].cdmphase[j].get_y()-cdmphase[j].get_y()));
        cdmephase[j].change_z(cdmephase[j].get_z()+(tab[listeessai[nbessai][i]].cdmphase[j].get_z()-cdmphase[j].get_x())*(tab[listeessai[nbessai][i]].cdmphase[j].get_z()-cdmphase[j].get_z()));
        }
    for (int ii=0;ii<3;ii++)
        for (int jj=0;jj<3;jj++)
            {
            tens1e(ii,jj)=tens1e(ii,jj)+(tab[listeessai[nbessai][i]].tens1(ii,jj)-tens1m(ii,jj))*(tab[listeessai[nbessai][i]].tens1(ii,jj)-tens1m(ii,jj));    
            tens2e(ii,jj)=tens2e(ii,jj)+(tab[listeessai[nbessai][i]].tens2(ii,jj)-tens2m(ii,jj))*(tab[listeessai[nbessai][i]].tens2(ii,jj)-tens2m(ii,jj));    
            tens3e(ii,jj)=tens3e(ii,jj)+(tab[listeessai[nbessai][i]].tens3(ii,jj)-tens3m(ii,jj))*(tab[listeessai[nbessai][i]].tens3(ii,jj)-tens3m(ii,jj));    
            }
 
    }
KDH.ecart=sqrt(KDH.ecart/(double)get_nbechantillon(nbessai));   
GDH.ecart=sqrt(GDH.ecart/(double)get_nbechantillon(nbessai));   
EDH.ecart=sqrt(EDH.ecart/(double)get_nbechantillon(nbessai));   
KCH.ecart=sqrt(KCH.ecart/(double)get_nbechantillon(nbessai));   
GCH.ecart=sqrt(GCH.ecart/(double)get_nbechantillon(nbessai));   
ECH.ecart=sqrt(ECH.ecart/(double)get_nbechantillon(nbessai));   
 
masse.ecart=sqrt(masse.ecart/(double)get_nbechantillon(nbessai));    
volume.ecart=sqrt(volume.ecart/(double)get_nbechantillon(nbessai));
cdme.change_x(sqrt(cdme.get_x()/(double)get_nbechantillon(nbessai)));
cdme.change_y(sqrt(cdme.get_y()/(double)get_nbechantillon(nbessai)));
cdme.change_z(sqrt(cdme.get_z()/(double)get_nbechantillon(nbessai)));
for (int j=0;j<nbph;j++)
{
    volumeph[j].ecart=sqrt(volumeph[j].ecart/(double)get_nbechantillon(nbessai));
    masseph[j].ecart=sqrt(masseph[j].ecart/(double)get_nbechantillon(nbessai));
    cdmephase[j].change_x(sqrt(cdmephase[j].get_x()/(double)get_nbechantillon(nbessai)));
    cdmephase[j].change_y(sqrt(cdmephase[j].get_y()/(double)get_nbechantillon(nbessai)));
    cdmephase[j].change_z(sqrt(cdmephase[j].get_z()/(double)get_nbechantillon(nbessai)));
}
for (int ii=0;ii<3;ii++)
    for (int jj=0;jj<3;jj++)
        {
        tens1e(ii,jj)=sqrt(tens1e(ii,jj)/(double)get_nbechantillon(nbessai));    
        tens2e(ii,jj)=sqrt(tens2e(ii,jj)/(double)get_nbechantillon(nbessai));    
        tens3e(ii,jj)=sqrt(tens3e(ii,jj)/(double)get_nbechantillon(nbessai));    
        }
 
            
char nom[5000];
std::string nomcumul=param.get_nom("fichiercumul");
strcpy(nom,nomcumul.c_str());
enregistrer(nbessai,nom);
}
 
 
 
int MSTRUCT_POLY_ANALYSE_RESULTAT::lire(void)
{
std::mutex mtx;
char nom[5000];
std::string nomcumul=param.get_nom("fichiercumul");
strcpy(nom,nomcumul.c_str());
strcat(nom,".txt");
mtx.lock();
FILE* in=fopen(nom,"rt");
    char mess[5000];
    fgets(mess,5000,in);
    while (feof(in)==false)
    {
    fgets(mess,5000,in);
    if (feof(in)==true) continue;
    char *p=mess;
    int ok=0;
    std::vector<double> valeurs;
    while (ok==0)
        {
        double val;
        sscanf(p,"%lf;" ,&val); 
        valeurs.push_back(val);
        p=strchr(p,';');
        if (p==NULL) ok=1;
        else p++;
        }
    
    MSTRUCT_POLY_RESULTAT res;
    res.nbcristaux=(int)(valeurs[0]);
    res.degre=(int)(valeurs[1]);
    res.dg=valeurs[2];
    res.KDH=valeurs[3];
    res.GDH=valeurs[4];
    res.EDH=valeurs[5];
    res.KCH=valeurs[6];
    res.GCH=valeurs[7];
    res.ECH=valeurs[8];
    res.tens1(0,0)=valeurs[9];
    res.tens1(1,0)=valeurs[10];
    res.tens1(2,0)=valeurs[11];
    res.tens1(0,1)=valeurs[12];
    res.tens1(1,1)=valeurs[13];
    res.tens1(2,1)=valeurs[14];
    res.tens1(0,2)=valeurs[15];
    res.tens1(1,2)=valeurs[16];
    res.tens1(2,2)=valeurs[17];
    res.tens2(0,0)=valeurs[18];
    res.tens2(1,0)=valeurs[19];
    res.tens2(2,0)=valeurs[20];
    res.tens2(0,1)=valeurs[21];
    res.tens2(1,1)=valeurs[22];
    res.tens2(2,1)=valeurs[23];
    res.tens2(0,2)=valeurs[24];
    res.tens2(1,2)=valeurs[25];
    res.tens2(2,2)=valeurs[26];
    res.tens3(0,0)=valeurs[27];
    res.tens3(1,0)=valeurs[28];
    res.tens3(2,0)=valeurs[29];
    res.tens3(0,1)=valeurs[30];
    res.tens3(1,1)=valeurs[31];
    res.tens3(2,1)=valeurs[32];
    res.tens3(0,2)=valeurs[33];
    res.tens3(1,2)=valeurs[34];
    res.tens3(2,2)=valeurs[35];
    res.nbphase=(int)(valeurs[36]);
    for (int i=0;i<res.nbphase;i++)
            {
            res.cristauxphase.push_back((int)valeurs[37+i*6]);
            res.volumephase.push_back(valeurs[38+i*6]);
            res.massephase.push_back(valeurs[39+i*6]);
            OT_VECTEUR_3D cdm(valeurs[40+i*6],valeurs[41+i*6],valeurs[42+i*6]);
            res.cdmphase.push_back(cdm);
            }
    res.volume=valeurs[42+(res.nbphase-1)*6+1];
    res.masse=valeurs[42+(res.nbphase-1)*6+2];
    OT_VECTEUR_3D cdm(valeurs[42+(res.nbphase-1)*6+3],valeurs[42+(res.nbphase-1)*6+4],valeurs[42+(res.nbphase-1)*6+5]);
    res.cdm=cdm;
    res.Kvoigt=valeurs[42+(res.nbphase-1)*6+6];
    res.Gvoigt=valeurs[42+(res.nbphase-1)*6+7];
    res.Evoigt=valeurs[42+(res.nbphase-1)*6+8];
    res.nuvoigt=valeurs[42+(res.nbphase-1)*6+9];
    res.Kreuss=valeurs[42+(res.nbphase-1)*6+10];
    res.Greuss=valeurs[42+(res.nbphase-1)*6+11];
    res.Ereuss=valeurs[42+(res.nbphase-1)*6+12];
    res.nureuss=valeurs[42+(res.nbphase-1)*6+13];
    tab.push_back(res);
    }
fclose(in);
mtx.unlock();
return tab.size();
}
 
void MSTRUCT_POLY_ANALYSE_RESULTAT::enregistrer(int nbessai,char* nom)
{
char nomout[5000];
sprintf(nomout,"%s%d.txt",nom,nbessai+1);
FILE* out=fopen(nomout,"wt");
fprintf(out,"----------------------------------------------------------\n");
fprintf(out,"nb cristaux = %d\n",tab[listeessai[nbessai][0]].nbcristaux);
fprintf(out,"nb phases = %d\n",tab[listeessai[nbessai][0]].nbphase);
fprintf(out,"degre = %d\n",tab[listeessai[nbessai][0]].degre);
fprintf(out,"nb echantillons = %d\n",get_nbechantillon(nbessai));
fprintf(out,"masse = %lf +- %lf \n",masse.moyenne,masse.ecart);
fprintf(out,"volume = %lf +- %lf \n",volume.moyenne,volume.ecart);
fprintf(out,"CDM : x=%lf +- %lf ,y=%lf +- %lf ,z=%lf +- %lf \n",cdmm.get_x(),cdme.get_x(),cdmm.get_y(),cdme.get_y(),cdmm.get_z(),cdme.get_z());
for (int i=0;i<tab[listeessai[nbessai][0]].nbphase;i++)
    {
    fprintf(out,"masse phase %d = %lf +- %lf \n",i+1,masseph[i].moyenne,masseph[i].ecart);
    fprintf(out,"volume phase %d  = %lf +- %lf \n",i+1,volumeph[i].moyenne,volumeph[i].ecart);
    fprintf(out,"CDM phase %d : x=%lf +- %lf ,y=%lf +- %lf ,z=%lf +- %lf \n",i+1,cdmphase[i].get_x(),cdmephase[i].get_x(),cdmphase[i].get_y(),cdmephase[i].get_y(),cdmphase[i].get_z(),cdmephase[i].get_z());
    }
fprintf(out,"----------------------------------------------------------\n");
fprintf(out,"KDH = %lf + %lf \n",KDH.moyenne,KDH.ecart);
fprintf(out,"GDH = %lf + %lf \n",GDH.moyenne,GDH.ecart);
fprintf(out,"EDH = %lf + %lf \n",EDH.moyenne,EDH.ecart);
fprintf(out,"KCH = %lf + %lf \n",KCH.moyenne,KCH.ecart);
fprintf(out,"GCH = %lf + %lf \n",GCH.moyenne,GCH.ecart);
fprintf(out,"ECH = %lf + %lf \n",ECH.moyenne,ECH.ecart);
fprintf(out,"----------------------------------------------------------\n");
fprintf(out,"tenseur orientation axe 1\n");
for (int ii=0;ii<3;ii++)
    {
    for (int jj=0;jj<3;jj++)
        fprintf(out,"%lf +- %lf  ",tens1m(ii,jj).get_x(),tens1e(ii,jj).get_x());
    fprintf(out,"\n");        
    }
fprintf(out,"----------------------------------------------------------\n");
fprintf(out,"tenseur orientation axe 2\n");
for (int ii=0;ii<3;ii++)
    {
    for (int jj=0;jj<3;jj++)
        fprintf(out,"%lf +- %lf  ",tens2m(ii,jj).get_x(),tens2e(ii,jj).get_x());
    fprintf(out,"\n");        
    }
fprintf(out,"----------------------------------------------------------\n");
fprintf(out,"tenseur orientation axe 3\n");
for (int ii=0;ii<3;ii++)
    {
    for (int jj=0;jj<3;jj++)
        fprintf(out,"%lf +- %lf  ",tens3m(ii,jj).get_x(),tens3e(ii,jj).get_x());
    fprintf(out,"\n");        
    }
fprintf(out,"----------------------------------------------------------\n");
 
 
        
 
 
 
 
for (int i=0;i<get_nbechantillon(nbessai);i++)
    {
    int n=listeessai[nbessai][i];    
    fprintf(out,";;;;;;;;;%d;%lf;%lf;%lf;%lf;%lf;%lf;;%lf;%lf;;",i+1,tab[n].KDH,tab[n].GDH,tab[n].EDH,tab[n].KCH,tab[n].GCH,tab[n].ECH,cumEDH[i],cumECH[i]);
    OT_TENSEUR t(3,3);
    t.identite();
    double2 cst(1./3.);
    t=cst*t;
    double ecartiso1=0.;
    double ecartiso2=0.;
    double ecartiso3=0.;
    for (int ii=0;ii<3;ii++)
        for (int jj=0;jj<3;jj++)
            {
            ecartiso1=ecartiso1+(tab[n].tens1(ii,jj).get_x()-t(ii,jj).get_x())*(tab[n].tens1(ii,jj).get_x()-t(ii,jj).get_x());
            ecartiso2=ecartiso2+(tab[n].tens2(ii,jj).get_x()-t(ii,jj).get_x())*(tab[n].tens2(ii,jj).get_x()-t(ii,jj).get_x());
            ecartiso3=ecartiso3+(tab[n].tens3(ii,jj).get_x()-t(ii,jj).get_x())*(tab[n].tens3(ii,jj).get_x()-t(ii,jj).get_x());
            }
    double ecartiso=1.-(sqrt(ecartiso1)+sqrt(ecartiso2)+sqrt(ecartiso3))/27.;
    if (i==0) cumiso[i]=ecartiso; else cumiso[i]=cumiso[i-1]+ecartiso;
    fprintf(out,"%lf;%lf;;",ecartiso,cumiso[i]/((double)i+1));
    fprintf(out,"%f;%f;%lf;%lf;%lf;%lf;%lf;%lf;;",tab[n].volume,tab[n].masse,tab[n].cdm.get_x(),tab[n].cdm.get_y(),tab[n].cdm.get_z(),cumcdm[i].get_x(),cumcdm[i].get_y(),cumcdm[i].get_z());
    for (int j=0;j<tab[listeessai[nbessai][0]].nbphase;j++)
        fprintf(out,"%lf;%lf;%lf;%lf;%lf;%lf;%lf;%lf;;",tab[n].volumephase[j],tab[n].massephase[j],tab[n].cdmphase[j].get_x(),tab[n].cdmphase[j].get_y(),tab[n].cdmphase[j].get_z(),cumcdmphase[j][i].get_x(),cumcdmphase[j][i].get_y(),cumcdmphase[j][i].get_z());
    fprintf(out,";;%lf;%lf;%lf;%lf;%lf;%lf",tab[n].Evoigt,tab[n].Ereuss,tab[n].Kvoigt,tab[n].Kreuss,tab[n].Gvoigt,tab[n].Greuss);
    fprintf(out,"\n");
        
    }
fclose(out);
enregistrer_gnu(nbessai,nom);
    
}
 
 
void MSTRUCT_POLY_ANALYSE_RESULTAT::enregistrer_gnu(int nbessai,char* nom)
    {
    double xmin=param.get_valeur("boitexmin");
        double ymin=param.get_valeur("boiteymin");
        double zmin=param.get_valeur("boitezmin");
        double xmax=param.get_valeur("boitexmax");
        double ymax=param.get_valeur("boiteymax");
        double zmax=param.get_valeur("boitezmax");
    int degre=tab[listeessai[nbessai][0]].degre;
    char typecalcul[500];
    if (degre==1) strcpy(typecalcul,"Linéaire");
    if (degre==2) strcpy(typecalcul,"Quadratique");
    double xcible=0.5*(xmin+xmax);
    double ycible=0.5*(ymin+ymax);
    double zcible=0.5*(zmin+zmax);
    char execute[5000];
    char nomdata[5000];
    char nomout[1000];
    char nomoutsvg[5000];
    //----------------------------------------K
    sprintf(nomdata,"%s%d.txt",nom,nbessai+1);
    sprintf(nomout,"%s%dK.gnu",nom,nbessai+1);
    sprintf(nomoutsvg,"%s%dK.svg",nom,nbessai+1);
    FILE *out=fopen((char*)nomout,"wt");  
        fprintf(out, "set terminal svg\n");
        fprintf(out, "set output \"%s\"\n", nomoutsvg);
        fprintf(out, "set datafile separator \";\"\n");
        fprintf(out, "set key on  outside bottom center horizontal\n");
        fprintf(out, "set xlabel 'Numéro d''essai'\n");
        fprintf(out, "set ylabel 'Module de compressibilité K (GPa)'\n");
        fprintf(out, "set y2label 'Écart à l''isotropie'\n");
        fprintf(out, "set pointsize 1\n");
        fprintf(out, "set autoscale x\n");
        //fprintf(out, "set autoscale y\n");
        fprintf(out,"stats '%s' using ($%d/1e9) name 'YX'\n",nomdata,37+9*tab[listeessai[nbessai][0]].nbphase);
    fprintf(out,"stats '%s' using ($%d/1e9) name 'YN'\n",nomdata,38+9*tab[listeessai[nbessai][0]].nbphase);
    fprintf(out,"delta = (YX_max - YN_min) * 0.05\n");
    fprintf(out,"ymin_zoom = YN_min - delta\n");
    fprintf(out,"ymax_zoom = YX_max + delta\n");
    fprintf(out,"set yrange [ymin_zoom:ymax_zoom]\n");
    fprintf(out, "set y2tics 0.5,0.1,1\n");
        fprintf(out, "set y2r [0.8:1]\n");
        fprintf(out, "set ytics nomirror\n");
        fprintf(out, "set xtics nomirror\n");
        fprintf(out, "#\n");
        fprintf(out, "d(x)=dd\n");
        fprintf(out, "fit d(x) '%s' u :($11/1e9) via dd\n",nomdata);
        fprintf(out, "c(x)=cc\n");
    fprintf(out, "fit c(x) '%s' u :($14/1e9) via cc\n",nomdata);
    fprintf(out, "is(x)=is\n");
    fprintf(out, "fit is(x) '%s' u :($21) via is\n",nomdata);
    fprintf(out, "#\n");
    fprintf(out, "set title \"Résultat Homogéinisation \\n %s\\n %d cristaux - Calcul %s\"\n",param.get_nom("nommateriau").c_str(), tab[listeessai[nbessai][0]].nbcristaux,typecalcul);
    fprintf(out, "plot '%s' u 10 :($11/1e9) title 'KDH' pointtype 7 lc rgb \"blue\", '%s' u 10 :($14/1e9) title 'KCH' pointtype 7 lc rgb \"red\",'%s' u 10 : ($%d/1e9) title 'Borne Voigt' w l lc rgb \"brown\" lw 2,'%s' u 10 : ($%d/1e9) title 'Borne Reuss' w l lc rgb \"orange\" lw 2,  '%s'  u 10:21 axis x1y2 title 'Isotropie' lc rgb \"green\" pointtype 3,c(x) title sprintf(\"moyenne %%.1f\",cc) lc rgb \"red\" dt 2 ,d(x) title sprintf(\"moyenne %%.1f\",dd) lc rgb \"blue\" dt 2, is(x) title sprintf(\"moyenne %%.2f\",is) axis x1y2 lc rgb \"green\" dt 2\n", nomdata,nomdata, nomdata,37+9*tab[listeessai[nbessai][0]].nbphase, nomdata,38+9*tab[listeessai[nbessai][0]].nbphase,nomdata);
        fclose(out);
    sprintf(execute,"gnuplot %s  > %s.log 2>&1",nomout,nomout);
    system(execute);
        //----------------------------------------G
    sprintf(nomdata,"%s%d.txt",nom,nbessai+1);
    sprintf(nomout,"%s%dG.gnu",nom,nbessai+1);
    sprintf(nomoutsvg,"%s%dG.svg",nom,nbessai+1);
    out=fopen((char*)nomout,"wt");  
        fprintf(out, "set terminal svg\n");
        fprintf(out, "set output \"%s\"\n", nomoutsvg);
        fprintf(out, "set datafile separator \";\"\n");
        fprintf(out, "set key on  outside bottom center horizontal\n");
        fprintf(out, "set xlabel 'Numéro d''essai'\n");
        fprintf(out, "set ylabel 'Module de Coulomb G (GPa)'\n");
        fprintf(out, "set y2label 'Écart à l''isotropie'\n");
        fprintf(out, "set pointsize 1\n");
        fprintf(out, "set autoscale x\n");
        //fprintf(out, "set autoscale y\n");
        fprintf(out,"stats '%s' using ($%d/1e9) name 'YX'\n",nomdata,39+9*tab[listeessai[nbessai][0]].nbphase);
    fprintf(out,"stats '%s' using ($%d/1e9) name 'YN'\n",nomdata,40+9*tab[listeessai[nbessai][0]].nbphase);
    fprintf(out,"delta = (YX_max - YN_min) * 0.05\n");
    fprintf(out,"ymin_zoom = YN_min - delta\n");
    fprintf(out,"ymax_zoom = YX_max + delta\n");
    fprintf(out,"set yrange [ymin_zoom:ymax_zoom]\n");
    fprintf(out, "set y2tics 0.5,0.1,1\n");
        fprintf(out, "set y2r [0.8:1]\n");
        fprintf(out, "set ytics nomirror\n");
        fprintf(out, "set xtics nomirror\n");
        fprintf(out, "#\n");
        fprintf(out, "d(x)=dd\n");
        fprintf(out, "fit d(x) '%s' u :($12/1e9) via dd\n",nomdata);
        fprintf(out, "c(x)=cc\n");
    fprintf(out, "fit c(x) '%s' u :($15/1e9) via cc\n",nomdata);
    fprintf(out, "is(x)=is\n");
    fprintf(out, "fit is(x) '%s' u :($21) via is\n",nomdata);
    fprintf(out, "#\n");
    fprintf(out, "set title \"Résultat Homogéinisation \\n %s\\n %d cristaux - Calcul %s\"\n",param.get_nom("nommateriau").c_str(), tab[listeessai[nbessai][0]].nbcristaux,typecalcul);
    fprintf(out, "plot '%s' u 10 :($12/1e9) title 'GDH' pointtype 7 lc rgb \"blue\", '%s' u 10 :($15/1e9) title 'GCH' pointtype 7 lc rgb \"red\",'%s' u 10 : ($%d/1e9) title 'Borne Voigt' w l lc rgb \"brown\" lw 2,'%s' u 10 : ($%d/1e9) title 'Borne Reuss' w l lc rgb \"orange\" lw 2,  '%s'  u 10:21 axis x1y2 title 'Isotropie' lc rgb \"green\" pointtype 3,c(x) title sprintf(\"moyenne %%.1f\",cc) lc rgb \"red\" dt 2 ,d(x) title sprintf(\"moyenne %%.1f\",dd) lc rgb \"blue\" dt 2, is(x) title sprintf(\"moyenne %%.2f\",is) axis x1y2 lc rgb \"green\" dt 2\n", nomdata,nomdata, nomdata,39+9*tab[listeessai[nbessai][0]].nbphase, nomdata,40+9*tab[listeessai[nbessai][0]].nbphase,nomdata);
        fclose(out);
    sprintf(execute,"gnuplot %s  > %s.log 2>&1",nomout,nomout);
    system(execute);
        //----------------------------------------E
    sprintf(nomdata,"%s%d.txt",nom,nbessai+1);
    sprintf(nomout,"%s%dE.gnu",nom,nbessai+1);
    sprintf(nomoutsvg,"%s%dE.svg",nom,nbessai+1);
    out=fopen((char*)nomout,"wt");  
        fprintf(out, "set terminal svg\n");
        fprintf(out, "set output \"%s\"\n", nomoutsvg);
        fprintf(out, "set datafile separator \";\"\n");
        fprintf(out, "set key on  outside bottom center horizontal\n");
        fprintf(out, "set xlabel 'Numéro d''essai'\n");
        fprintf(out, "set ylabel 'Module d''Young E (GPa)'\n");
        fprintf(out, "set y2label 'Écart à l''isotropie'\n");
        fprintf(out, "set pointsize 1\n");
        fprintf(out, "set autoscale x\n");
        //fprintf(out, "set autoscale y\n");
        fprintf(out,"stats '%s' using ($%d/1e9) name 'YX'\n",nomdata,35+9*tab[listeessai[nbessai][0]].nbphase);
    fprintf(out,"stats '%s' using ($%d/1e9) name 'YN'\n",nomdata,36+9*tab[listeessai[nbessai][0]].nbphase);
    fprintf(out,"delta = (YX_max - YN_min) * 0.05\n");
    fprintf(out,"ymin_zoom = YN_min - delta\n");
    fprintf(out,"ymax_zoom = YX_max + delta\n");
    fprintf(out,"set yrange [ymin_zoom:ymax_zoom]\n");
    fprintf(out, "set y2tics 0.5,0.1,1\n");
        fprintf(out, "set y2r [0.8:1]\n");
        fprintf(out, "set ytics nomirror\n");
        fprintf(out, "set xtics nomirror\n");
        fprintf(out, "#\n");
        fprintf(out, "d(x)=dd\n");
        fprintf(out, "fit d(x) '%s' u :($13/1e9) via dd\n",nomdata);
        fprintf(out, "c(x)=cc\n");
    fprintf(out, "fit c(x) '%s' u :($16/1e9) via cc\n",nomdata);
    fprintf(out, "is(x)=is\n");
    fprintf(out, "fit is(x) '%s' u :($21) via is\n",nomdata);
    fprintf(out, "#\n");
    fprintf(out, "set title \"Résultat Homogéinisation \\n %s\\n %d cristaux - Calcul %s\"\n",param.get_nom("nommateriau").c_str(), tab[listeessai[nbessai][0]].nbcristaux,typecalcul);
    fprintf(out, "plot '%s' u 10 :($13/1e9) title 'EDH' pointtype 7 lc rgb \"blue\", '%s' u 10 :($16/1e9) title 'ECH' pointtype 7 lc rgb \"red\",'%s' u 10 : ($%d/1e9) title 'Borne Voigt' w l lc rgb \"brown\" lw 2,'%s' u 10 : ($%d/1e9) title 'Borne Reuss' w l lc rgb \"orange\" lw 2,  '%s'  u 10:21 axis x1y2 title 'Isotropie' lc rgb \"green\" pointtype 3,c(x) title sprintf(\"moyenne %%.1f\",cc) lc rgb \"red\" dt 2 ,d(x) title sprintf(\"moyenne %%.1f\",dd) lc rgb \"blue\" dt 2, is(x) title sprintf(\"moyenne %%.2f\",is) axis x1y2 lc rgb \"green\" dt 2\n", nomdata,nomdata, nomdata,35+9*tab[listeessai[nbessai][0]].nbphase, nomdata,36+9*tab[listeessai[nbessai][0]].nbphase,nomdata);
        fclose(out);
    sprintf(execute,"gnuplot %s  > %s.log 2>&1",nomout,nomout);
    system(execute);
        //----------------------------------------E cumulé
    sprintf(nomout,"%s%dcum.gnu",nom,nbessai+1);
    sprintf(nomoutsvg,"%s%dcum.svg",nom,nbessai+1);
    out=fopen((char*)nomout,"wt");  
        fprintf(out, "set terminal svg\n");
        fprintf(out, "set output \"%s\"\n", nomoutsvg);
        fprintf(out, "set multiplot layout 2,1\n");
        fprintf(out, "set datafile separator \";\"\n");
        fprintf(out, "set key on  outside bottom center horizontal\n");
        fprintf(out, "set ylabel 'Module d''Young E (GPa)'\n");
        fprintf(out, "set pointsize 1\n");
        fprintf(out, "set autoscale x\n");
        fprintf(out, "set autoscale y\n");
        fprintf(out, "set ytics nomirror\n");
        fprintf(out, "set xtics nomirror\n");
        fprintf(out, "#\n");
        fprintf(out, "d(x)=dd\n");
        fprintf(out, "fit d(x) '%s' u :($13/1e9) via dd\n",nomdata);
        fprintf(out, "c(x)=cc\n");
    fprintf(out, "fit c(x) '%s' u :($16/1e9) via cc\n",nomdata);
    fprintf(out, "is(x)=is\n");
    fprintf(out, "fit is(x) '%s' u :($21) via is\n",nomdata);
    fprintf(out, "#\n");
    fprintf(out, "set title \"Résultat Homogéinisation - Moyenne des n premiers échantillons\\n %s\\n %d cristaux - Calcul %s\"\n",param.get_nom("nommateriau").c_str(), tab[listeessai[nbessai][0]].nbcristaux,typecalcul);
    fprintf(out, "set size 1,0.70\n");
    fprintf(out, "set origin 0,0.3\n");
    fprintf(out, "plot '%s' u 10 :($18/1e9) title 'EDH' w l lc rgb \"blue\", '%s' u 10 :($19/1e9) title 'ECH' w l lc rgb \"red\",c(x) title sprintf(\"moyenne %%.1f\",cc) lc rgb \"red\" dt 2 ,d(x) title sprintf(\"moyenne %%.1f\",dd) lc rgb \"blue\" dt 2\n", nomdata,nomdata);
        fprintf(out, "unset title\n");
        fprintf(out, "set size 1,0.3\n");
    fprintf(out, "set origin 0,0\n");
    //fprintf(out, "set key outside top center horizontal\n");
    fprintf(out, "set xlabel 'Numéro d''essai'\n");
        fprintf(out, "set ytics 0.9,0.02,1\n");
        fprintf(out, "stats '%s' using 22 nooutput\n",nomdata);
    fprintf(out,"if (STATS_min < 0.96) {\n");
    fprintf(out,"    set yrange [*:1]\n");
    fprintf(out,"} else {\n");
    fprintf(out,"    set yrange [0.96:1]\n"); 
    fprintf(out,"}\n");
        //fprintf(out, "set autoscale y\n");
        fprintf(out, "set ylabel 'Écart à l''isotropie'\n");
        fprintf(out, "plot  '%s'  u 10:22 w l axis x1y1 title 'Isotropie' lc rgb \"green\",is(x) title sprintf(\"moyenne %%.2f\",is) axis x1y1 lc rgb \"green\" dt 2\n", nomdata);
        fprintf(out, "unset multiplot\n");
        fclose(out);
    sprintf(execute,"gnuplot %s  > %s.log 2>&1",nomout,nomout);
    system(execute);
        //----------------------------------------Volume masse cumulé
    sprintf(nomoutsvg,"%s%dmv.svg",nom,nbessai+1);
    sprintf(nomout,"%s%dmv.gnu",nom,nbessai+1);
    out=fopen((char*)nomout,"wt");  
        fprintf(out, "set terminal svg\n");
        fprintf(out, "set output \"%s\"\n", nomoutsvg);
        fprintf(out, "set datafile separator \";\"\n");
    fprintf(out, "set key on  outside bottom center horizontal\n");
        fprintf(out, "set pointsize 1\n");
        fprintf(out, "set autoscale x\n");
        fprintf(out, "set autoscale y\n");
    fprintf(out, "set autoscale y2\n");
    fprintf(out, "set y2tics\n");
        fprintf(out, "set ytics nomirror\n");
        fprintf(out, "set title \"Volume et masse des échantillons \\n%s\\n %d cristaux - Calcul %s\"\n",param.get_nom("nommateriau").c_str(), tab[listeessai[nbessai][0]].nbcristaux,typecalcul);
    fprintf(out, "set xlabel 'Numéro d''essai'\n");
        fprintf(out, "set ylabel 'Volume'\n");
        fprintf(out, "set y2label 'Masse'\n");
    fprintf(out, "d(x)=dd\n");
    fprintf(out, "fit d(x) '%s' u : 25 via dd\n",nomdata);
    int nbphase=(int)param.get_valeur("Nbphase");
    std::vector<double> deno(nbphase),cible(nbphase);
    double nume=1.;
    for (int i=0;i<nbphase;i++)
        deno[i]=1.;
    for (int i=0;i<nbphase;i++)
        {
        cible[i]=param.get_valeur("Proportionphase",i);    
        double densite=param.get_valeur("Densite",i);    
        nume=nume*densite;
        for (int j=0;j<nbphase;j++)
            if (i==j) deno[j]=deno[j]*cible[i]; else deno[j]=deno[j]*densite;
        }
    double masse=0.;
    for (int i=0;i<nbphase;i++)
        masse=masse+deno[i];
    masse=nume/masse;
    fprintf(out,"bb=%lf\n",masse);
    fprintf(out,"b(x)=bb\n");
    fprintf(out, "plot '%s' u 10 : 24 title 'Volume' pointtype 6 lc rgb \"blue\", '%s' u 10 : 25 title 'Masse' axis x1y2 pointtype 6 lc rgb \"red\", d(x) title sprintf(\"moyenne masse %%.6f\",dd) axis x1y2 lc rgb \"red\" dt 2,b(x) title sprintf(\"masse cible %%.6f\",bb) axis x1y2 lc rgb \"green\" dt 2\n", nomdata,nomdata);
    fclose(out);
    sprintf(execute,"gnuplot %s  > %s.log 2>&1",nomout,nomout);
    system(execute);
        //----------------------------------------Centre de masse
    sprintf(nomoutsvg,"%s%dcdm.svg",nom,nbessai+1);
    sprintf(nomout,"%s%dcdm.gnu",nom,nbessai+1);
    out=fopen((char*)nomout,"wt");  
        fprintf(out, "set terminal svg\n");
        fprintf(out, "set output \"%s\"\n", nomoutsvg);
        fprintf(out, "set datafile separator \";\"\n");
    fprintf(out, "set key on  outside bottom center horizontal\n");
        fprintf(out, "set pointsize 1\n");
        fprintf(out, "set autoscale x\n");
        fprintf(out, "set autoscale y\n");
    fprintf(out, "set title \"Centre de masse des échantillons \\n %s\\n %d cristaux - Calcul %s\"\n",param.get_nom("nommateriau").c_str(), tab[listeessai[nbessai][0]].nbcristaux,typecalcul);
    fprintf(out, "set xlabel 'Numéro d''essai'\n");
        fprintf(out, "set ylabel 'Volume'\n");
    fprintf(out, "plot '%s' u 10 : 26 title 'CDMx' pointtype 6 lc rgb \"blue\", '%s' u 10 : 27 title 'CDMy' pointtype 6 lc rgb \"red\", '%s' u 10 : 28 title 'CDMz' pointtype 6 lc rgb \"green\", %lf notitle dt 2  lc rgb \"blue\", %lf notitle dt 2  lc rgb \"red\", %lf notitle dt 2  lc rgb \"green\"\n",nomdata, nomdata,nomdata,xcible,ycible,zcible);
    fclose(out);
    sprintf(execute,"gnuplot %s  > %s.log 2>&1",nomout,nomout);
    system(execute);
        //----------------------------------------Centre de masse cumulé
    sprintf(nomoutsvg,"%s%dcdmcum.svg",nom,nbessai+1);
    sprintf(nomout,"%s%dcdmcum.gnu",nom,nbessai+1);
    out=fopen((char*)nomout,"wt");  
        fprintf(out, "set terminal svg\n");
        fprintf(out, "set output \"%s\"\n", nomoutsvg);
        fprintf(out, "set datafile separator \";\"\n");
    fprintf(out, "set key on  outside bottom center horizontal\n");
        fprintf(out, "set pointsize 1\n");
        fprintf(out, "set autoscale x\n");
        fprintf(out, "set autoscale y\n");
    fprintf(out, "set title \"Position moyenne des CDM des n premiers échantillons\\n %s\\n %d cristaux - Calcul %s\"\n",param.get_nom("nommateriau").c_str(), tab[listeessai[nbessai][0]].nbcristaux,typecalcul);
    fprintf(out, "set xlabel 'Numéro d''essai'\n");
        fprintf(out, "set ylabel 'Volume'\n");
    fprintf(out, "plot '%s' u 10 : 29 title 'CDMx' w l lc rgb \"blue\", '%s' u 10 : 30 title 'CDMy' w l lc rgb \"red\", '%s' u 10 : 31 title 'CDMz' w l lc rgb \"green\", %lf notitle dt 2  lc rgb \"blue\", %lf notitle dt 2  lc rgb \"red\", %lf notitle dt 2  lc rgb \"green\"\n",nomdata, nomdata,nomdata,xcible,ycible,zcible);
    fclose(out); 
    sprintf(execute,"gnuplot %s  > %s.log 2>&1",nomout,nomout);
    system(execute);
    //----------------------------------------Par phase
    if (tab[listeessai[nbessai][0]].nbphase>1)
      for (int j=0;j<tab[listeessai[nbessai][0]].nbphase;j++)
        {
        sprintf(nomoutsvg,"%s%dph%dmv.svg",nom,nbessai+1,j+1);
        sprintf(nomout,"%s%dph%dmv.gnu",nom,nbessai+1,j+1);
        out=fopen((char*)nomout,"wt");  
        fprintf(out, "set terminal svg \n");
        fprintf(out, "set output \"%s\"\n", nomoutsvg);
        fprintf(out, "set datafile separator \";\"\n");
        fprintf(out, "set key on  outside bottom center horizontal\n");
        fprintf(out, "set pointsize 1\n");
        fprintf(out, "set autoscale x\n");
        fprintf(out, "set autoscale y\n");
        fprintf(out, "set autoscale y2\n");
        fprintf(out, "set y2tics\n");
        fprintf(out, "set ytics nomirror\n");
        fprintf(out, "set title \"Volume et masse de la phase %s des échantillons \\n%s\\n %d cristaux - Calcul %s\"\n",param.get_nom("Nomphase",j).c_str(),param.get_nom("nommateriau").c_str(), tab[listeessai[nbessai][0]].nbcristaux,typecalcul);
        fprintf(out, "set xlabel 'Numéro d''essai'\n");
        fprintf(out, "set ylabel 'Volume'\n");
        fprintf(out, "set y2label 'Masse'\n");
        fprintf(out, "d(x)=dd\n");
        fprintf(out, "fit d(x) '%s' u : %d via dd\n",nomdata,34+9*j);
        fprintf(out,"bb=%lf\n",cible[j]*masse);
        fprintf(out,"b(x)=%lf\n",cible[j]*masse);
        fprintf(out, "plot '%s' u 10 : %d title 'Volume' pointtype 6 lc rgb \"blue\", '%s' u 10 : %d title 'Masse' axis x1y2 pointtype 6 lc rgb \"red\", d(x) title sprintf(\"moyenne masse %%.6f\",dd) axis x1y2 lc rgb \"red\" dt 2,b(x) title sprintf(\"masse cible %%.6f\",bb) axis x1y2 lc rgb \"green\" dt 2\n", nomdata,33+9*j,nomdata,34+9*j);
        fclose(out);
        sprintf(execute,"gnuplot %s  > %s.log 2>&1",nomout,nomout);
        system(execute);
        sprintf(nomoutsvg,"%s%dph%dcdm.svg",nom,nbessai+1,j+1);
        sprintf(nomout,"%s%dph%dcdm.gnu",nom,nbessai+1,j+1);
        out=fopen((char*)nomout,"wt");  
        fprintf(out, "set terminal svg\n");
        fprintf(out, "set output \"%s\"\n", nomoutsvg);
        fprintf(out, "set datafile separator \";\"\n");
        fprintf(out, "set key on  outside bottom center horizontal\n");
        fprintf(out, "set pointsize 1\n");
        fprintf(out, "set autoscale x\n");
        fprintf(out, "set autoscale y\n");
        fprintf(out, "set title \"Centre de masse de la phase %s des échantillons \\n %s\\n %d cristaux - Calcul %s\"\n",param.get_nom("Nomphase",j).c_str(),param.get_nom("nommateriau").c_str(), tab[listeessai[nbessai][0]].nbcristaux,typecalcul);
        fprintf(out, "set xlabel 'Numéro d''essai'\n");
        fprintf(out, "set ylabel 'Volume'\n");
        fprintf(out, "plot '%s' u 10 : %d title 'CDMx' pointtype 6 lc rgb \"blue\", '%s' u 10 : %d title 'CDMy' pointtype 6 lc rgb \"red\", '%s' u 10 : %d title 'CDMz' pointtype 6 lc rgb \"green\",  %lf notitle dt 2  lc rgb \"blue\", %lf notitle dt 2  lc rgb \"red\", %lf notitle dt 2  lc rgb \"green\"\n",nomdata,35+9*j,nomdata,36+9*j,nomdata,37+9*j,xcible,ycible,zcible);
        fclose(out);
        sprintf(execute,"gnuplot %s  > %s.log 2>&1",nomout,nomout);
        system(execute);
        sprintf(nomoutsvg,"%s%dph%dcdmcum.svg",nom,nbessai+1,j+1);
        sprintf(nomout,"%s%dph%dcdmcum.gnu",nom,nbessai+1,j+1);
        out=fopen((char*)nomout,"wt");  
        fprintf(out, "set terminal svg\n");
        fprintf(out, "set output \"%s\"\n", nomoutsvg);
        fprintf(out, "set datafile separator \";\"\n");
        fprintf(out, "set key on  outside bottom center horizontal\n");
        fprintf(out, "set pointsize 1\n");
        fprintf(out, "set autoscale x\n");
        fprintf(out, "set autoscale y\n");
        fprintf(out, "set title \"Position moyenne des CDM de la phase %s des n premiers échantillons\\n %s\\n %d cristaux - Calcul %s\"\n",param.get_nom("Nomphase",j).c_str(),param.get_nom("nommateriau").c_str(), tab[listeessai[nbessai][0]].nbcristaux,typecalcul);
        fprintf(out, "set xlabel 'Numéro d''essai'\n");
        fprintf(out, "set ylabel 'Volume'\n");
        fprintf(out, "plot '%s' u 10 : %d title 'CDMx' w l lc rgb \"blue\", '%s' u 10 : %d title 'CDMy' w l lc rgb \"red\", '%s' u 10 : %d title 'CDMz' w l lc rgb \"green\", %lf notitle dt 2  lc rgb \"blue\", %lf notitle dt 2  lc rgb \"red\", %lf notitle dt 2  lc rgb \"green\"\n",nomdata,38+9*j, nomdata,39+9*j,nomdata,40+9*j,xcible,ycible,zcible);
        fclose(out);  
        sprintf(execute,"gnuplot %s  > %s.log 2>&1",nomout,nomout);
        system(execute);
        }
        }
 
int MSTRUCT_POLY_ANALYSE_RESULTAT::get_nbechantillon(int nbessai)
{
return listeessai[nbessai].size();
    
}
 
 
int MSTRUCT_POLY_ANALYSE_RESULTAT::get_degreechantillon(int nbessai)
{
return tab[listeessai[nbessai][0]].degre;
}
 
int MSTRUCT_POLY_ANALYSE_RESULTAT::get_nbcristauxechantillon(int nbessai)
{
return tab[listeessai[nbessai][0]].nbcristaux;
}
 
int MSTRUCT_POLY_ANALYSE_RESULTAT::get_nbechantillonratee(void)
{
return listeessairatee.size();
    
}