Invert jacobian in loopy

Do not call a c++ function for calculation of the inverse of the jacobian and the determinant but implement it directly in loopy.