--- optimisation/src/optimisation_topo.cpp	2009/07/21 15:00:12	199
+++ optimisation/src/optimisation_topo.cpp	2010/03/03 14:42:23	235
@@ -23,7 +23,6 @@
 #include "write_num.h"
 #include "ecriture_materiau.h"
 #include "affecter_materiau.h"
-#include "structure.h"
 
 int main(int argc, char **argv )
 {
@@ -40,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]    ;
@@ -69,32 +69,25 @@ int main(int argc, char **argv )
   ulong j                  ;
   ulong l                  ;
   ulong iter = 1           ;
-  ulong iter_max = 50      ;         //On sp�cifit le nombre d'it�ration maximum permit (par défaut 50)
   ulong _type              ;         //Entier correspondant au type d'�l�ment
   ulong fic_size1;
   ulong fic_size2;
   ulong fic_size3;
   ulong fic_size4;
 
-  vifl l1; vifl l2           ;          //Variables pour la dichotomie
-  vifl m= 0.18; vifl lmid    ;
+  vifl l1; vifl l2; vifl lmid;  //Variables pour la dichotomie
   vifl critere_densite = 0.0 ;
-  vifl E = 210e9             ;   	//module de young (par d�faut acier)
-  vifl nu = 0.3              ;  	//coefficient de poisson (par d�faut acier)
-  vifl rho = 7850.           ; 		//masse volumique (par d�faut acier)
-  vifl Vx = 0.               ;  	//volume de mati�re � optimiser (design)
-  vifl V1 = 0.               ;  	//volume de mati�re fixe (non-design)
-  vifl frac_vol = 0.6        ;  	//fraction volumique de la zone active (valeur par défaut)
-  vifl penal = 3.0           ;          //coefficient de p�nalisation (par défaut 3)
-  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 densite_min = 1e-3    ;          //Densit� minimale de 0.001
+  vifl E = 210e9             ;  //module de young (par defaut acier)
+  vifl nu = 0.3              ;  //coefficient de poisson (par defaut acier)
+  vifl rho = 7850.           ; 	//masse volumique (par defaut acier)
+  vifl Vx = 0.               ;  //volume de matiere fixe optimiser (design)
+  vifl V1 = 0.               ;  //volume de matiere fixe (non-design)
   vifl densite_max = 1.0     ;
+  vifl change = 1.           ;
   vifl m_item                ;
   vifl _volume               ;
   vifl m_distance            ;
-  vifl **m_V_distance        ;		//Tableau de pointeur vers les distances aux voisins de chaque éléments
+  vifl **m_V_distance        ;	//Tableau de pointeur vers les distances aux voisins de chaque elements
   vifl x1                    ;
   vifl x2                    ;
   vifl x3                    ;
@@ -106,7 +99,18 @@ int main(int argc, char **argv )
 
   ty_RSLT m_design           ;
 
-  //Acquisition des param�tres de l'�tude
+  ulong iter_max = 50           ;       //Nombre d'iteration maximum permit 
+  vifl critere_convergence = 0.1;       //Critère de convergence par défaut
+  vifl coeff_voisins = 2.5      ;       //Valeur par défaut du coefficient multipliant pour le calcul des voisins
+  vifl m = 0.2                  ;       //Limite supérieure m pour la dichotomie
+  vifl frac_vol = 0.6           ;  	//fraction volumique de la zone active f 
+  vifl penal = 3.0              ;       //coefficient de penalisation p 
+  vifl k = 3.0                  ;       //Coefficient k pour le calcul des poids Hv
+  vifl niveau = 150.            ;       //Nombre de niveaux de densites 
+  vifl seuil = 0.8              ;	//Seuil pour la visualisation des éléments conservés
+  vifl densite_min = 1e-3       ;       //Densite minimale (par défaut 0.001) 
+
+  //Acquisition des parametres de l'etude
   for (int i=0;i<argc;i++)
 	  {
 	  if (strcmp(argv[i],"-fichier")==0) strcpy(m_nom_fic_parametre,argv[i+1]);
@@ -114,6 +118,12 @@ 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]);
+	  if (strcmp(argv[i],"-k")==0) k=atof(argv[i+1]);
+	  if (strcmp(argv[i],"-rhomin")==0) densite_min=atof(argv[i+1]);
+	  if (strcmp(argv[i],"-coeffvoisins")==0) coeff_voisins=atof(argv[i+1]);
+	  if (strcmp(argv[i],"-critere")==0) critere_convergence=atof(argv[i+1]);
+	  if (strcmp(argv[i],"-m")==0) m=atof(argv[i+1]);
 	  }
   FILE* parametres=fopen(m_nom_fic_parametre,"rt");
   char message[100];
