From 8bcc2f1fd10f5c1efdf89c00d04af1f0de948741 Mon Sep 17 00:00:00 2001 From: Thomas White Date: Thu, 26 Jun 2014 16:34:59 +0200 Subject: Add Thin Ewald Sphere model --- tests/pr_p_gradient_check.c | 9 ++++++--- 1 file changed, 6 insertions(+), 3 deletions(-) (limited to 'tests/pr_p_gradient_check.c') diff --git a/tests/pr_p_gradient_check.c b/tests/pr_p_gradient_check.c index 8f7a66e4..14a2b579 100644 --- a/tests/pr_p_gradient_check.c +++ b/tests/pr_p_gradient_check.c @@ -174,6 +174,7 @@ static Crystal *new_shifted_crystal(Crystal *cr, int refine, double incr_val) return cr_new; } + static void calc_either_side(Crystal *cr, double incr_val, int *valid, long double *vals[3], int refine, PartialityModel pmodel) @@ -225,7 +226,6 @@ static void calc_either_side(Crystal *cr, double incr_val, } - static double test_gradients(Crystal *cr, double incr_val, int refine, const char *str, const char *file, PartialityModel pmodel, int quiet, int plot) @@ -450,7 +450,7 @@ int main(int argc, char *argv[]) rng = gsl_rng_alloc(gsl_rng_mt19937); - for ( i=0; i<2; i++ ) { + for ( i=0; i<3; i++ ) { UnitCell *rot; double val; @@ -459,9 +459,12 @@ int main(int argc, char *argv[]) if ( i == 0 ) { pmodel = PMODEL_SPHERE; STATUS("Testing flat sphere model:\n"); - } else { + } else if ( i == 1 ) { pmodel = PMODEL_GAUSSIAN; STATUS("Testing Gaussian model:\n"); + } else { + pmodel = PMODEL_THIN; + STATUS("Testing thin Ewald sphere model:\n"); } orientation = random_quaternion(rng); -- cgit v1.2.3