--- optimisation/src/optimisation_topo.cpp	2009/07/27 13:53:35	209
+++ optimisation/src/optimisation_topo.cpp	2009/11/10 19:01:30	218
@@ -39,6 +39,7 @@ int main(int argc, char **argv )
 
   char              m_nom_fic_parametre[co_lcc_max_nom_fic] ;
   char              m_nom_fic_etude[co_lcc_max_nom_fic]     ;
+  char		    m_nom_fic_densite[co_lcc_max_nom_fic]   ;
   char              m_nom_gm[co_lcc_max_nom_gm]             ;
   char              m_nom_champ[co_lcc_max_nom_champ_resu]  ;
   char              m_nom_fic_export[co_lcc_max_nom_fic]    ;
@@ -88,6 +89,7 @@ int main(int argc, char **argv )
   vifl k = 3.0               ;          //Coefficient k pour le calcul des poids Hv
   vifl niveau = 150.         ;          //Nombre de niveaux de densit�s (par défaut 150)
   vifl change = 1.           ;
+  vifl seuil = 0.8           ;		//Seuil pour la visualisation des éléments conservés (0.8 par défaut)
   vifl densite_min = 1e-3    ;          //Densit� minimale de 0.001
   vifl densite_max = 1.0     ;
   vifl m_item                ;
@@ -105,6 +107,8 @@ int main(int argc, char **argv )
 
   ty_RSLT m_design           ;
 
+
+
   //Acquisition des param�tres de l'�tude
   for (int i=0;i<argc;i++)
 	  {
@@ -113,6 +117,7 @@ int main(int argc, char **argv )
 	  if (strcmp(argv[i],"-penal")==0) penal=atof(argv[i+1]);
 	  if (strcmp(argv[i],"-niveau")==0) niveau=atof(argv[i+1]);
 	  if (strcmp(argv[i],"-nbiter")==0) iter_max=atoi(argv[i+1]);
+	  if (strcmp(argv[i],"-seuil")==0) seuil=atoi(argv[i+1]);
 	  }
   FILE* parametres=fopen(m_nom_fic_parametre,"rt");
   char message[100];
@@ -180,6 +185,10 @@ int main(int argc, char **argv )
   strcat(str,m_nom_fic_etude);
   strcat(str,"_evol.txt");
   strcpy( &m_nom_fic_compliance[0], str);
+  strcpy(str,"./");
+  strcat(str,m_nom_fic_etude);
+  strcat(str,"_densite.txt");
+  strcpy( &m_nom_fic_densite[0], str);
   strcpy( &m_nom_fic_groupm[0]     , "./groupes.mail"      );
 
   //----------------------------------------------------------------------------
@@ -741,9 +750,19 @@ int main(int argc, char **argv )
   }
   fclose(compliance_iter);
 
+  //On imprime aussi le vecteur densite final dans un fichier texte pour usage ultérieur
+  FILE *densite_final = fopen (m_nom_fic_densite, "wt" );
+  fprintf(densite_final,"Vecteur Densite\n");
+  for (i=_premier_element;i<=_dernier_element;i++)
+  {
+      fprintf(densite_final,"%lu %.12lf\n",i,densite[i]);
+  }
+  fclose(densite_final);
+
   //Impression des temps de calcul
   printf("Temps de calcul pour le voisinage : %lu secondes\n",temps_calcul_voisins);
   affiche_chrono();
+  printf("\n");
 
   printf( "\n\n-------------------------------------------------------\n" );
   printf(     "-- finitialisation des modules ------------------------\n" );