@@ -151,41 +161,41 @@ int main(int argc, char **argv )
 
   printf( "\nfait.\n" );
 
-  //Possibilité d'utiliser une seule chaîne Str (à modifier ultérieurement)
-  char str1[80];
-  char str2[80];
-  char str3[80];
-  char str4[80];
-  char str5[80];
-  char str6[80];
-
-  strcpy(str1,"./");
-  strcat(str1,m_nom_fic_etude);
-  strcat(str1,".export");
-  strcpy( &m_nom_fic_export[0], str1); 
-  strcpy(str2,"./");
-  strcat(str2,m_nom_fic_etude);
-  strcat(str2,".mail");
-  strcpy( &m_nom_fic_mail_init[0], str2);
-  strcpy(str3,"./");
-  strcat(str3,m_nom_fic_etude);
-  strcat(str3,".mail");
-  strcpy( &m_nom_fic_mail[0], str3);
-  strcpy(str4,"./");
-  strcat(str4,m_nom_fic_etude);
-  strcat(str4,"1.mail");
-  strcpy( &m_nom_fic_mail_iter[0], str4);
-  strcpy(str5,"./");
-  strcat(str5,m_nom_fic_etude);
-  strcat(str5,".comm");
-  strcpy( &m_nom_fic_comm[0], str5);
-  strcpy(str6,"./");
-  strcat(str6,m_nom_fic_etude);
-  strcat(str6,".resu");
-  strcpy( &m_nom_fic_resu[0], str6);
-  strcpy( &m_nom_fic_groupm[0]     , "./groupes.mail"      );
-  strcpy( &m_nom_fic_compliance[0] , "./compliance.txt "   );
+  char str[80];
 
+  strcpy(str,"./");
+  strcat(str,m_nom_fic_etude);
+  strcat(str,".export");
+  strcpy( &m_nom_fic_export[0], str); 
+  strcpy(str,"./");
+  strcat(str,m_nom_fic_etude);
+  strcat(str,".mail");
+  strcpy( &m_nom_fic_mail_init[0], str);
+  strcpy(str,"./");
+  strcat(str,m_nom_fic_etude);
+  strcat(str,".mail");
+  strcpy( &m_nom_fic_mail[0], str);
+  strcpy(str,"./");
+  strcat(str,m_nom_fic_etude);
+  strcat(str,"1.mail");
+  strcpy( &m_nom_fic_mail_iter[0], str);
+  strcpy(str,"./");
+  strcat(str,m_nom_fic_etude);
+  strcat(str,".comm");
+  strcpy( &m_nom_fic_comm[0], str);
+  strcpy(str,"./");
+  strcat(str,m_nom_fic_etude);
+  strcat(str,".resu");
+  strcpy( &m_nom_fic_resu[0], str);
+  strcpy(str,"./");
+  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"      );
 
   //----------------------------------------------------------------------------
   //----------------------------------------------------------------------------
@@ -290,7 +300,7 @@ int main(int argc, char **argv )
   //Allocation de la m�moire n�cessaire dans les tableaux
   ulong *compteur = (ulong*) malloc (sizeof(ulong)*(_dernier_element+1))  ;
   vifl *densite = (vifl*) malloc (sizeof(vifl)*(_dernier_element+1))      ;     //vecteur des variables densit�s pour chaque �l�ment
