如果M是一個密集的m×n矩陣並且v是一個n分量向量,那麼產品u = Mv是由u[i] = sum(M[i,j] * v[j], 1 <= j <= n)給出的m分量向量。一個簡單的實現該乘法是 allocate m-component vector u of zeroes
for i = 1:m
for j = 1:n
u[i] += M[i,j] * v[j]
e
我有一段可怕的時間試圖做這個看似簡單的函數。我不明白爲什麼這不起作用。 我的代碼: C = matrix(0 , dim(matrixA)[1] , dim(matrixB)[2] , TRUE)
mat_mult = function(matrixA, matrixB){
for (ii in 1:dim(matrixA)[1]) {
for (jj in 1:dim(matrix