testlist <- list(i1 = -359125535L, i2 = 974295192L, X1 = structure(c(1.59289758025936e+240, 9.83154296510061e-219, 1.76824674457715e-27, 1.06070182584727e+60, 3.79268133289658e-62, 3.81747343765256e-296, 3.76389905463183e-48, 5.5813785128595e+31, 2.02743218632367e-260, 1.48207328366427e-297, 4.25923363672778e+136, 3.88481413819747e-32, 5.30335823633862e-292, 1.02264341624485e+291, 1.49384549398107e-118, 3.36222681779754e+59, 0.458508263044358, 3.25447937900225e-268, 40579.8806533551, 1.32229538094513e-288, 5.91258509337665e+97, 7.1824273432924e-163, 2.27125932461535e-51, 1.46435670316606e+63, 2.32499781549928e+103, 3.89444106107916e-44, 4.69101996929299e-85, 1.97430230093398e-107, 9.75098920301577e-204, 2.83363026718099e-284, 1.1714692424335e-190, 5.36415234459078e-39, 3.24344121758973e-153, 5.26807814693046e-92, 7.59949478474065e+62, 2.15491391334949e+38, 2.24163538980291e+116, 4.62052484558903e+132, 2.15330029890863e+35, 3.15833893127071e+27, 5.50444336166662e-289, 5.38681834396875e-269), .Dim = 6:7), X2 = structure(c(1.21884463874721e-233, 4.69454500785274e-52), .Dim = 1:2)) result <- do.call(hetGP:::partial_d_dist_dX1_i1_i2_X2,testlist) str(result)