直接替换for循环常出错，因循环隐含顺序依赖、中间状态更新或条件分支，而NumPy向量化要求整块数组并行执行；强行用np.vectorize只是语法糖，性能更差；真正有效的是识别循环数学本质（映射/累积/掩码/查表等）并选用对应NumPy工具链。
为什么直接替换 for 循环常出错？
因为 for 循环隐含顺序依赖、中间状态更新或条件分支，而 NumPy 向量化要求所有操作在整块数组上并行执行。强行“把 for 换成 np.vectorize”不是真向量化，只是语法糖，性能甚至更差。
真正有效的重构，是识别原循环的数学本质：它是逐元素映射？累积计算？条件掩码？索引查表？还是滚动窗口？每种对应不同 NumPy 工具链。
逐元素计算：用广播 + 基础运算替代 for + append
常见于“对每个 x[i] 算一个表达式，存进 y[i]”。这时直接用数组运算，无需 np.vectorize 或 np.frompyfunc。
- 错误写法：
np.vectorize(lambda x: x**2 + 2*x + 1)(arr)—— 仍是 Python 函数调用，不加速 - 正确写法：
arr**2 + 2*arr + 1—— 利用 NumPy 广播，C 层实现 - 含条件时用
np.where：比如np.where(arr > 0, np.log(arr), 0)，避免for中if/else - 注意数据类型：
arr.astype(np.float64)再运算，防止整数除法截断或溢出
累积类逻辑：优先用 np.cumsum、np.cumprod、np.diff
循环中带 total += x[i] 或 prev = x[i] - prev 这类依赖前值的操作，不能简单广播，但多数有对应累积函数。
立即学习“Python免费学习笔记（深入）”；
-
for实现累加：res = []; s = 0; for v in arr: s += v; res.append(s)→ 直接np.cumsum(arr) - 反向累积（如后缀和）：
np.cumsum(arr[::-1])[::-1] - 一阶差分：
for i in range(1, len(arr)): diff.append(arr[i] - arr[i-1])→np.diff(arr) - 注意边界：这些函数默认 axis=0，多维时需显式指定
axis=1避免误操作
索引与查找类循环：用 np.searchsorted、布尔索引、np.take
循环遍历数组 A，在数组 B 中找匹配项、插值位置或按规则取值，这类最易卡在“不知道怎么摆脱 for”。
- 查找插入点（如二分查找）：
[np.searchsorted(boundaries, x) for x in values]→np.searchsorted(boundaries, values)（values可为数组） - 按条件取值：
[data[i] for i in indices if mask[i]]→data[mask]或data[indices[mask]] - 查表映射：
lookup = {0:10, 1:20}; [lookup[x] for x in keys]→ 先构造mapping = np.array([10,20])，再mapping[keys]（要求keys是合法整数索引） - 陷阱：
np.searchsorted要求boundaries有序；布尔索引返回副本，原地修改要用data[mask] *= 2这类赋值形式
最难的是嵌套循环或带状态的递推（比如卡尔曼滤波、动态规划），它们无法被单个 NumPy 函数替代，得用 Numba 加速或改写为 scipy.linalg 等专用算法。别硬向量化——先确认问题是否真的适合数组范式。
