Low-priority, but it would be nice if mvndst.f were modified to be thread-safe (so that mvn.pyf could release the GIL). I don't really know fortran, but at a glance this looks like a non-zero but not-ridiculous amount of work: there's a bit of cached data in COMMON and SAVE variables, but nothing that fundamentally needs to be single-threaded.