PyMC3 use Arviz which for continuous value use a kde function to smooth out the histogram and get the local maximum:
It might not work well if you have multi mode but otherwise we think it is quite good in most case
2 Likes