dMdQ is now obtained via a second interpolation through points from the direct spline derivative (in order to reduce oscillations). Direct derivative remains callable via get_dMdQ_direct