aboutsummaryrefslogtreecommitdiff
path: root/tests/pr_p_gradient_check.c
diff options
context:
space:
mode:
authorThomas White <taw@physics.org>2014-08-14 15:07:21 +0200
committerThomas White <taw@physics.org>2014-09-25 10:53:56 +0200
commit629934d82e202ea04b334c49efffe09aaa0f1c4e (patch)
treee12b6ca6bd72f76e8e687317eec8fc659d6c712c /tests/pr_p_gradient_check.c
parenta06a3f67f57de0bc85982976b9ea6d598598e014 (diff)
Remove "sphere", "thin" and "gaussian" partiality models, add "scgaussian"
Diffstat (limited to 'tests/pr_p_gradient_check.c')
-rw-r--r--tests/pr_p_gradient_check.c31
1 files changed, 4 insertions, 27 deletions
diff --git a/tests/pr_p_gradient_check.c b/tests/pr_p_gradient_check.c
index 82ed2b02..b224f9b4 100644
--- a/tests/pr_p_gradient_check.c
+++ b/tests/pr_p_gradient_check.c
@@ -71,21 +71,6 @@ static void scan_partialities(RefList *reflections, RefList *compare,
}
get_partial(refl2, &r1, &r2, &p, &clamp_low, &clamp_high);
- if ( clamp_low && clamp_high && (pmodel != PMODEL_SCSPHERE) ) {
- if ( !within_tolerance(p, 1.0, 0.001) ) {
-
- signed int h, k, l;
-
- get_indices(refl, &h, &k, &l);
-
- ERROR("%3i %3i %3i - double clamped but"
- " partiality not close to 1.0 (%5.2f)\n",
- h, k, l, p);
-
- }
- valid[i] = 0;
- }
-
vals[idx][i] = p;
i++;
}
@@ -451,26 +436,18 @@ int main(int argc, char *argv[])
rng = gsl_rng_alloc(gsl_rng_mt19937);
- for ( i=0; i<4; i++ ) {
+ for ( i=0; i<2; i++ ) {
UnitCell *rot;
double val;
PartialityModel pmodel;
if ( i == 0 ) {
- pmodel = PMODEL_SPHERE;
- STATUS("Testing flat sphere model:\n");
- } else if ( i == 1 ) {
- pmodel = PMODEL_GAUSSIAN;
- /* FIXME: Gradients for Gaussian model are not good */
- STATUS("NOT testing Gaussian model.\n");
- continue;
- } else if ( i == 2 ) {
- pmodel = PMODEL_THIN;
- STATUS("Testing Thin Ewald Sphere model:\n");
- } else if ( i == 3 ) {
pmodel = PMODEL_SCSPHERE;
STATUS("Testing SCSphere model:\n");
+ } else if ( i == 1 ) {
+ pmodel = PMODEL_SCGAUSSIAN;
+ STATUS("Testing SCGaussian model.\n");
} else {
ERROR("WTF?\n");
return 1;