How to model very complicated likelihood automatically?

Finally, I choose to use approximate bayesian computation (ABC), like using kernel density estimate (KDE) to calculate likelihood approximately.