Compute unnormalised PSD using random mixture of B-splines
qpsd(omega, k, v, w, u, z, degree, recompute, db.list)