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 Thomas White <taw27@cam.ac.uk>
* dtr - Diffraction Tomography Reconstruction
*
*/
#ifdef HAVE_CONFIG_H
#include <config.h>
#endif
#ifndef REFLECTION_H
#define REFLECTION_H
typedef enum {
REFLECTION_NORMAL, /* Measured - expressed as x, y, z position */
REFLECTION_CENTRAL, /* The central beam */
REFLECTION_GENERATED /* Generated and intensity-measured - expressed as h, k, l index */
} ReflectionType;
typedef struct reflection_struct {
double x;
double y;
double z;
double intensity;
signed int h;
signed int k;
signed int l;
ReflectionType type;
struct reflection_struct *next; /* MUST BE LAST in order for caching to work */
} Reflection;
typedef struct {
double x;
double y;
double z;
double modulus; /* Convenience */
} Vector;
typedef struct {
Vector a;
Vector b;
Vector c;
} Basis;
typedef struct rctx_struct {
Reflection *reflections;
Reflection *last_reflection;
Basis *basis;
} ReflectionContext;
extern ReflectionContext *reflection_init(void);
void reflection_clear(ReflectionContext *reflectionctx);
extern void reflection_add(ReflectionContext *reflectionctx, double x, double y, double z, double intensity, ReflectionType type);
extern void reflection_add_index(ReflectionContext *reflectionctx, signed int h, signed int k, signed int l, double intensity, ReflectionType type);
#include "control.h"
extern void reflection_add_from_dp(ControlContext *ctx, double x, double y, double tilt_degrees, double intensity);
extern void reflection_add_from_reciprocal(ControlContext *ctx, double x, double y, double tilt_degrees, double intensity);
extern void reflection_add_from_reflection(ReflectionContext *rctx, Reflection *r);
#endif /* REFLECTION_H */
|
