testlist <- list(I = c(NA, Inf, -1.09441869215336e+284, -5.64169870293727e-98, NaN, -Inf, 1.02522710329175e-279, 3.8030023523168e+24, NA), Mv = c(-Inf, -2.17099405775236e-33, -Inf, NA, NA), U1m = structure(c(1.08802942622396e-293, 6.65097473913016e-260, 9.93520077448275e+42, 1.08820604487168e+142, 1.74885170013241e+276, 3.68822443083645e-159, 1.48396048760851e-153, 3.92438875200275e-290, 1.01227910123508e-116, 7.76415668710948e-76, 3.0852705372561e+237, 2.88013821501196e-76, 3.83166867124718e+152, 6.49465960533195e-293, 2.34263610572941e-124, 6.63273697768014e-117, 2.33903696368172e-267, 1.15397297759478e+152, 5.59199246023954e+278, 7.43957891850816e+192, 1.72610087451397e-63, 1.63991633646346e+246, 2.48702195973594e-34, 1.05633173801447e-181, 7.89268557020173e-238, 5.48955018394356e+54, 3.69082222941437e-168, 5.36026497980026e-253, 1.41387477910997e-151, 5.12319547171943e-70, 6.05232307035114e-265, 4.45856749180119e+113, 8.01021890336983e+304, 4.79612622872741e-11, 6.29283118217647e+126, 1.71303509990658e-198, 4.81154875302327e+42, 6.13400938393855e-189, 5.72577012638986e+64, 4.51276261235999e+272, 4.2142340826561e-197, 1.05678586284277e-154, 1.9088637112238e+83, 1.26678222181216e+86, 3.91254592729957e+193, 3.64014866604686e+195, 8.24489692479263e-25, 6.08448796562516e-98, 8.13927445783608e+65, 1.46458388595101e+289, 1.18779940620107e-287, 374096.630477381, 1.42308123762198e-25, 1.98667977548534e-298, 9.00723533378294e-298, 2.08915060471778e+106, 3.88978931392951e+207, 8.86771684546731e+268, 2.17907772264764e-278, 4.78967389290926e+219, 4.93560832714263e-77, 3.99818599844377e+78, 7.62751612031997e+87, 5.52575563232104e-288, 1.14290399642585e-276, 2.61000445356826e-18, 6.11154982057764e-117, 4.31900718508964e-179, 2.83057484978328e-93, 9.35626819691237e-45, 8.01231174025142e+124, 6.08730340660145e+214), .Dim = 9:8), U2m = structure(c(6.55821997918708e-277, 9.60330917511895e-76, 5.61623761316617e+245, 4.67101646804148e-67, 1.96635951124372e+58, 3.80400403515519e-211, 3.94140934345603e+76, 1.60699812287607e+137, 2.76871756167236e+80, 1.99796556528912e-279, 9.58944044600216e-17, 6.70553219161172e+164, 5.65392662548674e+169, 1.16533480133947e-179, 1.96374930146633e+169, 4.16304518070251e+288 ), .Dim = c(4L, 4L)), Xm = structure(c(1.04677824199772e-134, 3.44329793513651e+299, 1.1629968206522e-18, 2.40926213705545e+109, 1.37104979282777e+308, 2.56265553086029e-82, 1.66130827925819e-208, 3.89441935319495e-145, 3.04402369120538e+65, 1.03273957420079e+210, 2.55658601588852e-308, 4.09882712914762e-305), .Dim = c(6L, 2L ))) result <- do.call(CMF:::p_updateMean,testlist) str(result)