diff options
| author | Thomas White <taw@physics.org> | 2010-10-29 18:19:52 +0200 |
|---|---|---|
| committer | Thomas White <taw@physics.org> | 2012-02-22 15:27:04 +0100 |
| commit | d16690a8e760d3cbe9a1285b6f6a8f4a0e8e727d (patch) | |
| tree | 2efdbb0850116b97f8e4e679e0436f02b41b0210 /src/facetron.c | |
| parent | cdf5df91471b4fd14c85fb06f3dab94f734750fb (diff) | |
facetron: Output counts
Diffstat (limited to 'src/facetron.c')
| -rw-r--r-- | src/facetron.c | 14 |
1 files changed, 8 insertions, 6 deletions
diff --git a/src/facetron.c b/src/facetron.c index 02d2603f..9bd518cf 100644 --- a/src/facetron.c +++ b/src/facetron.c @@ -192,14 +192,13 @@ static void refine_all(struct image *images, int n_total_patterns, static void estimate_full(struct image *images, int n_total_patterns, struct detector *det, const char *sym, - ReflItemList *obs, double *i_full, int nthreads) + ReflItemList *obs, double *i_full, unsigned int *cts, + int nthreads) { int i; - unsigned int *cts; struct integrate_args *tasks; pthread_mutex_t list_lock = PTHREAD_MUTEX_INITIALIZER; - cts = new_list_count(); clear_items(obs); tasks = malloc(n_total_patterns * sizeof(struct integrate_args)); @@ -250,6 +249,7 @@ int main(int argc, char *argv[]) int config_checkprefix = 1; struct detector *det; double *i_full; + unsigned int *cts; ReflItemList *obs; int i; int n_total_patterns; @@ -433,8 +433,10 @@ int main(int argc, char *argv[]) fclose(fh); free(prefix); + cts = new_list_count(); + /* Make initial estimates */ - estimate_full(images, n_total_patterns, det, sym, obs, i_full, + estimate_full(images, n_total_patterns, det, sym, obs, i_full, cts, nthreads); /* Iterate */ @@ -448,12 +450,12 @@ int main(int argc, char *argv[]) /* Re-estimate all the full intensities */ estimate_full(images, n_total_patterns, det, sym, obs, i_full, - nthreads); + cts, nthreads); } /* Output results */ - write_reflections(outfile, obs, i_full, NULL, NULL, NULL, 1.0); + write_reflections(outfile, obs, i_full, NULL, cts, NULL, 1.0); STATUS("Sigma(I) values in output file are not (yet) meaningful.\n"); /* Clean up */ |