diff options
| -rw-r--r-- | src/partialator.c | 27 |
1 files changed, 27 insertions, 0 deletions
diff --git a/src/partialator.c b/src/partialator.c index 00d1865b..b1e32c2f 100644 --- a/src/partialator.c +++ b/src/partialator.c @@ -689,6 +689,33 @@ int main(int argc, char *argv[]) /* Output results */ write_reflist(outfile, full); + /* Output split results */ + char tmp[1024]; + RefList *split; + Crystal *crystals1[n_crystals]; + Crystal *crystals2[n_crystals]; + int n_crystals1 = 0; + int n_crystals2 = 0; + for ( i=0; i<n_crystals; i++ ) { + if ( i % 2 ) { + crystals1[n_crystals1] = crystals[i]; + n_crystals1++; + } else { + crystals2[n_crystals2] = crystals[i]; + n_crystals2++; + } + } + snprintf(tmp, 1024, "%s1", outfile); + split = lsq_intensities(crystals1, n_crystals1, nthreads, + pmodel, min_measurements); + write_reflist(tmp, split); + reflist_free(split); + snprintf(tmp, 1024, "%s2", outfile); + split = lsq_intensities(crystals2, n_crystals2, nthreads, + pmodel, min_measurements); + write_reflist(tmp, split); + reflist_free(split); + /* Dump parameters */ FILE *fh; fh = fopen("partialator.params", "w"); |