diff options
author | taw27 <taw27@bf6ca9ba-c028-0410-8290-897cf20841d1> | 2008-01-18 18:24:18 +0000 |
---|---|---|
committer | taw27 <taw27@bf6ca9ba-c028-0410-8290-897cf20841d1> | 2008-01-18 18:24:18 +0000 |
commit | 61f8a2450fd944920523c0153ea1702f72d751fb (patch) | |
tree | ab09c913ac6208863f320e3d3260f655716d8110 | |
parent | 3a2ad78bf87b86ccc9a8ee7856a3ffbc261dcb7e (diff) |
Refinement diagnostics...
git-svn-id: svn://cook.msm.cam.ac.uk:745/diff-tomo/dtr@259 bf6ca9ba-c028-0410-8290-897cf20841d1
-rw-r--r-- | src/iprtest.c | 24 |
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; |