diff options
| author | Thomas White <taw@bitwiz.org.uk> | 2010-10-14 01:14:29 +0200 |
|---|---|---|
| committer | Thomas White <taw@physics.org> | 2012-02-22 15:27:03 +0100 |
| commit | 67d8cf120ba7e8140449379d0153cef93795c813 (patch) | |
| tree | 24bba9a620a3e91b1ea882e11be2b88a9145930a /src/compare_hkl.c | |
| parent | e9803e863335882f161295af184d91ebd5e9e02c (diff) | |
compare_hkl: Add Rinternal(I)
Diffstat (limited to 'src/compare_hkl.c')
| -rw-r--r-- | src/compare_hkl.c | 10 |
1 files changed, 8 insertions, 2 deletions
diff --git a/src/compare_hkl.c b/src/compare_hkl.c index 959f1fed..ed8636fb 100644 --- a/src/compare_hkl.c +++ b/src/compare_hkl.c @@ -293,6 +293,7 @@ int main(int argc, char *argv[]) char *bfile = NULL; char *sym = NULL; double scale, scale_r2, scale_rdig, R1, R2, R1i, Rdiff, pearson; + double scale_rintint; int i, ncom; ReflItemList *i1, *i2, *icommon; int config_shells = 0; @@ -451,11 +452,16 @@ int main(int argc, char *argv[]) STATUS("R1(I) = %5.4f %% (scale=%5.2e)\n", R1i*100.0, scale); Rdiff = stat_rdiff_ignore(ref1, ref2_transformed, icommon, &scale_rdig); - STATUS("Rdiff(F) = %5.4f %% (scale=%5.2e) (ignoring negative intensities)\n", + STATUS("Rint(F) = %5.4f %% (scale=%5.2e) (ignoring negative intensities)\n", Rdiff*100.0, scale); Rdiff = stat_rdiff_zero(ref1, ref2_transformed, icommon, &scale); - STATUS("Rdiff(F) = %5.4f %% (scale=%5.2e) (zeroing negative intensities)\n", + STATUS("Rint(F) = %5.4f %% (scale=%5.2e) (zeroing negative intensities)\n", + Rdiff*100.0, scale); + + Rdiff = stat_rdiff_intensity(ref1, ref2_transformed, icommon, + &scale_rintint); + STATUS("Rint(I) = %5.4f %% (scale=%5.2e)\n", Rdiff*100.0, scale); pearson = stat_pearson_i(ref1, ref2_transformed, icommon); |