testlist <- list(sigma = -1.30015477240878e-75, W = c(1.04181950239833e+246, -9.9710023554788e-183, 2.50968374253996e+259, NaN, NaN, -2.06282204313793e-63, -Inf, -2.2421709329862e-206, -84504.3030713751, Inf), x = -8.82200427882342e+72, X = c(NA, -11490814657114308608, -7.51550601878762e+234, -7.51550601878762e+234, NaN, -1.3314381453655e-138, -2.95191671611102e-222, -6.56582059567982e-301, 1.02435972144764e+277, -6.71856748114919e+63, -5.65535735314369e+111, 4.75813569835136e-257, -1.4954750043117e+32, -2.83574726323411e+224, 0)) result <- do.call(hetGP:::c1_mat32_cpp,testlist) str(result)