improved

Jean-Philippe Bernard
1 parent 1ad87749
Showing 1 changed file with 8 additions and 3 deletions Show diff stats
src/idl/dustem_brute_force_fit.pro
@@ -110,7 +110,7 @@ ENDIF
 ;TABLE_PMINS     DOUBLE    = Array[3]
 ;TABLE_PMAXS     DOUBLE    = Array[3]
  
-;stop
+stop
  
 ;=== These are the full arrays needed for calculations
 chi2s=dblarr(Nseds)
@@ -190,11 +190,16 @@ FOR i=0L,Nparams-1 DO BEGIN
 	weighted_params[i]=total(reform(1.d0*grid_param_values[*,i])*weights)/tot_weights
 ENDFOR
  
+;stop
 ;==== Check bumping to table edges.
 params_hit=intarr(Nparams)
 FOR i=0L,Nparams-1 DO BEGIN
-  IF params[i] EQ table_pmins[i] THEN params_hit[i]=-1
-  IF params[i] EQ table_pmaxs[i] THEN params_hit[i]=1
+	diff_min=abs(params[i]-table_pmins[i])/table_pmins[i]*100.
+	diff_max=abs(params[i]-table_pmaxs[i])/table_pmaxs[i]*100.
+	IF diff_min LE 1. THEN params_hit[i]=-1
+	IF diff_max LE 1. THEN params_hit[i]=1
+  ;IF params[i] LE table_pmins[i] THEN params_hit[i]=-1
+  ;IF params[i] GE table_pmaxs[i] THEN params_hit[i]=1
 ENDFOR
  
 ;==== Compute parameter uncertainties