我可以使这一行矩阵乘法和加法代码更有效吗？

Question

这一行代码目前占用了大约50％的处理时间。有没有更有效的方法来处理这条线？

o = o + sum(sum(im.*p))*p/sum(sum(p)));

变量（ im ， p ＆amp; o ）都是大小为1000x1000（ish）的二维矩阵。

• ' im '是一个具有常量值的固定矩阵
• ' p '是在每次迭代时都会更改的稀疏矩阵
• ' o '是一个矩阵，在每次重置之后会重置 1000（ish）代码块的迭代。

如果您需要更多信息，请告诉我，或者我需要澄清我所说的任何内容。感谢您提前给出答案/反馈。