--- optimisation/src/optimisation_topo.cpp	2009/11/10 19:01:30	218
+++ optimisation/src/optimisation_topo.cpp	2010/03/03 14:42:23	235
@@ -69,33 +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 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 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                    ;
@@ -107,9 +99,18 @@ int main(int argc, char **argv )
 
   ty_RSLT m_design           ;
 
+  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 param�tres de l'�tude
+  //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]);
@@ -118,6 +119,11 @@ int main(int argc, char **argv )
 	  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];
@@ -294,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))     ;
@@ -315,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++)
   	{
@@ -333,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///
@@ -357,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]);
@@ -391,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))
   {
 
         ///////////////////////////////////////
@@ -653,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);