<分区>
这是 C 代码:
btemp = (*beta)/(*beta_prev);
for (k=0; k<xsize*ysize; k++) {
parray[k] = zarray[k] + btemp*parray[k];
}
我正在用 Python 执行以下操作:
def saxpy(a, x, y):
return np.array([a * xi + yi for xi, yi in zip(x, y)], np.float32)
#...
btemp = beta / beta_prev
ptemp = saxpy(btemp, parray, zarray)
parray[:] = ptemp
在我的代码中,它似乎工作正常(zarray
和 parray
不断变化,因为它们在 while 循环中)。
但后来我做了:
btemp = beta / beta_prev
parray = saxpy(btemp, parray, zarray)
我的代码在循环中迭代几次后失败,它们不一样吗?