From e447c2efb128823c93358c51a7d8f23636740f68 Mon Sep 17 00:00:00 2001
From: Thomas White <taw@bitwiz.org.uk>
Date: Tue, 4 May 2010 07:38:31 -0700
Subject: process_hkl: Implement --scale option

---
 src/likelihood.c | 44 ++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 44 insertions(+)
 create mode 100644 src/likelihood.c

(limited to 'src/likelihood.c')

diff --git a/src/likelihood.c b/src/likelihood.c
new file mode 100644
index 00000000..e8c62010
--- /dev/null
+++ b/src/likelihood.c
@@ -0,0 +1,44 @@
+/*
+ * 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,
+                        unsigned int *new_counts)
+{
+	/* Placeholder... */
+}
+
+void scale_intensities(const double *model, double *new_pattern,
+                       const unsigned int *model_counts,
+                       unsigned int *new_counts)
+{
+	double s;
+	unsigned int i;
+
+	s = stat_scale_intensity(model, model_counts, new_pattern, new_counts);
+	printf("%f\n", s);
+
+	/* NaN -> abort */
+	if ( isnan(s) ) return;
+
+	/* Multiply the new pattern up by "s" */
+	for ( i=0; i<LIST_SIZE; i++ ) {
+		new_counts[i] *= s;
+	}
+}
-- 
cgit v1.2.3