Solved by albanD in post #2 
              
              
                If I remember correctly, directional derivatives are Jacobian vector products. 
So you can compute them with the trick shown in this gist.