compare/src/visumain.cpp

//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//  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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       visumain.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:57:20 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------

#include <iostream>
#include "vtkdisplay.hxx"
#include "ndata.hxx"
#include "mg_file.h"
#include <map>
#if _COMPARESTEP_
#include "occ_import.h"
#endif

void affiche(char* message)
{

    std::cout << message << std::endl;

}


struct groupe_couleur
{
  int num;
  std::map<unsigned long,unsigned long> A;
  std::map<unsigned long,unsigned long> B;
  bool est_dans_A(unsigned long val)
      {
      std::map<unsigned long,unsigned long>::iterator it=A.find(val);
      if (it!=A.end()) return true;
      return false;
      };
  bool est_dans_B(unsigned long val)
      {
      std::map<unsigned long,unsigned long>::iterator it=B.find(val);
      if (it!=B.end()) return true;
      return false;
      };
};


struct data_fichier
{
std::string nomA;  
std::string nomB;  
std::map<int,groupe_couleur> list_face_similaire;
std::map<int,groupe_couleur> list_face_identique;
std::map<int,groupe_couleur> list_face_localise;
std::map<int,groupe_couleur> list_arete_localise;
std::map<int,groupe_couleur> list_sommet_localise;
std::map<int,groupe_couleur> list_face_modifie;
std::map<int,groupe_couleur> list_arete_modifie;
};


void ini_couleur(std::vector<color> &tabcoul,int dim)
{
 int ecart=255/(dim+1);
 std::map<int,int> list;
 for (int i=0;i<dim*dim*dim;i++)
    {
    int r=rand()%dim+1;
    int g=rand()%dim+1;
    int b=rand()%dim+1;
    int c=r*ecart+255*g*ecart+255*255*b*ecart;
    std::pair<int,int> tmp(c,c);
    std::pair<std::map<int,int>::iterator,bool> ret=list.insert(tmp);
    if (ret.second==false) i--;
    else 
      {
      color cc(r*ecart,g*ecart,b*ecart);
      tabcoul.push_back(cc);  
      }
    } 
}


void ajouter_liste_couleur(std::map<int,groupe_couleur> &list,unsigned long valA,unsigned long valB)
{
static int num_couleur=0;
std::map<int,groupe_couleur>::iterator it;
for (it=list.begin();it!=list.end();it++)
    {
    bool res=it->second.est_dans_A(valA);
    if (res==true) 
        {
          it->second.A.insert(std::pair<unsigned long,unsigned long>(valA,valA));
          it->second.B.insert(std::pair<unsigned long,unsigned long>(valB,valB));
          return;
        }
          
    res=it->second.est_dans_B(valB);
    if (res==true) 
        {
          it->second.A.insert(std::pair<unsigned long,unsigned long>(valA,valA));
          it->second.B.insert(std::pair<unsigned long,unsigned long>(valB,valB));
          return;
        }
          
    }
groupe_couleur gc;
std::pair<std::map<int,groupe_couleur>::iterator,bool> ret=list.insert(std::pair<int,groupe_couleur>(num_couleur,gc));  
ret.first->second.num=num_couleur;
ret.first->second.A.insert(std::pair<unsigned long,unsigned long>(valA,valA));
ret.first->second.B.insert(std::pair<unsigned long,unsigned long>(valB,valB));
num_couleur++;
}


