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
|
/*
* utils.c
*
* Utility stuff
*
* (c) 2006-2009 Thomas White <thomas.white@desy.de>
*
* pattern_sim - Simulate diffraction patterns from small crystals
*
*/
#include <math.h>
#include <string.h>
#include <stdio.h>
#include "utils.h"
/* Angle between two vectors. Answer in radians */
double angle_between(double x1, double y1, double z1,
double x2, double y2, double z2)
{
double mod1 = modulus(x1, y1, z1);
double mod2 = modulus(x2, y2, z2);
return acos( (x1*x2 + y1*y2 + z1*z2) / (mod1*mod2) );
}
size_t skipspace(const char *s)
{
size_t i;
for ( i=0; i<strlen(s); i++ ) {
if ( (s[i] != ' ') && (s[i] != '\t') ) return i;
}
return strlen(s);
}
void chomp(char *s)
{
size_t i;
if ( !s ) return;
for ( i=0; i<strlen(s); i++ ) {
if ( (s[i] == '\n') || (s[i] == '\r') ) {
s[i] = '\0';
return;
}
}
}
void progress_bar(int val, int total, const char *text)
{
double frac;
int n, i;
char s[1024];
const int width = 50;
frac = (double)val/total;
n = (int)(frac*width);
for ( i=0; i<n; i++ ) s[i] = '=';
for ( i=n; i<width; i++ ) s[i] = '.';
s[width] = '\0';
printf("\r%s: |%s|", text, s);
if ( val == total ) printf("\n");
fflush(stdout);
}
double poisson_noise(double expected)
{
double L;
int k;
double p = 1.0;
L = exp(-expected);
do {
double r;
k++;
r = (double)random()/RAND_MAX;
p *= r;
} while ( p > L );
return k - 1;
}
|
