So I run rev/fwd tests for scal/fun/gamma_p and they pass.
Then the mixed test fails. I added prints to the point that I can tell
bvi->adj_ is zero when the
gamma_p_vv_vari is constructed but by the time the
chain() method is called it has turned into
I can see that
gamma_p(fvar<var>, fvar<var>) gets called in between but everything I have access to and everything calculated in it is fine.
I'm not sure what else would get called in between those two points.
Any suggestions what I could look into? This is my first time digging into the mixed mode stuff so it's a steep learning curve. I made quite a few changes in getting the rev/fwd tests to pass but when I instrument those with prints everything comes out fine and scanning the changes I don't do anything exotic but I'm clearly missing something.