void lire_fichier(char *nom,data_fichier &data)
{
FILE* in=fopen(nom,"rt");

while (!feof(in))
{
  char mess[255],mess2[255];
  fgets(mess,255,in);
  if (strlen(mess)>1) 
    if (mess[0]=='/') 
      if (mess[1]=='/') continue;
  if (strcmp(mess,"Model A:\n")==0)
    {
    fgets(mess,255,in);
    sscanf(mess,"%s",mess2);    
    data.nomA=mess2;
    }
  else if (strcmp(mess,"Model B:\n")==0)
    {
    fgets(mess,255,in);
    sscanf(mess,"%s",mess2);    
    data.nomB=mess2;
    }
  else
    {
    char mot1[255],mot2[255],mot3[255];
    int nb;
    int n=sscanf(mess,"%d %s %s %s",&nb,mot1,mot2,mot3);
    std::map<int,groupe_couleur> *list=NULL;
    if (n>2)
      {
        if (strcmp(mot1,"Similar")==0)
        if (strcmp(mot2,"faces")==0)
          list= &(data.list_face_similaire);
        if (strcmp(mot1,"Identical")==0)
        if (strcmp(mot2,"faces")==0)
          list= &(data.list_face_identique);       
        if (strcmp(mot1,"Localized")==0)
        if (strcmp(mot2,"identical")==0)
        if (strcmp(mot3,"faces")==0)
          list=&(data.list_face_localise);
        if (strcmp(mot1,"Localized")==0)
        if (strcmp(mot2,"identical")==0)
        if (strcmp(mot3,"edges")==0)
          list=&(data.list_arete_localise);
        if (strcmp(mot1,"Localized")==0)
        if (strcmp(mot2,"identical")==0)
        if (strcmp(mot3,"vertices")==0)
          list=&(data.list_sommet_localise);
        if (strcmp(mot1,"Localized")==0)
        if (strcmp(mot2,"modified")==0)
        if (strcmp(mot3,"faces")==0)
          list=&(data.list_face_modifie);
        if (strcmp(mot1,"Localized")==0)
        if (strcmp(mot2,"modified")==0)
        if (strcmp(mot3,"edges")==0)
          list=&(data.list_arete_modifie);
    if (list!=NULL)  
      for (int i=0;i<nb;i++)
        {
        fgets(mess,255,in);
        unsigned long valA,valB;
        sscanf(mess,"%lu %lu",&valA,&valB);     
        ajouter_liste_couleur(*list,valA,valB);
        }
      }
    }

  
}


fclose(in);
}


