testlist <- list(Mu1 = structure(c(4.38301413865151e-35, 2.8545821110746e-303, 358321389985497, 2.15202335026903e+42, 3.09544053665359e-229, 1.59383505298191e-118, 3.38579670155587e-138, 5.81679655092294e+295, 9.81465984779599e+122, 7.70799944571642e-280, 2.36180070899938e-148, 2.20063746926306e-84, 2.5263837895969e-299, 1.32750227584696e+85, 114449888591.725, 8486255.23406358, 1.7741571488086e+248, 5.19157220070939e-131, 7.80304816637635e-256, 2.2636440796552e-39, 1.49654517175992e+156, 3.36808417958313e-300, 1.31525791543e+104, 40637425520827891712, 2.41969872196061e+37, 3.00168066547143e-179, 1.20262486711107e+263, 3.97540886654159e-202, 1.37409419429572e-135, 0.0191068285870576, 8.50834915575324e-198, 2.32906384625588e-33, 4.48667040887042e+284, 4.9807763814948e-119, 4.30975676360895e+100, 5.7280198428626e+246, 1.10003333878246e-175, 4.10936008776511e-243, 7.12949415444687e-239, 1.68501758319177e-135, 2.21691734365394e-12, 3.08048967198988e+133, 4.85406985049335e-292, 9.1605854201672e+296, 4.7089102509115e+55, 1.49076689588462e-204, 1.00082758517919e+280, 9.03039594543831e-225, 26.2101255832183, 4.8399371213227e+86), .Dim = c(5L, 10L)), Mu2 = structure(c(6.24661546491917e+234, 1.77217125239726e+156, 6.19155847346587e+279, 2.97365548629946e+76, 2.00855676413494e-69, 1.7287788908324e-23, 9.05874833753716e-61, 4.58007613797977e+268, 7.95618330887762e-193, 4.75697157357305e-121, 7.79594966458731e+204, 6.87012279834584e+31, 6.02776133992854e+40, 6.96925489492716e-237, 2.08583575738941e-257, 2.87876358173649e+97, 6.27799658395281e-76, 2.69499027827958e+153, 2.35645726499149e+203, 1.51309052228304e+43, 1.06057783067153e-308, 2.35764783317463e-118, 3.67303558308306e-08, 3.65583966135793e+235, 9.86182652201338e-207, 7.38693948562522e+149, 3.62353771859118e+173, 3.55612572154597e-97, 5.33028074247052e-145, 9.69220548064053e+305, 5080902604941.94, 9.94199731557921e+50, 5.91539853818402e+153, 1.48571100223794e-249, 2.68921642919075e+70, 2.31547858529958e-197, 2.89451914904983e-50, 6.81765982188645e-65, 8.70937205221206e+132, 1.47001498138747e-227, 5.9419321090095e-84, 8.22352033542228e+195, 7.6739888029008e-232, 2.05529509741212e+222, 4.69581533203528e-31, 6.59713702361519e+244, 5.66911922652848e+83, 7.22739118489752e-50, 5.62070074749708e-128 ), .Dim = c(7L, 7L)), sigma = c(9.83043909548236e-82, NaN, NA, -1.50475398751523e+230, -3.2152055911576e+249, -2.89564090217024e-89, Inf, Inf, -4.29636395053306e-238, -2.21061027048667e-49, -2.27061266389401e-188, 0)) result <- do.call(hetGP:::Wijs_mat32_cpp,testlist) str(result)