BVH 比暴力遍历快得多，因其通过空间聚类构建二叉树，利用AABB剔除大幅减少相交检测次数，实测百万面模型单次射线求交从2–5ms降至0.03–0.1ms；关键在高剔除率而非建树速度。
为什么 BVH 比暴力遍历快得多
射线与三角形逐个求交，时间复杂度是 O(N)，N 是三角形数量。一个 100 万面的网格，单次射线测试就要做百万次浮点运算和分支判断，实际渲染或交互中根本扛不住。BVH 把三角形按空间聚类分层组织成二叉树，每次射线遍历先快速剔除整个子树（比如用 AABB::intersect 判断射线是否与包围盒相交），平均下来只需检测几十到几百个节点 —— 实测对百万面模型，BVH 可将单次射线求交从 2–5ms 降到 0.03–0.1ms。
关键不是“建树快”，而是“遍历剔除率高”。如果三角形空间分布极不均匀（比如全部挤在角落），或者 BVH 分割策略太粗糙（如只按坐标轴中位数切，不考虑三角形朝向），剔除效果会断崖式下降。
用 embree 快速接入 BVH 射线求交
embree 是 Intel 维护的工业级 BVH 库，支持 AVX2/AVX-512 加速，C++ 接口干净，无需手写 BVH 构建逻辑。它不强制你管理内存布局，但要求顶点和索引满足特定对齐与顺序约束：
- 顶点数组必须是
float[3]连续排列（即std::vector<:array>></:array>或float* v = new float[3*n];） - 三角形索引必须是
int32_t[3]每组（不是uint16_t，哪怕面数 - 创建场景前需调用
rtcInit()；，销毁后调用rtcExit()，否则内存泄漏且后续调用可能崩溃 - 射线结构体
RTCIntersectContext必须栈分配（不能 new），且每次调用rtcIntersect1()前要重置context.flags = RTC_INTERSECT_CONTEXT_FLAG_COHERENT（对拾取这种局部相关射线更高效）
RTCDevice device = rtcNewDevice(nullptr); RTCScene scene = rtcNewScene(device); rtcSetSceneFlags(scene, RTC_SCENE_FLAG_COMPACT | RTC_SCENE_FLAG_ROBUST); <p>// 假设 verts 和 indices 已就绪 RTCGeometry geom = rtcNewGeometry(device, RTC_GEOMETRY_TYPE_TRIANGLE); rtcSetSharedGeometryBuffer(geom, RTC_BUFFER_TYPE_VERTEX, 0, RTC_FORMAT_FLOAT3, verts, 0, 3 <em> sizeof(float), nverts); rtcSetSharedGeometryBuffer(geom, RTC_BUFFER_TYPE_INDEX, 0, RTC_FORMAT_UINT32, indices, 0, 3 </em> sizeof(uint32_t), ntris); rtcCommitGeometry(geom); rtcAttachGeometry(scene, geom); rtcCommitScene(scene);</p>
rtcIntersect1() 返回结果怎么解读
rtcIntersect1() 不返回布尔值，而是直接覆写你传入的 RTCRayHit 结构体。常见误判点：
立即学习“C++免费学习笔记（深入）”；
-
ray.tfar被改写为最近交点距离 —— 若无交点，它保持原值（通常初始化为std::numeric_limits<float>::max()</float>），**不能靠hit.geomID == RTC_INVALID_GEOMETRY_ID判定失败**，因为有些情况下 ID 仍有效但 tfar 未更新 -
hit.primID是三角形索引（对应indices数组的第 i 组），不是顶点索引；若需获取顶点坐标，得用indices[hit.primID * 3 + 0]等去查verts - 拾取时经常忽略射线方向归一化：embree 内部假设
ray.dir是单位向量，否则tmin/tmax范围计算错误，导致近处交点被裁掉 - 若模型做了世界变换，别把矩阵应用到每个顶点再传给 embree —— 更高效的做法是用
rtcSetGeometryTransform()设置局部到世界变换，让 embree 在 BVH 遍历时自动处理
自己实现 BVH 构建时最容易崩的三个地方
不用 embree 而手写 BVH（比如为了动态更新或教学），90% 的 crash 和慢性能来自这三点：
- 构建递归深度失控：未设节点数或三角形数阈值，小网格也分到叶子只剩 1 个三角形，树高爆炸。建议叶子节点控制在 4–16 个三角形，用
std::vector存子节点而非裸指针，避免栈溢出 - AABB 合并错误：对三角形顶点算包围盒时，用了
min(x0,x1,x2)但忘了初始化为FLT_MAX，导致初始值是 0，包围盒严重偏小 - 射线-AABB 相交测试用错公式：网上很多代码抄的是 “slab method” 但漏了对
ray.dir.x == 0的分支处理，导致除零或 NaN 传播；正确做法是用if (abs(dir.x) 跳过该轴，或统一用 <code>inv_dir.x = 1.0f / dir.x预计算倒数（需提前检查零）
BVH 的真正难点不在算法本身，而在于数值鲁棒性 —— 尤其当模型含退化三角形（三点共线）、超大坐标（如地理坐标系下 x=1e7）、或射线起点紧贴表面时，浮点误差会直接让交点漂移甚至消失。调试时务必打开 RTC_SCENE_FLAG_ROBUST，并用 rtcDebug() + assert 检查每层包围盒是否真包含其子节点。