void affiche_structure(color base,MG_GEOMETRIE* geo,bool A,std::vector<color> &tabcoul,data_container &data, std::map<int,groupe_couleur> &list_face, std::map<int,groupe_couleur> &list_arete, std::map<int,groupe_couleur> &list_sommet)
{
    unsigned long nb1=geo->get_mg_face(geo->get_nb_mg_face()-1)->get_id();
    unsigned long nb2=geo->get_mg_arete(geo->get_nb_mg_arete()-1)->get_id();
    unsigned long nb3=geo->get_mg_sommet(geo->get_nb_mg_sommet()-1)->get_id();
    int nbid=nb1;
    if (nb2>nbid) nbid=nb2;
    if (nb3>nbid) nbid=nb3;
    nbid=nbid+5;
    bool affiche[nbid];
    for (int i=0;i<nbid;i++) 
      affiche[i]=false;
    std::map<int,groupe_couleur>::iterator it;
    for (it=list_face.begin();it!=list_face.end();it++)
    {
    groupe_couleur gc=it->second;
    properties ptr=data.getproptriangles();
    ptr.c=color(tabcoul[gc.num].R,tabcoul[gc.num].G,tabcoul[gc.num].B,255);
    ptr.edgeon=false;
    data.setproptriangles(ptr);
    std::map<unsigned long,unsigned long>::iterator it2;
    std::map<unsigned long,unsigned long> *tmp;
    if (A==true) tmp=&(gc.A);
      else tmp=&(gc.B);
    
    for (it2=tmp->begin();it2!=tmp->end();it2++)
      {
      MG_FACE* face=geo->get_mg_faceid(it2->second);
      affiche[face->get_id()]=true;
      double xg=0.,yg=0.,zg=0.,deno=0.;
      int nbtri=face->get_lien_maillage()->get_nb();
      for (int i=0;i<nbtri;i++)
          {
          MG_TRIANGLE* tri=(MG_TRIANGLE*)face->get_lien_maillage()->get(i);
          double *xyz1=tri->get_noeud1()->get_coord();
          double *xyz2=tri->get_noeud2()->get_coord();
          double *xyz3=tri->get_noeud3()->get_coord();
          npoint p1(xyz1[0],xyz1[1],xyz1[2]);
          npoint p2(xyz2[0],xyz2[1],xyz2[2]);
          npoint p3(xyz3[0],xyz3[1],xyz3[2]);
          triangle tr;
          tr.pts[0]=p1;
          tr.pts[1]=p2;
          tr.pts[2]=p3;
          data.add_triangle(tr);
          OT_VECTEUR_3D vec1(xyz1,xyz2);
          OT_VECTEUR_3D vec2(xyz1,xyz3);
          OT_VECTEUR_3D sur=vec1&vec2;
          xg=xg+0.3333333333*(xyz1[0]+xyz2[0]+xyz3[0])*sur.get_longueur();
          yg=yg+0.3333333333*(xyz1[1]+xyz2[1]+xyz3[1])*sur.get_longueur();
          zg=zg+0.3333333333*(xyz1[2]+xyz2[2]+xyz3[2])*sur.get_longueur();
          deno=deno+sur.get_longueur();
          }
      xg=xg/deno;
      yg=yg/deno;
      zg=zg/deno;
      std::map<double,MG_TRIANGLE*> classetri;
      for (int i=0;i<nbtri;i++)
          {
          MG_TRIANGLE* tri=(MG_TRIANGLE*)face->get_lien_maillage()->get(i);
          double *xyz1=tri->get_noeud1()->get_coord();
          double *xyz2=tri->get_noeud2()->get_coord();
          double *xyz3=tri->get_noeud3()->get_coord();  
          double x=0.3333333333*(xyz1[0]+xyz2[0]+xyz3[0]);
          double y=0.3333333333*(xyz1[1]+xyz2[1]+xyz3[1]);
          double z=0.3333333333*(xyz1[2]+xyz2[2]+xyz3[2]);
          OT_VECTEUR_3D vec(x-xg,y-yg,z-zg);
          classetri.insert(std::pair<double,MG_TRIANGLE*>(vec.get_longueur2(),tri)); 
          }
      if (nbtri>0)
        {
        MG_TRIANGLE *tri=classetri.begin()->second;
        double *xyz1=tri->get_noeud1()->get_coord();
        double *xyz2=tri->get_noeud2()->get_coord();
        double *xyz3=tri->get_noeud3()->get_coord();  
        npoint p(0.3333333333*(xyz1[0]+xyz2[0]+xyz3[0]),0.3333333333*(xyz1[1]+xyz2[1]+xyz3[1]),0.3333333333*(xyz1[2]+xyz2[2]+xyz3[2]));
        point pp;
        pp.pts=p;
        char mess[255];
        sprintf(mess,"F%lu",face->get_id());
        pp.info=mess;   
        color cc=color(255-ptr.c.R,255-ptr.c.G,255-ptr.c.B);
        std::pair<point,color> tmp(pp,cc);
        data.add_text(2,tmp);    
        }
          
        }
    }
    properties ptr=data.getproptriangles();
    ptr.c=base;
    ptr.edgeon=false;
    data.setproptriangles(ptr);
    LISTE_MG_FACE::iterator itface;
     for (MG_FACE* face=geo->get_premier_face(itface);face!=NULL;face=geo->get_suivant_face(itface))
     {
     if (affiche[face->get_id()]==true)  continue;
     affiche[face->get_id()]=true;
     int nbtri=face->get_lien_maillage()->get_nb();
     double xg=0.,yg=0.,zg=0.,deno=0.;
      for (int i=0;i<nbtri;i++)
          {
          MG_TRIANGLE* tri=(MG_TRIANGLE*)face->get_lien_maillage()->get(i);
          double *xyz1=tri->get_noeud1()->get_coord();
          double *xyz2=tri->get_noeud2()->get_coord();
          double *xyz3=tri->get_noeud3()->get_coord();
          npoint p1(xyz1[0],xyz1[1],xyz1[2]);
          npoint p2(xyz2[0],xyz2[1],xyz2[2]);
          npoint p3(xyz3[0],xyz3[1],xyz3[2]);
          triangle tr;
          tr.pts[0]=p1;
          tr.pts[1]=p2;
          tr.pts[2]=p3;
          data.add_triangle(tr);
          OT_VECTEUR_3D vec1(xyz1,xyz2);
          OT_VECTEUR_3D vec2(xyz1,xyz3);
          OT_VECTEUR_3D sur=vec1&vec2;
          xg=xg+0.3333333333*(xyz1[0]+xyz2[0]+xyz3[0])*sur.get_longueur();
          yg=yg+0.3333333333*(xyz1[1]+xyz2[1]+xyz3[1])*sur.get_longueur();
          zg=zg+0.3333333333*(xyz1[2]+xyz2[2]+xyz3[2])*sur.get_longueur();
          deno=deno+sur.get_longueur();
          }
    xg=xg/deno;
    yg=yg/deno;
    zg=zg/deno;
    std::map<double,MG_TRIANGLE*> classetri;
    for (int i=0;i<nbtri;i++)
          {
          MG_TRIANGLE* tri=(MG_TRIANGLE*)face->get_lien_maillage()->get(i);
          double *xyz1=tri->get_noeud1()->get_coord();
          double *xyz2=tri->get_noeud2()->get_coord();
          double *xyz3=tri->get_noeud3()->get_coord();  
          double x=0.3333333333*(xyz1[0]+xyz2[0]+xyz3[0]);
          double y=0.3333333333*(xyz1[1]+xyz2[1]+xyz3[1]);
          double z=0.3333333333*(xyz1[2]+xyz2[2]+xyz3[2]);
          OT_VECTEUR_3D vec(x-xg,y-yg,z-zg);
          classetri.insert(std::pair<double,MG_TRIANGLE*>(vec.get_longueur2(),tri)); 
          }
    if (nbtri>0)
      {
      MG_TRIANGLE *tri=classetri.begin()->second;
      double *xyz1=tri->get_noeud1()->get_coord();
      double *xyz2=tri->get_noeud2()->get_coord();
      double *xyz3=tri->get_noeud3()->get_coord();  
      npoint p(0.3333333333*(xyz1[0]+xyz2[0]+xyz3[0]),0.3333333333*(xyz1[1]+xyz2[1]+xyz3[1]),0.3333333333*(xyz1[2]+xyz2[2]+xyz3[2]));
      point pp;
      pp.pts=p;
      char mess[255];
      sprintf(mess,"F%lu",face->get_id());
      pp.info=mess;     
      color cc=color(255-ptr.c.R,255-ptr.c.G,255-ptr.c.B);
      std::pair<point,color> tmp(pp,cc);
      data.add_text(2,tmp);    
      }
    }
    
      
    for (it=list_arete.begin();it!=list_arete.end();it++)
    {
    groupe_couleur gc=it->second;
    properties ptr=data.getproplines();
    ptr.c=tabcoul[gc.num];
    ptr.thickness=5;
    data.setproplines(ptr);
    std::map<unsigned long,unsigned long>::iterator it2;
    std::map<unsigned long,unsigned long> *tmp;
    if (A==true) tmp=&(gc.A);
      else tmp=&(gc.B);
    
    for (it2=tmp->begin();it2!=tmp->end();it2++)
      {
      MG_ARETE* are=geo->get_mg_areteid(it2->second);
      affiche[are->get_id()]=true;
      int nbseg=are->get_lien_maillage()->get_nb();
      double xg=0.,yg=0.,zg=0.,deno=0.;
      for (int i=0;i<nbseg;i++)
          {
          MG_SEGMENT* seg=(MG_SEGMENT*)are->get_lien_maillage()->get(i);
          double *xyz1=seg->get_noeud1()->get_coord();
          double *xyz2=seg->get_noeud2()->get_coord();
          npoint p1(xyz1[0],xyz1[1],xyz1[2]);
          npoint p2(xyz2[0],xyz2[1],xyz2[2]);
          line li;
          li.pts[0]=p1;
          li.pts[1]=p2;
          data.add_line(li);
          OT_VECTEUR_3D vec(xyz1,xyz2);
          xg=xg+0.5*(xyz1[0]+xyz2[0])*vec.get_longueur();
          yg=yg+0.5*(xyz1[1]+xyz2[1])*vec.get_longueur();
          zg=zg+0.5*(xyz1[2]+xyz2[2])*vec.get_longueur();
          deno=deno+vec.get_longueur();
          }
      xg=xg/deno;
      yg=yg/deno;
      zg=zg/deno;
      std::map<double,MG_SEGMENT*> classeseg;
      for (int i=0;i<nbseg;i++)
          {
           MG_SEGMENT* seg=(MG_SEGMENT*)are->get_lien_maillage()->get(i);
            double *xyz1=seg->get_noeud1()->get_coord();
            double *xyz2=seg->get_noeud2()->get_coord();
            double x=0.5*(xyz1[0]+xyz2[0]);
            double y=0.5*(xyz1[1]+xyz2[1]);
            double z=0.5*(xyz1[2]+xyz2[2]);
            OT_VECTEUR_3D vec(x-xg,y-yg,z-zg);
          classeseg.insert(std::pair<double,MG_SEGMENT*>(vec.get_longueur2(),seg)); 
          }
      if (nbseg>0)
        {
        MG_SEGMENT* seg=classeseg.begin()->second;
        double *xyz1=seg->get_noeud1()->get_coord();
        double *xyz2=seg->get_noeud2()->get_coord();
        npoint p(0.5*(xyz1[0]+xyz2[0]),0.5*(xyz1[1]+xyz2[1]),0.5*(xyz1[2]+xyz2[2]));
        point pp;
        pp.pts=p;
        char mess[255];
        sprintf(mess,"E%lu",are->get_id());
        pp.info=mess;
        color cc=color(ptr.c.R,ptr.c.G,ptr.c.B);
        std::pair<point,color> tmp(pp,cc);
        data.add_text(1,tmp);  
        }
      }
    }
    ptr=data.getproplines();
    ptr.c=color(255-base.R,255-base.G,255-base.B);
    ptr.thickness=2;
    data.setproplines(ptr);
    LISTE_MG_ARETE::iterator itare;
     for (MG_ARETE* are=geo->get_premier_arete(itare);are!=NULL;are=geo->get_suivant_arete(itare))
     {
     if (affiche[are->get_id()]==true)  continue;
     affiche[are->get_id()]=true;
     int nbseg=are->get_lien_maillage()->get_nb();
     double xg=0.,yg=0.,zg=0.,deno=0.;
     for (int i=0;i<nbseg;i++)
          {
          MG_SEGMENT* seg=(MG_SEGMENT*)are->get_lien_maillage()->get(i);
          double *xyz1=seg->get_noeud1()->get_coord();
          double *xyz2=seg->get_noeud2()->get_coord();
          npoint p1(xyz1[0],xyz1[1],xyz1[2]);
          npoint p2(xyz2[0],xyz2[1],xyz2[2]);
          line li;
          li.pts[0]=p1;
          li.pts[1]=p2;
          data.add_line(li);
          OT_VECTEUR_3D vec(xyz1,xyz2);
          xg=xg+0.5*(xyz1[0]+xyz2[0])*vec.get_longueur();
          yg=yg+0.5*(xyz1[1]+xyz2[1])*vec.get_longueur();
          zg=zg+0.5*(xyz1[2]+xyz2[2])*vec.get_longueur();
          deno=deno+vec.get_longueur();
          }
     xg=xg/deno;
     yg=yg/deno;
     zg=zg/deno;
     std::map<double,MG_SEGMENT*> classeseg;
     for (int i=0;i<nbseg;i++)
          {
           MG_SEGMENT* seg=(MG_SEGMENT*)are->get_lien_maillage()->get(i);
            double *xyz1=seg->get_noeud1()->get_coord();
            double *xyz2=seg->get_noeud2()->get_coord();
            double x=0.5*(xyz1[0]+xyz2[0]);
            double y=0.5*(xyz1[1]+xyz2[1]);
            double z=0.5*(xyz1[2]+xyz2[2]);
            OT_VECTEUR_3D vec(x-xg,y-yg,z-zg);
          classeseg.insert(std::pair<double,MG_SEGMENT*>(vec.get_longueur2(),seg)); 
          }
      if (nbseg>0)
        {
        MG_SEGMENT* seg=classeseg.begin()->second;
        double *xyz1=seg->get_noeud1()->get_coord();
        double *xyz2=seg->get_noeud2()->get_coord();
        npoint p(0.5*(xyz1[0]+xyz2[0]),0.5*(xyz1[1]+xyz2[1]),0.5*(xyz1[2]+xyz2[2]));
        point pp;
        pp.pts=p;
        char mess[255];
        sprintf(mess,"E%lu",are->get_id());
        pp.info=mess;
        color cc=color(ptr.c.R,ptr.c.G,ptr.c.B);
        std::pair<point,color> tmp(pp,cc);
        data.add_text(1,tmp);   
        }
    }
     
    
    for (it=list_sommet.begin();it!=list_sommet.end();it++)
    {
    groupe_couleur gc=it->second;
    properties ptr=data.getproppoints();
    ptr.c=tabcoul[gc.num];
    ptr.pointsize=15;
    data.setproppoints(ptr);
    std::map<unsigned long,unsigned long>::iterator it2;
    std::map<unsigned long,unsigned long> *tmp;
    if (A==true) tmp=&(gc.A);
      else tmp=&(gc.B);
    
    for (it2=tmp->begin();it2!=tmp->end();it2++)
      {
      MG_SOMMET* som=geo->get_mg_sommetid(it2->second);
      affiche[som->get_id()]=true;
      int nbpoint=som->get_lien_maillage()->get_nb();
      for (int i=0;i<nbpoint;i++)
          {
          MG_NOEUD* nd=(MG_NOEUD*)som->get_lien_maillage()->get(i);
          double *xyz1=nd->get_coord();
          npoint p1(xyz1[0],xyz1[1],xyz1[2]);
          data.add_point(p1);
          point pp;
          pp.pts=p1;
          char mess[255];
          sprintf(mess,"V%lu",som->get_id());
          pp.info=mess;
          color cc=color(ptr.c.R,ptr.c.G,ptr.c.B);
          std::pair<point,color> tmp(pp,cc);
          data.add_text(0,tmp);
          }
      }
    }
    ptr=data.getproppoints();
    ptr.c=color(255-base.R,255-base.G,255-base.B);
    ptr.pointsize=4;
    
    data.setproppoints(ptr);
    LISTE_MG_SOMMET::iterator itsom;
     for (MG_SOMMET* som=geo->get_premier_sommet(itsom);som!=NULL;som=geo->get_suivant_sommet(itsom))
     {
     if (affiche[som->get_id()]==true)  continue;
     affiche[som->get_id()]=true;
     int nbpoint=som->get_lien_maillage()->get_nb();
      for (int i=0;i<nbpoint;i++)
          {
          MG_NOEUD* nd=(MG_NOEUD*)som->get_lien_maillage()->get(i);
          double *xyz1=nd->get_coord();
          npoint p1(xyz1[0],xyz1[1],xyz1[2]);
          data.add_point(p1);
          point pp;
          pp.pts=p1;
          char mess[255];
          sprintf(mess,"V%lu",som->get_id());
          pp.info=mess;
          color cc=color(ptr.c.R,ptr.c.G,ptr.c.B);
          std::pair<point,color> tmp(pp,cc);
          data.add_text(0,tmp);  
          }
      
    }
      
   
};


