Refinement diagnostics...

git-svn-id: svn://cook.msm.cam.ac.uk:745/diff-tomo/dtr@259 bf6ca9ba-c028-0410-8290-897cf20841d1

1 files changed, 23 insertions, 1 deletions

diff --git a/src/iprtest.c b/src/iprtest.c
index 6ff7651..4a06a4a 100644
--- a/src/iprtest.c
+++ b/src/iprtest.c
@@ -23,6 +23,7 @@
 #include "reproject.h"
 #include "refine.h"
 #include "utils.h"
+#include "mapping.h"
 
 typedef struct {
 	
@@ -110,6 +111,27 @@ int main(int argc, char *argv[]) {
 	GtkWidget *window;
 	GtkWidget *drawing;
 	IPRTestContext *ctx;
+	double dlx, dly, dlz;
+	
+	mapping_rotate(1.0, 0.0, 0.0, &dlx, &dly, &dlz, 0.0, 0.0);
+	printf("omega=0.0, tilt=0.0   => image 'x' = { %f %f %f }\n", dlx, dly, dlz);
+	mapping_rotate(0.0, 1.0, 0.0, &dlx, &dly, &dlz, 0.0, 0.0);
+	printf("omega=0.0, tilt=0.0   => image 'y' = { %f %f %f }\n", dlx, dly, dlz);
+
+	mapping_rotate(1.0, 0.0, 0.0, &dlx, &dly, &dlz, deg2rad(90.0), 0.0);
+	printf("omega=90.0, tilt=0.0   => image 'x' = { %f %f %f }\n", dlx, dly, dlz);
+	mapping_rotate(0.0, 1.0, 0.0, &dlx, &dly, &dlz, deg2rad(90.0), 0.0);
+	printf("omega=90.0, tilt=0.0   => image 'y' = { %f %f %f }\n", dlx, dly, dlz);
+
+	mapping_rotate(1.0, 0.0, 0.0, &dlx, &dly, &dlz, 0.0, deg2rad(90.0));
+	printf("omega=0.0, tilt=90.0   => image 'x' = { %f %f %f }\n", dlx, dly, dlz);
+	mapping_rotate(0.0, 1.0, 0.0, &dlx, &dly, &dlz, 0.0, deg2rad(90.0));
+	printf("omega=0.0, tilt=90.0   => image 'y' = { %f %f %f }\n", dlx, dly, dlz);
+
+	mapping_rotate(1.0, 0.0, 0.0, &dlx, &dly, &dlz, deg2rad(90.0), deg2rad(90.0));
+	printf("omega=90.0, tilt=90.0   => image 'x' = { %f %f %f }\n", dlx, dly, dlz);
+	mapping_rotate(0.0, 1.0, 0.0, &dlx, &dly, &dlz, deg2rad(90.0), deg2rad(90.0));
+	printf("omega=90.0, tilt=90.0   => image 'y' = { %f %f %f }\n", dlx, dly, dlz);
 	
 	gtk_init(&argc, &argv);
 	
@@ -117,7 +139,7 @@ int main(int argc, char *argv[]) {
 	
 	ctx->image = malloc(sizeof(ImageRecord));
 	ctx->image->image = NULL;
-	ctx->image->tilt = deg2rad(25.0);
+	ctx->image->tilt = deg2rad(90.0);
 	ctx->image->omega = deg2rad(90.0);
 	ctx->image->slop = deg2rad(0.0);
 	ctx->image->fmode = FORMULATION_PIXELSIZE;