Verlet积分+距离约束迭代（PBD）因跳过速度与力计算、仅修正位置而稳定高效；其核心是用prev_pos和curr_pos做显式预测，再通过多轮边投影强制归位，避免位置漂移导致的发散或抖动。
用 Verlet 积分 + 距离约束迭代（PBD）最直接
纯物理模拟里硬解牛顿方程容易发散，而布料对“质点间距离不准”又极度敏感。Verlet 积分天然稳定，配合距离约束的投影修正（Position-Based Dynamics），能快速收敛且不依赖步长。关键是：不求解速度，只修正位置，跳过质量、力、加速度的中间计算。
实操建议：
- 每个质点存
prev_pos和curr_pos，用curr_pos += (curr_pos - prev_pos) + dt * dt * gravity做显式 Verlet 预测 - 把每根弹簧（边）当作一个距离约束：目标长度是初始静止时的
rest_length - 每次迭代遍历所有边，对端点
p1、p2做一次投影：delta = p2 - p1，dist = length(delta)，若dist > 1e-6，则corr = (dist - rest_length) / dist * 0.5，再更新p1 -= delta * corr，p2 += delta * corr - 通常 3–8 次全局约束迭代就足够；更多次提升精度但边际收益递减
为什么不能只用单次弹簧力（Hooke 定律）
直接算力 → 加速度 → 更新速度 → 更新位置，会遇到三个硬伤：显式积分在大步长或软弹簧下爆炸；隐式需要解大型稀疏矩阵；而布料常含上百个质点+上千条边，实时性崩盘。
典型错误现象：
- 质点飞出屏幕，控制台没报错但
curr_pos值溢出到1e30 - 网格像果冻一样高频抖动，尤其在拉扯后回弹阶段
- 调高阻尼系数（
damping）后运动变粘滞，但抖动仍在——说明问题不在阻尼，而在位置漂移未被约束
根本原因：力模型放任位置误差累积，而距离约束是“强制归位”，属于几何修复而非动力学逼近。
边约束顺序和权重影响稳定性
约束迭代不是数学意义上的精确求解，顺序和权重决定收敛路径。无序遍历边会导致不同帧结果微异，动画闪烁；等权处理刚性边和边界边，会让固定点附近形变异常。
推荐做法：
- 按拓扑层级分组迭代：先处理固定顶点相连的边（如衣领锚点），再处理内部结构边，最后处理剪裁/缝合边
- 给固定点连接的边设更高约束强度（例如乘
2.0），即corr *= stiffness，其中stiffness是 per-edge 参数 - 避免在单次迭代中重复修改同一质点多次后才写回内存——用双缓冲数组或原子操作，否则出现“半更新”撕裂
GPU 加速时注意数据布局与访存模式
迁移到 CUDA 或 Vulkan Compute 时，常见性能坑不在 kernel 复杂度，而在内存访问。每个线程处理一条边，但 p1 和 p2 可能跨 cache line，反复读写同个质点引发 bank conflict。
关键调整：
- 质点数组必须是 AOS（
struct Particle { float x,y,z; float px,py,pz; }[]）而非 SOA，否则边迭代需 gather 多次 - 把边索引数组（
std::vector<:pair>>）转成两个 int 数组edge_i[]和edge_j[]，利于 coalesced read - 约束迭代轮数不宜上 GPU —— 通常只做 1 轮 GPU 投影 + CPU 补 2–3 轮，因分支少的 CPU 循环比 GPU warp divergence 更快
真正难的是多约束耦合：比如三角面片面积约束 + 边长约束 + 弯曲约束同时存在时，迭代顺序和 relaxation 系数得调半天。别指望一套参数通吃所有布料类型。
C++免费学习笔记（深入）：立即使用
在学习笔记中，你将探索 C++ 的入门与实战技巧！
