From 6255a45e28b32ee8116418d849db68395ec4e096 Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Fri, 5 Dec 2025 10:46:23 +0100
Subject: Compute maf again

---
 src/mathfunc.h | 5 +++++
 1 file changed, 5 insertions(+)

(limited to 'src/mathfunc.h')

diff --git a/src/mathfunc.h b/src/mathfunc.h
index 641d0a3..b408ea7 100644
--- a/src/mathfunc.h
+++ b/src/mathfunc.h
@@ -24,6 +24,7 @@
 // #include "Eigen/Dense"
 #include "gsl/gsl_matrix.h"
 #include "gsl/gsl_vector.h"
+#include <tuple>
 
 #define CONDITIONED_MAXRATIO 2e+6 // based on http://mathworld.wolfram.com/ConditionNumber.html
 #define EIGEN_MINVALUE 1e-10
@@ -77,4 +78,8 @@ unsigned char Double02ToUchar(const double dosage);
 // void uchar_matrix_get_row(const vector<vector<unsigned char>> &X,
 //                          const size_t i_row, Eigen::VectorXd &x_row);
 
+std::tuple<size_t, size_t, size_t> compute_ratio(size_t ni_total, const gsl_vector *gs);
+double compute_maf(size_t ni_total, size_t ni_test, size_t n_miss, const double *gs);
+
+
 #endif
-- 
cgit 1.4.1