Archean Eon
- Fix nargs and nf order mixup in Jacobian
- add more complex test
- fix tile cov by nobs
- move partial derivative to subfunction
- try threading, but same speed, and would only work with NumPy anyway
- TODO: try pthreads or winthreads if safe to release Python GIL