aboutsummaryrefslogtreecommitdiff
path: root/src/templates.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/templates.c')
-rw-r--r--src/templates.c183
1 files changed, 0 insertions, 183 deletions
diff --git a/src/templates.c b/src/templates.c
deleted file mode 100644
index 1ad1c95a..00000000
--- a/src/templates.c
+++ /dev/null
@@ -1,183 +0,0 @@
-/*
- * templates.c
- *
- * Handle templates
- *
- * (c) 2006-2009 Thomas White <thomas.white@desy.de>
- *
- * pattern_sim - Simulate diffraction patterns from small crystals
- *
- */
-
-
-#define _GNU_SOURCE 1
-#include <stdlib.h>
-#include <assert.h>
-#include <math.h>
-#include <stdio.h>
-
-#include "templates.h"
-#include "cell.h"
-#include "image.h"
-#include "utils.h"
-#include "relrod.h"
-
-
-struct _templatelist
-{
- int n_templates;
- struct template *templates;
-};
-
-
-struct template_feature
-{
- float x;
- float y;
- struct template_feature *next;
-};
-
-
-struct template
-{
- float omega;
- float tilt;
- struct template_feature *features;
-};
-
-
-static int template_add(TemplateList *list, struct template *template)
-{
- if ( list->templates ) {
- list->templates = realloc(list->templates,
- (list->n_templates+1)*sizeof(struct template));
- } else {
- assert(list->n_templates == 0);
- list->templates = malloc(sizeof(struct template));
- }
-
- /* Copy the data */
- list->templates[list->n_templates] = *template;
-
- list->n_templates++;
-
- return list->n_templates - 1;
-}
-
-
-static TemplateList *template_list_new()
-{
- TemplateList *list;
-
- list = malloc(sizeof(TemplateList));
-
- list->n_templates = 0;
- list->templates = NULL;
-
- return list;
-}
-
-
-TemplateList *generate_templates(UnitCell *cell, struct image params)
-{
- TemplateList *list;
- double omega, tilt;
-
- list = template_list_new();
-
- omega = deg2rad(40.0);
-
- //for ( omega=deg2rad(-180); omega<deg2rad(180); omega+=deg2rad(1) ) {
-
- params.omega = omega;
-
- for ( tilt=0; tilt<deg2rad(180); tilt+=deg2rad(1) ) {
-
- struct template t;
- struct template_feature *tfc;
- int nrefl, i;
-
- t.omega = omega;
- t.tilt = tilt;
- t.features = malloc(sizeof(struct template_feature));
- t.features->next = NULL;
- tfc = t.features;
-
- params.tilt = tilt;
- get_reflections(&params, cell, 0.01e9);
-
- nrefl = image_feature_count(params.rflist);
- for ( i=0; i<nrefl; i++ ) {
-
- struct imagefeature *f;
-
- f = image_get_feature(params.rflist, i);
-
- tfc->x = f->x;
- tfc->y = f->y;
- tfc->next = malloc(sizeof(struct template_feature));
- tfc = tfc->next;
-
- }
-
- template_add(list, &t);
-
- image_feature_list_free(params.rflist);
-
- printf("Generating templates... %+5.2f %+5.2f\r",
- rad2deg(omega), rad2deg(tilt));
-
- }
- //}
-
- printf("Generating templates... done \n");
-
- return list;
-}
-
-
-static int try_template(struct image *image, struct template template)
-{
- int fit = 0;
- struct template_feature *f;
-
- f = template.features;
- while ( f != NULL ) {
-
- int x, y;
-
- x = f->x;
- y = f->y; /* Discards digits after the decimal point */
-
- fit += image->data[y*image->width+x];
-
- f = f->next;
-
- }
-
- return fit;
-}
-
-
-int try_templates(struct image *image, TemplateList *list)
-{
- int i;
- int fit_max = 0;
- int i_max = 0;
-
- for ( i=0; i<list->n_templates; i++ ) {
-
- int fit;
-
- fit = try_template(image, list->templates[i]);
- if ( fit > fit_max ) {
- fit_max = fit;
- i_max = i;
- }
-
- }
- image->omega = list->templates[i_max].omega;
- image->tilt = list->templates[i_max].tilt;
-
- return 0;
-}