powered by
Compute a PSD in the Bernstein-Dirichlet parametrization.
qpsd(omega, v, w, k, beta_basis_k, epsilon = 1e-20)