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
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
|
/*
* symmetry.c
*
* Symmetry
*
* (c) 2006-2010 Thomas White <taw@physics.org>
*
* Part of CrystFEL - crystallography with a FEL
*
*/
#ifdef HAVE_CONFIG_H
#include <config.h>
#endif
#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include "utils.h"
/* Conditions for a reflection to be in the asymmetric unit cell */
#define COND_1(h, k, l) (1)
#define COND_6MMM(h, k, i, l) ( (h>=0) && (k>=0) && (l>=0) && ((h>k)||((h==0)&&(k==0))) )
/* Macros for checking the above conditions and returning if satisfied */
#define CHECK_COND_FOURIDX(h, k, i, l, cond) \
if ( COND_##cond((h), (k), (i), (l)) ) { \
*hp = (h); *kp = (k); *lp = (l); \
return; \
}
#define CHECK_COND_THREEIDX(h, k, l, cond) \
if ( COND_##cond((h), (k), (l)) ) { \
*hp = (h); *kp = (k); *lp = (l); \
return; \
}
/* Abort macro if no match found */
#define SYM_ABORT \
ERROR("No match in %s for %i %i %i\n", sym, h, k, l); \
abort();
void get_asymm(signed int h, signed int k, signed int l,
signed int *hp, signed int *kp, signed int *lp,
const char *sym)
{
if ( strcmp(sym, "1") == 0 ) {
CHECK_COND_THREEIDX(h, k, l, 1);
SYM_ABORT;
}
if ( strcmp(sym, "6/mmm") == 0 ) {
const signed int i = h+k;
CHECK_COND_FOURIDX(h, k, i, l, 6MMM);
CHECK_COND_FOURIDX(h, i, k, l, 6MMM);
CHECK_COND_FOURIDX(k, h, i, l, 6MMM);
CHECK_COND_FOURIDX(k, i, h, l, 6MMM);
CHECK_COND_FOURIDX(i, h, k, l, 6MMM);
CHECK_COND_FOURIDX(i, k, h, l, 6MMM);
CHECK_COND_FOURIDX(h, k, i, -l, 6MMM);
CHECK_COND_FOURIDX(h, i, k, -l, 6MMM);
CHECK_COND_FOURIDX(k, h, i, -l, 6MMM);
CHECK_COND_FOURIDX(k, i, h, -l, 6MMM);
CHECK_COND_FOURIDX(i, h, k, -l, 6MMM);
CHECK_COND_FOURIDX(i, k, h, -l, 6MMM);
CHECK_COND_FOURIDX(-h, -k, -i, l, 6MMM);
CHECK_COND_FOURIDX(-h, -i, -k, l, 6MMM);
CHECK_COND_FOURIDX(-k, -h, -i, l, 6MMM);
CHECK_COND_FOURIDX(-k, -i, -h, l, 6MMM);
CHECK_COND_FOURIDX(-i, -h, -k, l, 6MMM);
CHECK_COND_FOURIDX(-i, -k, -h, l, 6MMM);
CHECK_COND_FOURIDX(-h, -k, -i, -l, 6MMM);
CHECK_COND_FOURIDX(-h, -i, -k, -l, 6MMM);
CHECK_COND_FOURIDX(-k, -h, -i, -l, 6MMM);
CHECK_COND_FOURIDX(-k, -i, -h, -l, 6MMM);
CHECK_COND_FOURIDX(-i, -h, -k, -l, 6MMM);
CHECK_COND_FOURIDX(-i, -k, -h, -l, 6MMM);
SYM_ABORT; /* Should never reach here */
}
ERROR("Unknown point group '%s'\n", sym);
abort();
}
|