-  ulong *densite_active = (ulong*) malloc (sizeof(ulong)*(_dernier_element+1));  //vecteur qui d�termine si la densit� est active
+  ulong *densite_active = (ulong*) malloc (sizeof(ulong)*(_dernier_element+1)); //vecteur qui d�termine si la densit� est active
   vifl *densite_new = (vifl*) malloc (sizeof(vifl)*(_dernier_element+1))  ;     //vecteur des variables densit�s pour la nouvelle it�ration
   vifl *volume = (vifl*) malloc (sizeof(vifl)*(_dernier_element+1))       ;     //vecteur des volumes pour chaque �l�ment
   vifl *Denergie = (vifl*) malloc (sizeof(vifl)*(_dernier_element+1))     ;
@@ -311,7 +321,7 @@ int main(int argc, char **argv )
 	compteur[i] = 1;
     	maille_niveau[i]=0.;
   	}
-
+  vifl volume_total = 0.;
   //Calcul du volume de design et du volume total du maillage 3D
   for (i=_premier_element;i<=_dernier_element;i++)
   	{
@@ -329,9 +339,9 @@ int main(int argc, char **argv )
 		densite[i] = 1.;
 		V1 = V1 + _volume;
 		}		
-  	
+  	volume_total = volume_total + _volume;
   	}
-
+  printf( "Volume total : %lf\n",volume_total);
   ///////////////////////////////////////
   ///////Calcul du voisinage et//////////
   //des pond�rations de r�gulation///
@@ -353,7 +363,7 @@ int main(int argc, char **argv )
 		                      	  (ptr_o_element_com *) &m_element    );
 	    	m_element->ren_volume   (                        &_volume     );
 	    	volume[i] = _volume                                            ;
-	    	distance_ref[i] = 2.5*pow((12.*_volume/co_rc_2),co_1_3)        ;
+	    	distance_ref[i] = coeff_voisins*pow((12.*_volume/co_rc_2),co_1_3);
 		m_ptr_etude->ren_nbr_voisins_sphere( (ulong)      i            ,
                                        		    (vifl *) &distance_ref[i]  ,
 						   (ulong *) &m_nbr_voisins[i]);
@@ -387,7 +397,7 @@ int main(int argc, char **argv )
     	Ctot[i] = 0.0;
   	}
 
-  while ((fabs(change) > 0.1) && (iter <= iter_max))
+  while ((fabs(change) > critere_convergence) && (iter <= iter_max))
   {
 
         ///////////////////////////////////////
@@ -649,7 +659,7 @@ int main(int argc, char **argv )
 
         //D�termination des multiplicateurs de Lagrange par dichotomie
         l1= 0. ;
-        l2= 1e8 ;
+        l2= 1e8;
         while ((l2-l1) > 1e-12)
         {
         lmid = 0.5*(l2+l1);
@@ -746,27 +756,37 @@ 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" );
 
-  //Libération de la mémoire utilisée
-//   free( compteur        );
-//   free( densite         );
-//   free( densite_active  );
-//   free( densite_new     );
-//   free( volume          );
-//   free( Denergie        );
-//   free( distance_ref    );
-//   free( Ctot            );
-//   free( nmailleniv      );
-//   free( maille_niveau   );
-//   free( m_V_distance    );
-//   free( m_nbr_voisins   );
-//   free( m_V_nos_voisins );
+  //Libération de la mémoire utilisée dans les tableaux
+   free( compteur        );
+   free( densite         );
+   free( densite_active  );
+   free( densite_new     );
+   free( volume          );
+   free( Denergie        );
+   free( distance_ref    );
+   free( Ctot            );
+   free( nmailleniv      );
+   free( maille_niveau   );
+   free( m_V_distance    );
+   free( m_nbr_voisins   );
+   free( m_V_nos_voisins );
 
   finit_etude_com( ) ;