R：应用于nXm矩阵的每一列，返回n*n*m数组

在矩阵/nXm数组的每一列上应用一个返回n*n矩阵的函数，并返回一个nXm数组，合适的函数是什么？我可以让它与列表一起工作：

dm <- structure(c(0.205292439079469, 0.222642769676055, 0.222786614980172, 
              0.210917067409095, 0.123955673451174, 0.10591537361648, 0.0887665882561804, 
              0.0725263089185602, 0.174088645938512, 0.171936521036433, 0.159619467190219, 
              0.143492610130578, 0.23303504847129, 0.269389294469438, 0.279103519495381, 
              0.240824583373471, 0.192637398624569, 0.201551950015116, 0.195748345926477, 
              0.181448686311511), .Dim = 4:5)
sapply(1:ncol(dm), function(x){dm[, x] %*% t(dm[, x])^-1}, simplify = F)

dm使用
simplify=“array”
而不是
simplify=F

sapply(1:ncol(dm), function(x){dm[, x] %*% t(dm[, x])^-1}, simplify = "array")