Реализация logsumexp на C?

Кто-нибудь знает о цифровой библиотеке C с открытым исходным кодом, которая предоставляет функцию logsumexp ?

Функция logsumexp (a) вычисляет сумму экспонент log (e ^ {a_1} + ... e ^ {a_n}) компонентов массива a, избегая числового переполнения.