射线与球体相交的快速解法是将射线方程代入球面方程得二次方程at²+bt+c=0，其中a=dot(d,d)，b=2·dot(d,oc)，c=dot(oc,oc)−r²；判别式决定交点数量，t>1e-4才有效，需防零向量和浮点误差。
射线与球体相交的数学原理怎么快速落地
直接用向量代数解二次方程，别手推公式——核心是把射线 r(t) = origin + t * direction 代入球面方程 |p - center|² = radius²，整理后得到关于 t 的标准二次方程 a*t² + b*t + c = 0。其中：
-
a = dot(direction, direction)（通常归一化后为 1） b = 2.0f * dot(direction, origin - center)c = dot(origin - center, origin - center) - radius * radius
判别式 discriminant = b*b - 4.0f*a*c 决定是否相交：小于 0 无交点；等于 0 相切（一个交点）；大于 0 有两个交点（进入和离开球体）。
为什么 t 值必须大于 0 才算有效交点
射线是有方向、从起点单向延伸的，t 意味着交点在射线起点“后面”，属于反向延长线，不是实际射线的一部分。很多初学者漏判这个条件，导致检测到“身后”的球体。
- 即使
discriminant >= 0，也必须检查t1 > 0或t2 > 0 - 若只关心是否相交（不区分进出），取较小正根即可：
if (t1 > 0.0f) return true;或if (t2 > 0.0f) return true; - 若需最近交点，选两个正根中更小的那个；若都为负，返回 false
C++ 实现时要注意浮点精度和零向量陷阱
direction 未归一化时，a 不为 1，但不影响求根逻辑；真正危险的是 direction 几乎为零向量——此时 a ≈ 0，二次方程退化为一次方程，直接除零或误判。
立即学习“C++免费学习笔记（深入）”；
- 加一道前置检查：
if (a （避免数值不稳定） - 计算
t时推荐用稳定公式：t1 = (-b - sqrt(discriminant)) / (2.0f * a)，t2 = (-b + sqrt(discriminant)) / (2.0f * a) - 比较
t与 0 时用t > 1e-4f而非t > 0.0f，防止浮点误差导致本该命中却判定为错失
bool intersectRaySphere(const vec3& origin, const vec3& direction,
const vec3& center, float radius) {
vec3 oc = origin - center;
float a = dot(direction, direction);
if (a < 1e-8f) return false;
float b = 2.0f * dot(direction, oc);
float c = dot(oc, oc) - radius * radius;
float discriminant = b*b - 4.0f*a*c;
if (discriminant < 0.0f) return false;
float sqrt_d = sqrt(discriminant);
float t1 = (-b - sqrt_d) / (2.0f * a);
float t2 = (-b + sqrt_d) / (2.0f * a);
return (t1 > 1e-4f) || (t2 > 1e-4f);
}要不要提前归一化 direction
归一化能让 a == 1，简化计算并让 t 具有真实距离意义（单位长度），但每次调用都做归一化有开销。实际工程中更常见的是：如果射线方向已知是单位向量（如相机视线、归一化后的方向向量），就跳过归一化；否则，在调用方保证输入合法，或在函数内加 assert。
- 不做归一化：性能略优，
t是沿原始方向的参数距离，非欧氏距离 - 做归一化：语义清晰，
t就是真实距离，适合需要精确距离的场景（如阴影计算、碰撞响应） - 别在函数里无条件
normalize(direction)——它会修改原值或产生临时对象，徒增开销
真正容易被忽略的不是公式本身，而是对 t 的物理解释和浮点边界处理。写完先用原点在球心、方向直指球心的 case 测一遍，再试个 origin 在球内的情况——这时候 t1 ，应返回 true。
