Automatic computation of three-loop two-point functions in large momentum expansion