Alpha diversity metrics more robust to edge cases
Better Brillouin calculations

1/ use gamma function instead of factorial to handle decimal counts
2/ use dedicated log(gamma()) function from Math::GSL::SF for speed