add callback function to SolveLiouvilleEq