Added implementation of second derivatives of DR bar stop masses w.r.t. vu