#include <cmath>
#include "ld.h"

template<class T>
  constexpr T ln2_v
    = static_cast<T>(0.693147180559945309417232121458176568L);

double ld(double x) { return log(x) / ln2_v<double>; }