diff options
| -rw-r--r-- | src/post-refinement.c | 10 |
1 files changed, 0 insertions, 10 deletions
diff --git a/src/post-refinement.c b/src/post-refinement.c index ab2d4425..b2e4e2e8 100644 --- a/src/post-refinement.c +++ b/src/post-refinement.c @@ -56,7 +56,6 @@ static double partiality_gradient(double r, double profile_radius) /* Calculate degree of penetration */ q = (r + profile_radius)/(2.0*profile_radius); - //STATUS("q=%f\n", q); /* dp/dq */ dpdq = 6.0*(q-pow(q, 2.0)); @@ -108,15 +107,11 @@ static double gradient(struct image *image, int k, Reflection *refl, double r) xl = hs*asx + ks*bsx + ls*csx; yl = hs*asy + ks*bsy + ls*csy; zl = hs*asz + ks*bsz + ls*csz; - //STATUS("%3i %3i %3i\n", hs, ks, ls); ds = 2.0 * resolution(image->indexed_cell, hs, ks, ls); tt = angle_between(0.0, 0.0, 1.0, xl, yl, zl+1.0/image->lambda); azix = angle_between(1.0, 0.0, 0.0, xl, yl, 0.0); aziy = angle_between(0.0, 1.0, 0.0, xl, yl, 0.0); - //STATUS("d*=%.2e, 2theta=%.2f deg, azi=%.2f deg\n", - // ds, rad2deg(tt), rad2deg(azi)); - get_partial(refl, &r1, &r2, &p, &clamp_low, &clamp_high); /* Calculate the gradient of partiality wrt excitation error. */ @@ -127,9 +122,6 @@ static double gradient(struct image *image, int k, Reflection *refl, double r) if ( clamp_high == 0 ) { g += partiality_gradient(r2, r); } - //STATUS("clamp status low=%i high=%i\n", clamp_low, clamp_high); - //STATUS("excitation errors %e %e\n", r1, r2); - //STATUS("dp/dr = %e\n", g); /* For many gradients, just multiply the above number by the gradient * of excitation error wrt whatever. */ @@ -187,8 +179,6 @@ static void apply_cell_shift(UnitCell *cell, int k, double shift) &bsx, &bsy, &bsz, &csx, &csy, &csz); - //STATUS("Shifting by %e\n", shift); - switch ( k ) { case REF_ASX : asx += shift; break; |