}
}
-/*
- * A(x - 0.5)^4 + B(x - 0.5)^3 + C(x - 0.5)^2 + D(x - 0.5)
- * A = 32,
- * B = -6
- * C = -7
- * D = 3
- * y = 32(x - 0.5)^4 - 6(x - 0.5)^3 - 7(x - 0.5)^2 + 3(x - 0.5)
- */
-static gdouble
-bayes_normalize_prob (gdouble x)
-{
- const gdouble a = 32, b = -6, c = -7, d = 3;
- gdouble xx, x2, x3, x4;
-
- xx = x - 0.5;
- x2 = xx * xx;
- x3 = x2 * xx;
- x4 = x3 * xx;
- return a*x4 + b*x3 + c*x2 + d*xx;
-}
gboolean
bayes_init (rspamd_mempool_t *pool, struct rspamd_classifier *cl)
* we need to rescale it to display correctly
*/
rspamd_snprintf (sumbuf, 32, "%.2f%%", (final_prob - 0.5) * 200.);
- final_prob = bayes_normalize_prob (final_prob);
+ final_prob = rspamd_normalize_probability (final_prob, 0.5);
g_assert (st != NULL);
cur = g_list_prepend (NULL, sumbuf);
rspamd_task_insert_result (task,
return map;
}
+
+/*
+ * A(x - 0.5)^4 + B(x - 0.5)^3 + C(x - 0.5)^2 + D(x - 0.5)
+ * A = 32,
+ * B = -6
+ * C = -7
+ * D = 3
+ * y = 32(x - 0.5)^4 - 6(x - 0.5)^3 - 7(x - 0.5)^2 + 3(x - 0.5)
+ */
+gdouble
+rspamd_normalize_probability (gdouble x, gdouble bias)
+{
+ const gdouble a = 32, b = -6, c = -7, d = 3;
+ gdouble xx, x2, x3, x4;
+
+ xx = x - bias;
+ x2 = xx * xx;
+ x3 = x2 * xx;
+ x4 = x3 * xx;
+
+ return a*x4 + b*x3 + c*x2 + d*xx;
+}
gpointer rspamd_shmem_xmap (const char *fname, guint mode,
gsize *size);
+/**
+ * Normalize probabilities using polynomial function
+ * @param x probability (bias .. 1)
+ * @return
+ */
+gdouble rspamd_normalize_probability (gdouble x, gdouble bias);
+
#define PTR_ARRAY_FOREACH(ar, i, cur) if (ar != NULL) for ((i) = 0, (cur) = g_ptr_array_index((ar), 0); (i) < (ar)->len; (cur) = g_ptr_array_index((ar), (i)), ++(i))
#endif