testlist <- list(p = c(-4.43405758149991e+267, -1.3239616079755e-49, -1.28037702066458e-25,  2.06308058365872e+285, 1.47411360534017e+187, 1.02108350714091e-250,  -5.07223918751024e+195, 1.27789398170398e-261, 3.1149493670469e+30,  2.03377499380837e-295, -1.63465835470276e-257, 5.55856889378762e-198,  4.71098492449163e+56, -9.58604350664938e-76), theta = -6.40300377992993e+218,      x = c(-621726556L, -1331494887L, 1854878590L, 314985932L,      1361933119L, -1262797933L, -1571830262L, -555281357L, 1176038391L,      -1125445590L, 2051319948L, -555281357L, 600611799L, -2015925719L,      640758675L, -2045256755L, 818332908L, -1294723540L, 415152419L,      -1290633258L, -568676763L, -250049000L, 1237429200L, -645489399L,      1298265599L, 889924430L, -907775620L, -152022060L))
result <- do.call(ExactMultinom::multinom_test_cpp,testlist)
str(result)