↰ Return to documentation for file (include/networkit/auxiliary/MissingMath.hpp)
/*
* MissingMath.hpp
*
* Created on: 21.03.2013
* Author: cls
*/
#ifndef NETWORKIT_AUXILIARY_MISSING_MATH_HPP_
#define NETWORKIT_AUXILIARY_MISSING_MATH_HPP_
#include <cassert>
#include <cinttypes>
#include <cmath>
#include <stdexcept>
namespace Aux {
namespace MissingMath {
inline int64_t binomial(int64_t n, int64_t k) {
if (k == 0)
return 1;
if (2 * k > n)
return binomial(n, n - k);
int64_t b = n - k + 1;
for (int64_t i = 2; i <= k; ++i) {
b = b * (n - k + i);
b = b / i;
}
return b;
}
inline double log_b(double x, double b) {
if (x == 0) {
throw std::domain_error("log(0) is undefined");
}
assert(std::log(b) != 0);
return std::log(x) / std::log(b);
}
} /* namespace MissingMath */
} /* namespace Aux */
#endif // NETWORKIT_AUXILIARY_MISSING_MATH_HPP_