@inbounds for Dirichlet estimation