Add tests/centering_check

3 files changed, 120 insertions, 3 deletions

diff --git a/Makefile.am b/Makefile.am
index 44e707ea..953a04df 100644
--- a/Makefile.am
+++ b/Makefile.am
@@ -7,12 +7,13 @@ bin_PROGRAMS = src/pattern_sim src/process_hkl src/get_hkl src/indexamajig \
                src/partialator src/check_hkl src/partial_sim
 
 noinst_PROGRAMS = tests/list_check tests/integration_check \
-                  tests/pr_gradient_check tests/symmetry_check
+                  tests/pr_gradient_check tests/symmetry_check \
+                  tests/centering_check
 
 TESTS = tests/list_check tests/first_merge_check tests/second_merge_check \
         tests/third_merge_check tests/fourth_merge_check \
         tests/integration_check tests/pr_gradient_check \
-        tests/symmetry_check
+        tests/symmetry_check tests/centering_check
 
 EXTRA_DIST += tests/first_merge_check tests/second_merge_check \
               tests/third_merge_check tests/fourth_merge_check
@@ -80,10 +81,12 @@ tests_integration_check_SOURCES = tests/integration_check.c
 
 tests_symmetry_check_SOURCES = tests/symmetry_check.c
 
-
 tests_pr_gradient_check_SOURCES = tests/pr_gradient_check.c \
                                   src/post-refinement.c
 
+tests_centering_check_SOURCES = tests/centering_check.c libcrystfel/src/cell.c \
+                                libcrystfel/src/cell-utils.c
+
 INCLUDES = -I$(top_srcdir)/libcrystfel/src -I$(top_srcdir)/data
 
 EXTRA_DIST += src/dw-hdfsee.h src/hdfsee.h src/render_hkl.h \
diff --git a/tests/.gitignore b/tests/.gitignore
index c2a45f5a..16acffc3 100644
--- a/tests/.gitignore
+++ b/tests/.gitignore
@@ -5,4 +5,5 @@ gpu_sim_check
 integration_check
 pr_gradient_check
 symmetry_check
+centering_check
 .dirstamp
diff --git a/tests/centering_check.c b/tests/centering_check.c
new file mode 100644
index 00000000..0dc79a5f
--- /dev/null
+++ b/tests/centering_check.c
@@ -0,0 +1,113 @@
+/*
+ * centering_check.c
+ *
+ * Check that centering of cells works
+ *
+ * Copyright © 2012 Thomas White <taw@physics.org>
+ *
+ * This file is part of CrystFEL.
+ *
+ * CrystFEL is free software: you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License as published by
+ * the Free Software Foundation, either version 3 of the License, or
+ * (at your option) any later version.
+ *
+ * CrystFEL is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with CrystFEL.  If not, see <http://www.gnu.org/licenses/>.
+ *
+ */
+
+#ifdef HAVE_CONFIG_H
+#include <config.h>
+#endif
+
+
+#include <stdlib.h>
+#include <stdio.h>
+#include <stdarg.h>
+
+#include <cell.h>
+#include <cell-utils.h>
+
+
+static int check_cell(UnitCell *cell)
+{
+	int err = 0;
+
+	if ( !cell_is_sensible(cell) ) {
+		ERROR("Warning: Unit cell parameters are not sensible.\n");
+		err = 1;
+	}
+
+	if ( !bravais_lattice(cell) ) {
+		ERROR("Warning: Unit cell is not a conventional Bravais"
+		      " lattice.\n");
+		err = 1;
+	}
+
+	if ( !right_handed(cell) ) {
+		ERROR("Warning: Unit cell is not right handed.\n");
+		err = 1;
+	}
+
+	return err;
+}
+
+
+static int check_centering()
+{
+	UnitCell *cell;
+	UnitCell *n;
+	int fail = 0;
+
+	cell = cell_new();
+
+	cell_set_parameters(cell, 10e-10, 20e-10, 30e-10,
+	                    deg2rad(90.0), deg2rad(90.0), deg2rad(100.0));
+	cell_set_centering(cell, 'C');
+	cell_set_lattice_type(cell, L_MONOCLINIC);
+	cell_set_unique_axis(cell, 'c');
+	if ( check_cell(cell) ) fail = 1;
+	n = uncenter_cell(cell);
+	if ( check_cell(n) ) fail = 1;
+
+	return fail;
+}
+
+
+
+int main(int argc, char *argv[])
+{
+	int fail = 0;
+
+	/* Triclinic P */
+	check_centering(10e-10, 20e-10, 30e-10, 100.0, 120.0, 140.0,
+	                L_TRICLINIC, 'P', '*',
+	                );
+	/* Monoclinic P */
+	/* Monoclinic A */
+	/* Monoclinic B */
+	/* Monoclinic C */
+	/* Orthorhombic P */
+	/* Orthorhombic A */
+	/* Orthorhombic B */
+	/* Orthorhombic C */
+	/* Orthorhombic I */
+	/* Orthorhombic F */
+	/* Tetragonal P */
+	/* Tetragonal I */
+	/* Rhombohedral R */
+	/* Hexagonal P */
+	/* Hexagonal H (PDB-speak for rhombohedral) */
+	/* Cubic P */
+	/* Cubic I */
+	/* Cubic F */
+
+
+	return fail;
+}