testlist <- list(n = 644268439L, p = 123721430L, thetamat = structure(c(1.05414702709612e-68, 1.86446765666332e-211, 449058805123160512, 2.72550962279315e-214, 6.53679281395345e-70, 2.42643425911725e+241, 1.65686753146763e+24 ), .Dim = c(1L, 7L))) result <- do.call(flam:::calcsum,testlist) str(result)