working likelihood function in scipy_dace