adjusted eps in setBayesModel and keep const in dprofile

include/fonction.h

@@ -109,7 +109,7 @@ double  dist(struct point A, struct point B);
 double  dlumcosmo1(double z);
 void      do_itos(double **im, struct pixlist *pl, int npl, double dlsds, double zs, double **source, double **erreur, int **imult);
 double  d_poisson(double xm, int *idum);
-double  d_profil(double x, double y, struct galaxie *gal);
+double  d_profil(double x, double y, const struct galaxie *gal);
 double  d_random(int *idum);
 double  dratio(double zl, double zs);
 void dratio_gal(struct galaxie *arclet, double zl);

src/d_profil.c

@@ -23,7 +23,7 @@
  *      alpha=1/(e*l)
  */
 
-double d_profil(double x, double y, struct galaxie *gal)
+double d_profil(double x, double y, const struct galaxie *gal)
 {
     double xx, yy, xxx, yyy;
     double res;
@@ -33,8 +33,8 @@ double d_profil(double x, double y, struct galaxie *gal)
     res = 0;
 
     //jrichard
-    //Recompute E.b according to E.a and eps
-    gal->E.b=gal->E.a*(1-gal->eps)/(1+gal->eps);
+    //Test E.b according to E.a and eps
+    //fprintf(stderr,"%f\n",gal->E.a*(1-gal->eps)/(1+gal->eps)-gal->E.b);
 
     if (gal->c == 's')
     {

src/readBayesModels.c

@@ -294,6 +294,9 @@ void    setBayesModel( long int iVal, long int nVal, double **array)
                 if( sblock[i][ipx] != 0 )
                     o_set_source( &source[i], ipx, method(nVal, iParam++, array));
 
+            if((sblock[i][SEPS] != 0 )|| (sblock[i][SA] !=0))
+                source[i].E.b = source[i].E.a * (1. - source[i].eps) / (1. + source[i].eps);
+
        }
     }
 
@@ -845,7 +848,6 @@ static void convertArray(long int nVal, double **array )
                     {
                         if( ipx == STHETA )
                             array[nParam][iVal] = array[nParam][iVal] * DTR;
-    
                         nParam ++;
                     }
     }