poly-dk-solve.h

#include "poly.h"
#include <complex>

std::vector<std::complex<double> > 
DK(Poly const & ply, const double tol=1e-10);