Added bitinformation.h

src/Bitrounding.cc

@@ -13,6 +13,7 @@
 #include "pmlist.h"
 #include "param_conversion.h"
 #include "progress.h"
+#include "bitinformation.h"
 
 struct BitroundParams
 {
@@ -32,16 +33,6 @@ struct VarStat
   int nsbMax = -1000;
 };
 
-#define NBITS 32
-  
-struct MutualInformation
-{
-  double M[NBITS] = { };
-};
-
-void signed_exponent(float *A, size_t n);
-MutualInformation bitinformation(float *A, size_t n);
-
 static int
 get_keepbits(MutualInformation &bitInfo, double inflevel)
 {

src/Makefile.am

@@ -7,6 +7,7 @@ libcdo_la_SOURCES =		after_dvtrans.cc          \
 				afterburner.h             \
 				afterburnerlib.cc         \
 				bitinformation.cc         \
+				bitinformation.h          \
 				cdi_lockedIO.cc           \
 				cdi_lockedIO.h            \
 				cdi_uuid.h                \

src/bitinformation.cc

@@ -18,18 +18,12 @@
 #endif
 
 #include "statistic.h" // cdo::normal_inv()
-
-constexpr int NBITS = 32; // Number of bits in type `float`
+#include "bitinformation.h"
 
 struct BitpairCounters
 {
   int C[NBITS][2][2] = { };
 };
-  
-struct MutualInformation
-{
-  double M[NBITS] = { };
-};
 
 /*
   p₁ = binom_confidence(n::Int,c::Real)

src/bitinformation.h

+#ifndef BITINFORMATION_H
+#define BITINFORMATION_H
+
+constexpr int NBITS = 32; // Number of bits in type `float`
+
+struct MutualInformation
+{
+  double M[NBITS] = { };
+};
+
+void signed_exponent(float *A, size_t n);
+MutualInformation bitinformation(float *A, size_t n);
+
+#endif