powered by
matrix-ified version of logsum to avoid needing apply()
logsum_matrix(x)
numeric matrix
rowwise sums