Added executable to compute Wasserstein distance.
Powering to 1/q is done in output. Maybe, should be changed in the
function itself?
#ifndef __NUMBER_TRAITS_H__
#define __NUMBER_TRAITS_H__
template<class NumberType_>
class number_traits
{
public:
typedef NumberType_ NumberType;
static NumberType& normalize(NumberType& n) { return n; }
};
#endif