Tests

Author: vladtcvs

src/c_modules/bayes/libbayes/bayes.c

@@ -311,7 +311,6 @@ void bayes_posterior_free(struct bayes_posterior_ctx_s *ctx)
 
 bool bayes_posterior_init(struct bayes_posterior_ctx_s *ctx, int num_dim)
 {
-    bayes_posterior_free(ctx);
     ctx->N = num_dim;
     ctx->f_integration = calloc(num_dim, sizeof(double));
     ctx->f_estimation = calloc(num_dim, sizeof(double));

src/c_modules/bayes/libbayes/bayes.h

@@ -46,7 +46,7 @@ struct bayes_posterior_ctx_s
  * @param ctx - context
  * @param num_frames - number of measurements
  * @param F - frames {F_a}, a=0..num_frames-1
- * @param f - {f_i} for which we find p({p_i}|{F_a})
+ * @param f - {f_i} for which we find p({f_i}|{F_a})
  * @param d - constant coefficients {d_i}, for dark, sky
  * @param v - spatial coefficient {v_i}, for flat
  * @param K - contribution coefficients {K_i,a}

src/c_modules/bayes/libbayes/tests/bayes.cc

@@ -0,0 +1,31 @@
+#include <gtest/gtest.h>
+#include <bayes.h>
+
+extern "C" double 
+
+TEST(dim_1, aposteriori) {
+    struct bayes_posterior_ctx_s ctx;
+    bayes_posterior_init(&ctx, 1);
+
+    uint64_t F[] = {1, 1, 1, 1, 1};  // Samples
+    double f[1] = {1};          // target for which find aposteriori
+    double d[1] = {0};          // darks
+    double v[1] = {1};          // flats
+    double Ks[1] = {1};         // K
+    double lows[1] = {0};
+    double highs[1] = {10};
+    double dl = 0.1;
+
+    auto apriori = [](const double *f, int N, const void *param) -> double
+    {
+        if (f[0] < 0)
+            return 0;
+        if (f[0] > 10)
+            return 0;
+        return 1.0/10;
+    };
+
+    double posterior = bayes_posterior(&ctx, sizeof(F)/sizeof(F[0]), F, f, d, v, Ks, apriori, NULL, lows, highs, dl);
+
+    bayes_posterior_free(&ctx);
+}