testlist <- list(theta = -2.52946195952421e+272, X1 = structure(c(1.71148770016929e-281, 1.46427674783955e+106, 1.3542881972378e+29, 3.91390114687515e-174, 5.5673644502049e-101, 7.38964498720712e+217, 1.36772606122782e+231, 1.18584814771421e-43, 1.04807931268534e-276, 1.21309430294063e+44 ), .Dim = c(1L, 10L)), X2 = structure(c(2.6274655749146e-112, 2.1627450010805e+106, 3.17209819179256e-24, 4.71441644263619e-27, 9.60785603800552e+25, 2.50882253234433e+271, 6.76295055483189e-185, 5.61562155049885e+188, 197567212180128, 1.15495685922267e-17, 1.50330193289056e+124, 1.24847770256251e-104, 1.9993322318126e+70, 4.60872253324786e-207, 2.47275449035714e+138, 2.12642390863418e-298, 7.1270469118521e+189, 8.14672456735504e+192, 1.73489904104165e-281, 2.70222976862149e+273, 6.14051220802483e+296, 30726.3674282008, 9.16380403656561e+64, 4.33006413061718e+242, 2.78664690945945e+265, 1.32901526215985e+287, 7.99188017639118e-230, 4.79328888514346e+109, 5.50380088837092e+243, 3.18053274134181e-172, 2.4420672371035e+230, 1.16562980053814e-39, 8.03576423005481e-269, 1.50879474723509e+132, 5.59694155795305e+256, 2.1040720961032e-220, 5.75412964504279e-19, 2.74929572251345e+92, 2.77961927026427e+293, 2.12175947511584e-82, 0.0308611642966599, 3.6127916780184e+133, 4.91401329653232e+234, 1.14078753704952e-55, 3.35967612129997e-123, 6.51052126746242e+213, 4.62038269452154e-145, 2.632502336894e-54, 1.5836458842161e+277, 7.90969629580482e+172, 2.5474821229675e+99, 5.67586778516112e-170, 3.52146147104697e+155, 6.60240777347802e+23, 6.04556062159806e+71, 5.08320729914457e+144, 2.92944783296954e-185, 4.81022598021046e-192, 1.35685606751012e+126, 2.55473483078402e+111, 4.88178903281954e-281, 2.44242688994123e+204, 6.28534284811382e+288, 8.11869628667905e+48, 9.91902940251281e+232, 0.0010484111775654, 4.26171479501292e-305, 2.72556404441845e-08, 4.32811565732498e+42, 9.17168661110121e-260, 1.36108442509638e+35, 2.12655370812234e+65, 7.17297304452373e+105, 1.24680702951631e+257, 1.08533257888063e-210, 3.93310951550094e+210, 1.46739066966128e+193, 5.11487409623387e-58, 8.27576892034027e-54, 6.90145472036555e-189), .Dim = c(10L, 8L))) result <- do.call(hetGP:::d_matern5_2_2args_theta_k_iso,testlist) str(result)