slerp中cos_theta接近±1时需特殊处理，因theta精度丢失且sin(theta)趋近0导致除零；此时应转为lerp+归一化，并先做符号校正；t超范围需截断或取模以避免绕远。
SLERP 公式里 cos_theta 接近 ±1 时为什么必须特殊处理
因为直接套用标准 SLERP 公式 slerp(q0, q1, t) = (q0 * sin((1-t)*theta) + q1 * sin(t*theta)) / sin(theta) 时，当 q0 和 q1 几乎同向（cos_theta ≈ 1）或反向（cos_theta ≈ -1），theta = acos(cos_theta) 会失去精度，sin(theta) 趋近于 0，导致除零或数值震荡。
实际中更稳妥的做法是：当 abs(cos_theta) > 0.9995f（阈值可调）时，退化为线性插值（LERP）再归一化——不是因为“简单”，而是此时球面路径和直线路径在单位球面上几乎重合，LERP+normalize 的误差远小于数值崩溃的风险。
- 检查前先对
q1做符号校正：if (dot(q0, q1) ，确保走短弧 - 阈值选
0.9995f而非1.0f，避开浮点acos在端点的梯度爆炸区 - 退化分支中必须归一化，因为 LERP 结果一般不在单位球面上
用 std::acos 计算 theta 前要先 clamp cos_theta
std::acos 只接受 [-1, 1] 区间输入，但浮点运算误差常使点积结果略超范围（如 1.00000012 或 -1.00000012），直接传入会返回 NaN，后续全盘失效。
必须手动截断：
float cos_theta = dot(q0, q1); cos_theta = std::clamp(cos_theta, -1.0f, 1.0f); // C++17；C++14 可用 fmaxf(fminf(cos_theta, 1.0f), -1.0f) float theta = std::acos(cos_theta);
注意：不要用 abs(cos_theta) 替代 clamp——反向四元数需靠前面的符号校正处理，这里只保定义域安全。
手写四元数乘法和点积时，w 分量顺序不能错
多数数学库（如 glm、Eigen）按 (x, y, z, w) 存储，但部分图形 API（如 DirectX Math）或旧代码可能用 (w, x, y, z)。SLERP 实现中若点积写成 q0.x*q1.x + q0.y*q1.y + q0.z*q1.z + q0.w*q1.w，而实际内存布局是 wxyz，就会把 w 当成第一个分量参与计算，结果完全错误。
- 统一用命名成员访问（
q.x,q.y…）而非数组下标，避免索引偏移 - 点积函数单独抽离，加注释说明分量顺序假设
- 如果对接外部库（如
XMVECTOR），务必查文档确认其XMVectorGetX/Y/Z/W的映射规则
插值参数 t 超出 [0,1] 时的行为要明确控制
SLERP 本身支持 t 或 <code>t > 1（即外推），但多数动画系统期望边界截断。不处理会导致旋转“绕远”或突变——比如 t=1.1 时本该停在 q1，却继续往 q1 之后的方向走。
建议策略：
- 调试期保留外推能力，方便观察轨迹
- 生产环境默认
t = std::clamp(t, 0.0f, 1.0f) - 若需循环动画（如
t持续增长），应先做t = fmod(t, 1.0f)，再 clamp，否则fmod可能返回负余数
真正麻烦的是多段关键帧拼接时的切线连续性——SLERP 本身不保证角速度恒定，更别说高阶导数。需要平滑过渡得上 Squad 或样条拟合，那是另一层问题了。
C++免费学习笔记（深入）：立即使用
在学习笔记中，你将探索 C++ 的入门与实战技巧！
