aboutsummaryrefslogtreecommitdiff
path: root/src/reflections.h
blob: 8fafe4bdade15cf9f41e4a8c5e7dc47d365379f6 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
/*
 * reflections.h
 *
 * Data structures in 3D space
 *
 * (c) 2007-2008 Thomas White <taw27@cam.ac.uk>
 *
 *  dtr - Diffraction Tomography Reconstruction
 *
 */
 
#ifndef REFLECTION_H
#define REFLECTION_H

#ifdef HAVE_CONFIG_H
#include <config.h>
#endif

typedef enum {
	REFLECTION_NORMAL,	/* Measured - expressed as x, y, z position */
	REFLECTION_GENERATED,	/* Generated and intensity-measured - expressed as h, k, l index */
	REFLECTION_MARKER,
	REFLECTION_VECTOR_MARKER_1,
	REFLECTION_VECTOR_MARKER_2,
	REFLECTION_VECTOR_MARKER_3
} ReflectionType;

typedef struct reflection_struct {

	double		x;
	double		y;
	double		z;
	double 		intensity;

	signed int	h;
	signed int	k;
	signed int	l;

	ReflectionType	type;
	
	/* Stuff used when finding bases */
	int		found;		/* This reflection has been used in the seed-finding process */
	int		lfom;		/* Line FoM for this seed */

	struct reflection_struct *next;	/* MUST BE LAST in order for caching to work */

} Reflection;

typedef struct reflectionlist_struct {

	Reflection *reflections;
	Reflection *last_reflection;
	
	unsigned int n_reflections;
	unsigned int list_capped;
	
} ReflectionList;

extern ReflectionList *reflectionlist_new(void);
extern void reflectionlist_clear(ReflectionList *reflectionlist);
extern void reflectionlist_clear_markers(ReflectionList *reflectionlist);
extern void reflectionlist_free(ReflectionList *reflectionlist);

extern Reflection *reflection_add(ReflectionList *reflectionlist, double x, double y, double z, double intensity, ReflectionType type);
extern int reflection_is_easy(Reflection *reflection);

extern double reflectionlist_largest_g(ReflectionList *reflectionlist);
extern Reflection *reflectionlist_find_nearest(ReflectionList *reflectionlist, double x, double y, double z);
extern Reflection *reflectionlist_find_nearest_longer_unknown(ReflectionList *reflectionlist, double x, double y, double z, double min_distance);
extern Reflection *reflectionlist_find_nearest_type(ReflectionList *reflectionlist, double x, double y, double z, ReflectionType type);
extern Reflection *reflectionlist_find(ReflectionList *reflectionlist, signed int h, signed int k, signed int l);

#include "basis.h"
extern ReflectionList *reflection_list_from_cell(Basis *basis);
extern void reflection_list_from_new_cell(ReflectionList *ordered, Basis *basis);

#endif	/* REFLECTION_H */