/usr/include/pbseq/alignment/statistics/StatUtilsImpl.hpp is in libblasr-dev 0~20161219-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 | #ifndef _BLASR_STAT_UTILS_IMPL_HPP_
#define _BLASR_STAT_UTILS_IMPL_HPP_
#include <vector>
#include <time.h>
inline int Choose(int a, int b) {
// A quick choose for small values.
assert(a >= b);
if (a < FactorialTableLength and b < FactorialTableLength) {
return (FactorialTable[a] /
(FactorialTable[b] * FactorialTable[a-b]));
} else {
//
// Less quick choose for numerically stable values.
//
unsigned long al =a, bl=b, fact=0;
fact = al;
//
// The full value of b! must be able to be computed.
//
if (b >= FactorialTableLength) {
return 0;
}
while (al > bl) {
--al;
if (ULONG_MAX / al < fact) {
// there is overflow in the computation of a!/(a-b)!
return 0;
}
else {
fact = fact * al;
}
}
return fact / FactorialTable[b];
}
}
template<typename T>
void MeanVar(std::vector<T> &values, float &mean, float &var) {
T sum = 0;
T sumsq = 0;
size_t i;
if (values.size() == 0) {
mean = 0; var = 0;
return;
}
for (i = 0; i < values.size(); i++) {
sum += values[i];
sumsq += values[i]*values[i];
}
mean = (1.0*sum) / values.size();
var = sumsq / values.size() - (mean*mean);
}
#endif
|