It seems no need for c++ code. Calculating jacobians manually is time consuming. Can we automatically calculate it?