I wrote a few different calculations of
d(gamma_p)/da and I’d like to document where the different versions perform best and where I should use limits/approximations. I know I have at least one where the current version does as badly as 1.0e-3 for error and the new version is much better.
I was trying to modify Vincent’s Mathematica function for computing gradients. For
d(gamma_p(a,z)/da I’d like to get it to return the point the gradient was calculated at along with the gradient itself to make testing less confusing in C++. So the table I’m hoping for has the columns:
a z d(gamma_p)/da d(gamma_p)/dz ... ... .... ...
Any chance I could get a hand with getting this to work? So far I’ve only managed to screw it up.