aboutsummaryrefslogtreecommitdiff
path: root/src/likelihood.c
blob: 417ebea3d3709c6f7049834fc0a074c30526eb6c (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
/*
 * likelihood.c
 *
 * Likelihood maximisation
 *
 * (c) 2006-2010 Thomas White <taw@physics.org>
 *
 * Part of CrystFEL - crystallography with a FEL
 *
 */


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

#include "statistics.h"
#include "utils.h"

void detwin_intensities(const double *model, double *new_pattern,
                        const unsigned int *model_counts,
                        ReflItemList *items)
{
	/* Placeholder... */
}

void scale_intensities(const double *model, double *new_pattern,
                       const unsigned int *model_counts,
                       ReflItemList *items, double f0, int f0_valid)
{
	double s;
	unsigned int i;
	unsigned int *new_counts;

	new_counts = items_to_counts(items);

	s = stat_scale_intensity(model, model_counts, new_pattern, new_counts);
	if ( f0_valid ) printf("%f %f\n", s, f0);

	/* NaN -> abort */
	if ( isnan(s) ) return;

	/* Multiply the new pattern up by "s" */
	for ( i=0; i<LIST_SIZE; i++ ) {
		new_pattern[i] *= s;
	}

	free(new_counts);
}