In case anyone finds it useful, here's the laguerre function as a numpy extension.
https://github.com/rm-rf-etc/ehlers

It's much faster. I'll try to get this published as a module so others can get it using pip.