diff options
author | Thomas White <taw@physics.org> | 2015-03-10 13:51:58 +0100 |
---|---|---|
committer | Thomas White <taw@physics.org> | 2015-03-12 16:37:25 +0100 |
commit | 5e34b26a7965eafbec680148367697a371de35e8 (patch) | |
tree | 985d401d9f864b58db42b7531e926f044dec0b7c /src | |
parent | e4a0136db78d7d4204fb8c2d712de65c7a176973 (diff) |
Remove debugging stuff
Diffstat (limited to 'src')
-rw-r--r-- | src/hrs-scaling.c | 36 |
1 files changed, 0 insertions, 36 deletions
diff --git a/src/hrs-scaling.c b/src/hrs-scaling.c index 0225d625..52850fac 100644 --- a/src/hrs-scaling.c +++ b/src/hrs-scaling.c @@ -109,8 +109,6 @@ static void run_scale_job(void *vwargs, int cookie) double c0, c1, cov00, cov01, cov11, chisq; double G, B; int r; - int debug = 0; - FILE *fh; /* If this crystal's scaling was dodgy, it doesn't contribute to the * merged intensities */ @@ -124,13 +122,6 @@ static void run_scale_job(void *vwargs, int cookie) return; } - if ( wargs->crystal_number == 0 ) { - char fn[256]; - snprintf(fn, 255, "scale-debug.%i", wargs->crystal_number); - debug = 1; - fh = fopen(fn, "w"); - } - for ( refl = first_refl(crystal_get_reflections(cr), &iter); refl != NULL; refl = next_refl(refl, iter) ) @@ -177,10 +168,6 @@ static void run_scale_job(void *vwargs, int cookie) x[n] = res*res; y[n] = log(Ihl/Ih); w[n] = 1.0;//log(get_esd_intensity(refl)*corr); - if ( debug ) { - fprintf(fh, "%i %10.2e %10.2f %10.2f %.2f %.2f %.2f\n", - n, x[n], y[n], w[n], Ihl, Ih, res/1e9); - } n++; } @@ -195,12 +182,6 @@ static void run_scale_job(void *vwargs, int cookie) G = 1.0/exp(c0); B = -c1/2.0; - if ( debug ) { - printf("intercept = %e, gradient = %e\n", c0, c1); - printf("Scale factor = %f\n", G); - printf("Relative B = %.2f A^2\n", B * 1e20); - fclose(fh); - } free(x); free(y); @@ -284,8 +265,6 @@ static void run_merge_job(void *vwargs, int cookie) Reflection *refl; RefListIterator *iter; double G, B; - FILE *fh; - int debug = 0; /* If this crystal's scaling was dodgy, it doesn't contribute to the * merged intensities */ @@ -294,13 +273,6 @@ static void run_merge_job(void *vwargs, int cookie) G = crystal_get_osf(cr); B = crystal_get_Bfac(cr); - if ( wargs->crystal_number == 0 ) { - char fn[256]; - snprintf(fn, 255, "merge-debug.%i", wargs->crystal_number); - debug = 1; - fh = fopen(fn, "w"); - } - for ( refl = first_refl(crystal_get_reflections(cr), &iter); refl != NULL; refl = next_refl(refl, iter) ) @@ -364,18 +336,11 @@ static void run_merge_job(void *vwargs, int cookie) den += 1.0; red++; - if ( debug ) { - fprintf(fh, "%.2f %.2f %f %e\n", res/1e9, - Ihl*G*exp(2.0*B*res), - G, B); - } - set_temp1(f, num); set_temp2(f, den); set_redundancy(f, red); unlock_reflection(f); } - if ( debug ) fclose(fh); } @@ -658,7 +623,6 @@ RefList *scale_intensities(Crystal **crystals, int n, } } norm_sf = total_sf / n_sf; - STATUS("norm = %f\n", norm_sf); for ( j=0; j<n; j++ ) { crystal_set_osf(crystals[j], crystal_get_osf(crystals[j])/norm_sf); |