testlist <- list(kt = 6.86175458771796e-165, X1 = structure(c(4.23557263671769e+195, 3.18466649861537e-103, 7.31623821489758e+231, 1.3893309810172e+294, 2.19737395714509e+120, 9.2595878678592e+301, 3.60082179961981e+73, 2.85854714673497e+162, 1.42543584611964, 9.96648148479037e+34, 3.86557585615097e-229, 3.85336711706673e+289, 4.17262911081223e+89, 9.5382644766677e-192, 1.68502814229884e+176, 8.32691060340956e-305, 2.35891172264555e-259, 4630.12425726243, 4.13491655520763e-66, 1.00143370341692e+149, 3.7469580425089e-227, 2.69300117472165e+236, 3.93093426192471e+192, 1.28714826811561e-211, 1.23641996983647e+269, 2.28776119140287e-200, 1.11356704442287e+73), .Dim = c(3L, 9L )), X2 = structure(c(1.83024402409462e+136, 5.50849136971673e-288, 4.30669059016305e-223, 4.46159392815599e+120, 2.88949314227748e-138, 7.86488204969752e+173, 2.98462041071517e-91, 2.90022730647209e-116, 6.27059809901118e-144, 1.07277125891001e-76, 1.98022092873024e+295, 55832.8564277184, 1.6385455562026e+118, 2.16870944877799e-12, 6.99283681783266e-155, 2.06656250657752e+245, 2.99279745548802e+131, 2.50653594286376e-164, 2.94218586521686e-25, 1.83124574864265e-19 ), .Dim = c(10L, 2L))) result <- do.call(hetGP:::d_matern5_2_2args_kthetag,testlist) str(result)