testlist <- list(Dx = structure(c(3.50419964461635e+307, 33464411539275352, 4.29892085338343e+147, 1.57238173865769e-23, 4.189376687691e+248, 2.17947419739359e+53, 2.59904368756192e-308, 5.65536745566999e-108, 5.04874306892903e+303, 2.25606437858434e-76, 1.59372707332277e-40, 3.07563721277347e-202, 1.9626495832104e-202, 1.46424493236939e+297, 4.69622589948052e-161, 2.45130177721813e-08, 4.60626137522111e-18, 3.52432344811232e-179, 1.17957554364568e-66, 2.81256767457907e+167 ), .Dim = c(10L, 2L))) result <- do.call(energy:::D_center,testlist) str(result)