阵列产品的矢量化

Question

如何对此循环进行矢量化？

让# example 1 tryCatch( hello_world(), message = function(x) { cat("the message is: ", x$message, "\n") } ) # example 2 tryCatch( hello_world(), message = function(x) { cat("the message is: ", x$message, "\n") hello_world() } ) =真实的大小数组（2,3）

b =真实的大小数组（7,2）

M =真实的大小数组（n，3）

y和b是不变的，M是＆＃34;大＆＃34; （在我的问题中，n> 80000）。

我想优化（矢量化）这个循环：

n

帮助将不胜感激

Answer 1

首先p（K）= r＆＃39; * M错误

（b * y（k，:)）＆＃39;是（2X1），M是（7X2））

所以我认为你的意思是p（k）= r＆＃39; * M＆＃39 ;;

然后

b是（2X3）

是（nX3）

r = b * y＆＃39;是（2Xn）

R＆＃39;是（nX2）

M＆＃39;是（2X7）

P = R＆＃39; * M＆＃39;是（nX7）

所以你不需要使用for循环，只需写：

P =（B * Y＆＃39;）＆＃39; M＆＃39; 要么 P =（M （B * Y＆＃39;））＆＃39;