aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorThomas White <taw@physics.org>2014-05-23 16:22:49 +0200
committerThomas White <taw@physics.org>2014-05-23 16:22:49 +0200
commit40eece704c37040e0333f2de22d35feeb09a30fb (patch)
treecac03d0bccbb58ffbda8ec30fa0aefd78dc67684
parentd316a9f1edfc96a9a8d23043a3b22692080c0cc3 (diff)
Remove remnants of old convergence criterion
-rw-r--r--src/hrs-scaling.c14
1 files changed, 2 insertions, 12 deletions
diff --git a/src/hrs-scaling.c b/src/hrs-scaling.c
index b7aee812..0b95cea2 100644
--- a/src/hrs-scaling.c
+++ b/src/hrs-scaling.c
@@ -62,7 +62,6 @@ struct scale_queue_args
RefList *reference;
Crystal **crystals;
int n_started;
- double max_shift;
PartialityModel pmodel;
};
@@ -70,7 +69,6 @@ struct scale_queue_args
struct scale_worker_args
{
Crystal *crystal;
- double shift;
RefList *reference;
PartialityModel pmodel;
};
@@ -101,7 +99,6 @@ static void run_scale_job(void *vwargs, int cookie)
double num = 0.0;
double den = 0.0;
double g;
- const double G = crystal_get_osf(cr);
for ( refl = first_refl(crystal_get_reflections(cr), &iter);
refl != NULL;
@@ -132,25 +129,21 @@ static void run_scale_job(void *vwargs, int cookie)
g = num / den;
if ( !isnan(g) && !isinf(g) ) {
crystal_set_osf(cr, g);
- wargs->shift = fabs((g/G)-1.0);
} else {
- wargs->shift = 0.0;
}
}
static void finalise_scale_job(void *vqargs, void *vwargs)
{
- struct scale_queue_args *qargs = vqargs;
struct scale_worker_args *wargs = vwargs;
- if ( wargs->shift > qargs->max_shift ) qargs->max_shift = wargs->shift;
free(wargs);
}
-static double iterate_scale(Crystal **crystals, int n, RefList *reference,
- int n_threads, PartialityModel pmodel)
+static void iterate_scale(Crystal **crystals, int n, RefList *reference,
+ int n_threads, PartialityModel pmodel)
{
struct scale_queue_args qargs;
@@ -159,13 +152,10 @@ static double iterate_scale(Crystal **crystals, int n, RefList *reference,
qargs.reference = reference;
qargs.n_started = 0;
qargs.crystals = crystals;
- qargs.max_shift = 0.0;
qargs.pmodel = pmodel;
run_threads(n_threads, run_scale_job, create_scale_job,
finalise_scale_job, &qargs, n, 0, 0, 0);
-
- return qargs.max_shift;
}