maybe we should package this up and submit it to scipy? It could be a 'low-noise' flag in the next scipy release.