Added executable to compute Wasserstein distance.
Powering to 1/q is done in output. Maybe, should be changed in the
function itself?
#include "grid2D.h"
#include <iostream>
int main()
{
Grid2D grid(40,60);
int i = 0;
for (int x = 0; x < 40; ++x)
for (int y = 0; y < 60; ++y)
{
grid(x,y) = i++;
}
std::cout << grid(20,30) << std::endl;
std::cout << grid << std::endl;
}