aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--src/post-refinement.c10
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;