Image, feature and unit cell infrastructure

Brought across from DTR and Synth3D

1 files changed, 82 insertions, 0 deletions

diff --git a/src/sim-main.c b/src/sim-main.c
new file mode 100644
index 00000000..38ad45bc
--- /dev/null
+++ b/src/sim-main.c
@@ -0,0 +1,82 @@
+/*
+ * sim-main.c
+ *
+ * Simulate test data
+ *
+ * (c) 2006-2009 Thomas White <thomas.white@desy.de>
+ *
+ * template_index - Indexing diffraction patterns by template matching
+ *
+ */
+
+
+#ifdef HAVE_CONFIG_H
+#include <config.h>
+#endif
+
+#include <stdarg.h>
+#include <stdlib.h>
+#include <stdio.h>
+#include <string.h>
+#include <unistd.h>
+
+#include "image.h"
+#include "relrod.h"
+#include "cell.h"
+#include "utils.h"
+
+
+static void main_show_help(const char *s)
+{
+	printf("Syntax: %s <file1.h5> <file2.h5> [...]\n\n", s);
+	printf("Index diffraction patterns\n\n");
+	printf("  -h              Display this help message\n");
+}
+
+
+int main(int argc, char *argv[])
+{
+	int c;
+	ImageList *list;
+	UnitCell *cell;
+	struct image image;
+
+	while ((c = getopt(argc, argv, "h")) != -1) {
+
+		switch ( c ) {
+
+			case 'h' : {
+				main_show_help(argv[0]);
+				return 0;
+			}
+
+		}
+
+	}
+
+	printf("Generating test data...\n");
+	list = image_list_new();
+	image.width = 512;
+	image.height = 512;
+	image.tilt = 0.0;
+	image.omega = 0.0;
+	image.fmode = FORMULATION_CLEN;
+	image.x_centre = 128;
+	image.y_centre = 128;
+	image.camera_len = 1.0; /* one metre */
+	image.resolution = 5120;
+	image.lambda = 0.6e-9;
+	image.data = malloc(512*512*2);
+	image_add(list, &image);
+
+	cell = cell_new_from_parameters(1.0,
+	                                1.0,
+	                                1.0,
+	                                deg2rad(90.0),
+	                                deg2rad(90.0),
+	                                deg2rad(90.0));
+
+	get_reflections(&image, cell);
+
+	return 0;
+}