int main(int argc,char **argv)
{       
    color base=color(0,0,0);
    char chemin[1000];
    sprintf(chemin,"%s/.vtkdisplay",getenv("HOME"));
    FILE *in=fopen(chemin,"rt");
    if (in!=NULL)
        {
        char mess[255];
        fgets(mess,255,in);  
        int R,G,B;
        sscanf(mess,"%d %d %d",&R,&G,&B);
        base=color(R,G,B);
        fclose(in);
        }
    printf("******************************\n");
    printf("*   Comparison visualization *\n");
    printf("*            by              *\n");
    printf("*  Jean-Christophe Cuillière *\n");
    printf("*            and             *\n");
    printf("*       Vincent Francois     *\n");
    printf("*          ERICCA-UQTR       *\n");
    printf("******************************\n\n\n");
    GESTIONVERSION v;
    char version[500];
    v.print(version);
    affiche(version);
    if (argc!=4) 
        {
        printf("Visucompare syntax error - Bad number arguments \n visucompare.exe A|B S|M|I|L nom_fichier_resultat_comparaison\n\n");
        return 0;
          
        }
    if (!((strcmp(argv[1],"A")==0) || (strcmp(argv[1],"B")==0)))
        {
        printf("Visucompare syntax error - first argument not valid : %s\n\n",argv[1]);
        return 0;
          
        }
    
    if (!( (strcmp(argv[2],"S")==0) || (strcmp(argv[2],"M")==0) || (strcmp(argv[2],"I")==0)|| (strcmp(argv[2],"L")==0) ))
        {
        printf("Visucompare syntax error - Second argument not valid : %s\n\n",argv[2]);
        return 0;
          
        }
    
    
    std::map<int,groupe_couleur> list_vide;
    std::vector<color> tabcoul;
    ini_couleur(tabcoul,10);
    
    data_fichier df;
    lire_fichier(argv[3],df);
    bool A;
    #if _COMPARESTEP_
      MG_GESTIONNAIRE gest=NULL;
      if (strcmp(argv[1],"A")==0) 
        {
        printf("Reading file A\n");
        OCC_IMPORT occimport;
        MG_GEOMETRIE* geo1=occimport.importer(gest,(char*)df.nomA.c_str(),FICHIERSTEP,1.,1e-6);
        A=true;
        }
      else if (strcmp(argv[1],"B")==0) 
        {
        printf("Reading file B\n");
        OCC_IMPORT occimport;
        MG_GEOMETRIE* geo2=occimport.importer(gest,(char*)df.nomB.c_str(),FICHIERSTEP,1.,1e-6);
        A=false;
        } 
    #else
      MG_FILE *gest=NULL;
      if (strcmp(argv[1],"A")==0) {printf("Reading file A\n");gest=new MG_FILE((char*)df.nomA.c_str());A=true;}
        else if (strcmp(argv[1],"B")==0) {printf("Reading file B\n");gest=new MG_FILE((char*)df.nomB.c_str());A=false;}
    #endif
    printf("  Number of color group :\n");
    if (strcmp(argv[2],"S")==0) printf("      Similar faces      : %d Visible\n",(int)df.list_face_similaire.size()); else printf("      Similar faces      : %d\n",(int)df.list_face_similaire.size());
    if (strcmp(argv[2],"I")==0) printf("      Identical faces    : %d Visible\n",(int)df.list_face_identique.size()); else printf("      Identical faces    : %d\n",(int)df.list_face_identique.size());
    if (strcmp(argv[2],"L")==0) printf("      Localized faces    : %d Visible\n",(int)df.list_face_localise.size()); else printf("      Localized faces    : %d\n",(int)df.list_face_localise.size());
    if (strcmp(argv[2],"L")==0) printf("      Localized edges    : %d Visible\n",(int)df.list_arete_localise.size()); else  printf("      Localized faces    : %d\n",(int)df.list_face_localise.size()); 
    if (strcmp(argv[2],"L")==0) printf("      Localized vertices : %d Visible\n",(int)df.list_sommet_localise.size()); else  printf("      Localized faces    : %d\n",(int)df.list_face_localise.size());
    if (strcmp(argv[2],"M")==0) printf("      Modified faces     : %d Visible\n",(int)df.list_face_modifie.size()); else printf("      Modified faces     : %d \n",(int)df.list_face_modifie.size());
    if (strcmp(argv[2],"M")==0) printf("      Modified edges     : %d Visible\n",(int)df.list_arete_modifie.size()); else printf("      Modified edges     : %d\n",(int)df.list_arete_modifie.size());
    
    #ifdef _COMPARESTEP_
    MG_GEOMETRIE* geo=gest.get_mg_geometrie(0);
    #else
    MG_GEOMETRIE* geo=gest->get_mg_geometrie(0);
    #endif
    data_container data;
    
    
    if (strcmp(argv[2],"M")==0) affiche_structure(base,geo,A,tabcoul,data,df.list_face_modifie,df.list_arete_modifie,list_vide);
    if (strcmp(argv[2],"I")==0) affiche_structure(base,geo,A,tabcoul,data,df.list_face_identique,list_vide,list_vide);
    if (strcmp(argv[2],"S")==0) affiche_structure(base,geo,A,tabcoul,data,df.list_face_similaire,list_vide,list_vide);
    if (strcmp(argv[2],"L")==0) affiche_structure(base,geo,A,tabcoul,data,df.list_face_localise,df.list_arete_localise,df.list_sommet_localise);

    printf("\n\n");
    printf("Command on keyboard : \n");
    printf("    e   : End of vizualisation\n");
    printf("    r   : Reset Camera\n");
    printf("    1   : Vertex Numbers visible on/off\n");
    printf("    4   : Edge Numbers visible on/off\n");
    printf("    7   : Face Numbers visible on/off\n");
    printf("    v   : Vertices visible on/off\n");
    printf("    l   : Edges visible on/off\n");
    printf("    t   : Faces visible on/off\n");
    printf("    j   : Save camera on current directory\n");
    printf("    k   : Load camera from current directory\n");
    printf("    w   : Wireframe\n");
    printf("    s   : Modelframe\n");
    printf("    -   : Decrease text size\n");
    printf("    +   : Increase text size\n");
    char titre [255];
    strcpy(titre,"Visucompare File ");
    strcat(titre,argv[1]);
    strcat(titre," ");
    strcat(titre,argv[2]);
    strcat(titre," ");
    strcat(titre,argv[3]);
    
      
    #ifndef _COMPARESTEP_
    delete gest;  
    #endif
    vtkdisplay d(base,titre);
    
    
    d.init_data(data);
    d.display();
    
    return 0;
}


#pragma package(smart_init)
Revision:	1157
Committed:	Thu Jun 13 22:18:27 2024 UTC (12 months ago) by francois
File size:	22942 byte(s)
Log Message:	compatibilité Ubuntu 22.04 Suppression des refeences à Windows Ajout d'une banière