aboutsummaryrefslogtreecommitdiff
path: root/src/compare_hkl.c
diff options
context:
space:
mode:
authorThomas White <taw@bitwiz.org.uk>2010-05-04 07:38:31 -0700
committerThomas White <taw@bitwiz.org.uk>2010-05-04 08:56:17 -0700
commite447c2efb128823c93358c51a7d8f23636740f68 (patch)
tree68fb800a56d2e06cc6f3610b58a26ecc9bb6cb57 /src/compare_hkl.c
parent07d406ce05760f633f72c2c3683dd3466423321a (diff)
process_hkl: Implement --scale option
Diffstat (limited to 'src/compare_hkl.c')
-rw-r--r--src/compare_hkl.c7
1 files changed, 2 insertions, 5 deletions
diff --git a/src/compare_hkl.c b/src/compare_hkl.c
index bb07d17e..e7752078 100644
--- a/src/compare_hkl.c
+++ b/src/compare_hkl.c
@@ -53,7 +53,6 @@ int main(int argc, char *argv[])
double scale, R;
unsigned int *c1;
unsigned int *c2;
- unsigned int *cjoint;
int i;
int nc1, nc2, ncom;
@@ -134,18 +133,16 @@ int main(int argc, char *argv[])
}
}
- cjoint = new_list_count();
nc1 = 0;
nc2 = 0;
ncom = 0;
for ( i=0; i<IDIM*IDIM*IDIM; i++ ) {
- cjoint[i] = c1[i] && c2[i];
nc1 += c1[i];
nc2 += c2[i];
- ncom += cjoint[i];
+ ncom += c1[i] && c2[i];
}
STATUS("%i,%i reflections: %i in common\n", nc1, nc2, ncom);
- R = stat_r2(ref1, ref2, cjoint, IDIM*IDIM*IDIM, &scale);
+ R = stat_r2(ref1, c1, ref2, c2, &scale);
STATUS("R2 = %5.4f %% (scale=%5.2f)\n", R*100.0, scale);
if ( outfile != NULL ) {