I'm working on adding automatic differentiation to `delicatessen` (to compute the variance exactly instead of approximating it). Still a work in progress but if you have time to test it out, it would help me out a lot

https://github.com/pzivich/Delicatessen/pull/24