- 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
0.1
- adds covariance to output - allows `__covariance__` to be passed as input - uses estimate Jacobian based on central finite difference method