Jacobian_inverse function returns J inverse transpose.

[jacob-moore22]

So, two options, change the function name to note that it returns the transpose, or modify the function so it returns the transpose of the transpose, so J^{-1}.