Newton SolverVBD CUDA 加速升级刚-软接触与可变形弹性求解性能优化解析【免费下载链接】newtonAn open-source, GPU-accelerated physics simulation engine built upon NVIDIA Warp, specifically targeting roboticists and simulation researchers.项目地址: https://gitcode.com/GitHub_Trending/newton9/newton本篇技术指南聚焦 Newton基于 NVIDIA Warp 的 GPU 加速物理仿真引擎在 changelog 条目 4141 中记录的一项核心性能变更SolverVBD在 CUDA 上对刚体-软体rigid-soft接触与可变形弹性deformable elasticity求解的加速以及伴随而来的软接触候选对按 shape 排序的底层优化。读完本文你将掌握该变更的适用场景、底层实现原理涉及 solver_vbd.py、collide.py 等源码、性能特征边界以及它可能带来的浮点结果差异与应对方式。变更概述一次面向 CUDA 大规模接触场景的求解加速changelog 片段 4141.changed.md 记录了本次变更的核心事实Speed upSolverVBDrigid-soft contact and deformable elasticity on CUDA, most in scenes with many contacts or worlds. Exact values can shift within floating-point tolerance relative to earlier releases.同组的另一片段 4141.changed.1.md 补充了配套优化Sort soft-contact candidate pairs by shape on every device. Results can shift within floating-point tolerance, and when contact capacity overflows the retained subset can differ. Faster on scenes with many shapes or worlds; a little slower on small dense-contact scenes.概括而言本次变更包含两个层面求解器内核加速SolverVBD的刚-软接触求解与可变形弹性求解在 CUDA 设备上被提速在接触数量多或 world并行世界数量多的场景收益最明显。候选对排序优化软接触候选对在所有设备上都改为按 shape 稳定排序让连续候选对共享同一 shape 的变换/缩放/SDF 数据与类型分发分支提升内存访问与分支效率。两者共同的代价是浮点结果可能相对旧版本发生微小偏移在接触容量溢出的极端情况下保留的接触子集可能不同——这些都属于可接受的数值容忍范围内变化。SolverVBD 是什么VBD AVBD 的统一求解器要理解这次加速改了什么首先需要明确SolverVBD的定位。在 solver_vbd.py 的类文档中SolverVBD被定义为粒子布料、软体、MPM 柔体采用 Vertex Block DescentVBD算法刚体关节、接触采用 Augmented Vertex Block DescentAVBD算法二者通过同一求解器统一处理并支持粒子-刚体耦合integrate_with_external_rigid_solver可切换为单向耦合模式。该求解器支持粒子自接触、刚体-刚体接触、以及刚体-粒子刚-软接触三大类交互其公开构造参数中与本次加速直接相关的包括参数默认值与本次变更的关系particle_enable_tile_solveTrue是否使用 tile API 加速粒子弹性求解仅 CUDA 生效对应“可变形弹性在 CUDA 上加速”rigid_body_particle_contact_buffer_size256每个刚体可跟踪的刚-软接触列表容量粒子 全表面 edge/facerigid_contact_historyFalse是否跨步持久化接触数值状态warm startiterations10每步 VBD 迭代次数deterministicNone是否启用原子发射内核的确定性模式继承wp.config.deterministic从源码结构看本次“刚-软接触”加速主要落在SolverVBD内部的 body-particle 接触路径上涉及的内核包括 rigid_vbd_kernels.py 中的accumulate_body_particle_contacts_per_body、build_body_particle_contact_lists、update_duals_body_particle_contacts等以及 particle_vbd_kernels.py 中的gather_particle_body_contact_force_and_hessian、build_particle_body_contact_adjacency_active等“可变形弹性”加速则对应粒子侧solve_elasticity/ tile 版弹性求解内核。加速原理一软接触候选对按 shape 稳定排序4141 变更的另一半在 collide.py 的_world_compatible_pairs函数中得到印证。该函数负责生成软体特征particle/edge/face与 shape 之间的候选对其注释明确描述了排序的动机collide.pyPairs are stably sorted by shape index so consecutive candidates process the same shape: on CUDA a warp then reads one shapes transform/scale/SDF data and takes one type-dispatch branch.也就是说候选对在构造阶段host 端通过np.argsort(s_idx, kindstable)按 shape 索引稳定排序collide.py排序后CUDA 上一个 warp 内连续的线程会处理同一个 shape的接触候选因此可以复用同一份 shape 变换、缩放与 SDF 数据减少重复读取只走一次类型分发分支减少分支分歧提升指令效率每个接触记录会保存其候选 tidthread id因此下游映射不依赖候选顺序——这正是结果可以安全排序、且排序不会破坏接触映射正确性的关键设计。性能特征边界根据 changelog 的表述该排序优化具有明确的适用边界更快shape 多或 world 多的场景如大规模多世界仿真、含大量软体对象的场景略慢小规模高密度接触场景shape 数量少、接触候选密集时排序开销占比上升。因此在评估性能回归时应优先在多 shape、多 world 的典型生产场景中测量而不是只跑单 shape 的小场景。加速原理二CUDA 上的 tile 弹性求解与软接触路径SolverVBD在 CUDA 上的弹性加速依赖particle_enable_tile_solve默认True参数。其实现要点见 solver_vbd.py构造时根据模型静态数据四面体/三角形/边刚度特化生成tile 弹性内核make_solve_elasticity_tile(include_triangles, include_tets, two_particles_per_warp)。纯四面体模型会编译出不包含三角形/边代码路径的专用内核特化依赖元素材料因此修改三角形或边刚度后需要重建求解器否则特化不会更新__init__注释明确提示tile 求解仅 CUDA 生效在 CPU 设备上会自动禁用源码在 debug 级别输出 Tiled solve requires model.devicecuda.刚-软接触侧的加速还包括粒子-刚体接触邻接表的构建build_particle_body_contact_adjacency_active、每刚体接触列表的累积accumulate_body_particle_contacts_per_body等内核它们共同支撑“刚性 body 与软性 particle/edge/face 全表面接触”的高效求解。使用前提与注意事项SolverVBD需要模型调用ModelBuilder.color()或set_coloring()为粒子和刚体生成颜色分组否则构造时会抛出ValueError刚-软接触缓冲rigid_body_particle_contact_buffer_size等是固定预分配的运行期不会动态扩容设置过小会导致接触溢出源码中通过check_contact_overflow等机制检测设置过大会增加显存占用若开启rigid_contact_history进行接触数值 warm start需要CollisionPipeline使用contact_matchinglatest或sticky且在图捕获graph capture场景中应按CollisionPipeline→SolverVBD的顺序构造。数值变化与兼容性浮点容差内结果可偏移本次变更最重要的用户可见影响是结果数值的微小变化浮点容差偏移由于接触内核的计算顺序、归约方式原子累加、tile 分块发生变化接触力、弹性形变等数值可能相对旧版本在浮点容差内偏移接触容量溢出的子集差异当接触数超过预分配容量时保留的接触子集可能不同排序改变了截断的选择顺序这属于文档明确声明的可接受行为确定性SolverVBD支持通过deterministic参数或继承wp.config.deterministic对原子发射内核启用确定性模式需要逐 bit 可复现结果的场景应显式配置并注意确定性模式下可能需要更大的记录缓冲deterministic_max_records。验证与回归建议对依赖旧版精确数值的测试建议以物理量级与行为一致性如接触力误差在容差内、稳定站立/堆叠行为不变而非逐 bit 相等作为回归标准对比新旧性能时构造“多 shape、多 world、多接触”的压力场景同时记录小规模密集接触场景作为对照以覆盖两条性能边界。相关源码与测试定位若需深入验证或继续调研本次变更可在仓库中定位以下文件solver_vbd.pySolverVBD求解器主体含全部参数解析、tile 弹性内核特化与刚-软接触状态初始化rigid_vbd_kernels.py刚体 AVBD 与 body-particle 软接触内核接触列表构建、力/Hessian 累积、双变量更新particle_vbd_kernels.py粒子 VBD 弹性求解、自接触与粒子-刚体接触力/Hessian 内核collide.py_world_compatible_pairs实现软接触候选对按 shape 稳定排序np.argsort(..., kindstable)vbd_coupling_kernels.py粒子-刚体耦合接触力收割内核用于代理刚体proxy body耦合场景test_solver_vbd.py求解器行为与正确性的测试覆盖如软体-刚体接触、弹性形变相关用例CHANGELOG.md发行说明中 4141 条目的最终渲染位置由 Towncrier 在发布分支合并生成。小结SolverVBD的本次 CUDA 加速changelog 4141通过两条技术路线实现一是软接触候选对按 shape 稳定排序提升 CUDA warp 级数据复用与分支效率二是粒子弹性 tile 求解等内核层面的加速。它显著利好多 shape、多 world、多接触的规模化场景同时以浮点容差内的结果偏移和“小密集场景略慢”为代价。对于正在使用 NewtonSolverVBD做大规模软体-刚体耦合仿真的开发者建议在升级后以物理行为一致性为回归基准并在多接触压力场景中重新测量性能。【免费下载链接】newtonAn open-source, GPU-accelerated physics simulation engine built upon NVIDIA Warp, specifically targeting roboticists and simulation researchers.项目地址: https://gitcode.com/GitHub_Trending/newton9/newton创作声明：本文部分内容由AI辅助生成（AIGC），仅供参考 
长河君 · 2026/9/18 12:12:23 · 阅读 1.2k Newton SolverVBD CUDA 加速升级：刚-软接触与可变形弹性求解性能优化解析
Newton SolverVBD CUDA 加速升级刚-软接触与可变形弹性求解性能优化解析【免费下载链接】newtonAn open-source, GPU-accelerated physics simulation engine built upon NVIDIA Warp, specifically targeting roboticists and simulation researchers.项目地址: https://gitcode.com/GitHub_Trending/newton9/newton本篇技术指南聚焦 Newton基于 NVIDIA Warp 的 GPU 加速物理仿真引擎在 changelog 条目 4141 中记录的一项核心性能变更SolverVBD在 CUDA 上对刚体-软体rigid-soft接触与可变形弹性deformable elasticity求解的加速以及伴随而来的软接触候选对按 shape 排序的底层优化。读完本文你将掌握该变更的适用场景、底层实现原理涉及 solver_vbd.py、collide.py 等源码、性能特征边界以及它可能带来的浮点结果差异与应对方式。变更概述一次面向 CUDA 大规模接触场景的求解加速changelog 片段 4141.changed.md 记录了本次变更的核心事实Speed upSolverVBDrigid-soft contact and deformable elasticity on CUDA, most in scenes with many contacts or worlds. Exact values can shift within floating-point tolerance relative to earlier releases.同组的另一片段 4141.changed.1.md 补充了配套优化Sort soft-contact candidate pairs by shape on every device. Results can shift within floating-point tolerance, and when contact capacity overflows the retained subset can differ. Faster on scenes with many shapes or worlds; a little slower on small dense-contact scenes.概括而言本次变更包含两个层面求解器内核加速SolverVBD的刚-软接触求解与可变形弹性求解在 CUDA 设备上被提速在接触数量多或 world并行世界数量多的场景收益最明显。候选对排序优化软接触候选对在所有设备上都改为按 shape 稳定排序让连续候选对共享同一 shape 的变换/缩放/SDF 数据与类型分发分支提升内存访问与分支效率。两者共同的代价是浮点结果可能相对旧版本发生微小偏移在接触容量溢出的极端情况下保留的接触子集可能不同——这些都属于可接受的数值容忍范围内变化。SolverVBD 是什么VBD AVBD 的统一求解器要理解这次加速改了什么首先需要明确SolverVBD的定位。在 solver_vbd.py 的类文档中SolverVBD被定义为粒子布料、软体、MPM 柔体采用 Vertex Block DescentVBD算法刚体关节、接触采用 Augmented Vertex Block DescentAVBD算法二者通过同一求解器统一处理并支持粒子-刚体耦合integrate_with_external_rigid_solver可切换为单向耦合模式。该求解器支持粒子自接触、刚体-刚体接触、以及刚体-粒子刚-软接触三大类交互其公开构造参数中与本次加速直接相关的包括参数默认值与本次变更的关系particle_enable_tile_solveTrue是否使用 tile API 加速粒子弹性求解仅 CUDA 生效对应“可变形弹性在 CUDA 上加速”rigid_body_particle_contact_buffer_size256每个刚体可跟踪的刚-软接触列表容量粒子 全表面 edge/facerigid_contact_historyFalse是否跨步持久化接触数值状态warm startiterations10每步 VBD 迭代次数deterministicNone是否启用原子发射内核的确定性模式继承wp.config.deterministic从源码结构看本次“刚-软接触”加速主要落在SolverVBD内部的 body-particle 接触路径上涉及的内核包括 rigid_vbd_kernels.py 中的accumulate_body_particle_contacts_per_body、build_body_particle_contact_lists、update_duals_body_particle_contacts等以及 particle_vbd_kernels.py 中的gather_particle_body_contact_force_and_hessian、build_particle_body_contact_adjacency_active等“可变形弹性”加速则对应粒子侧solve_elasticity/ tile 版弹性求解内核。加速原理一软接触候选对按 shape 稳定排序4141 变更的另一半在 collide.py 的_world_compatible_pairs函数中得到印证。该函数负责生成软体特征particle/edge/face与 shape 之间的候选对其注释明确描述了排序的动机collide.pyPairs are stably sorted by shape index so consecutive candidates process the same shape: on CUDA a warp then reads one shapes transform/scale/SDF data and takes one type-dispatch branch.也就是说候选对在构造阶段host 端通过np.argsort(s_idx, kindstable)按 shape 索引稳定排序collide.py排序后CUDA 上一个 warp 内连续的线程会处理同一个 shape的接触候选因此可以复用同一份 shape 变换、缩放与 SDF 数据减少重复读取只走一次类型分发分支减少分支分歧提升指令效率每个接触记录会保存其候选 tidthread id因此下游映射不依赖候选顺序——这正是结果可以安全排序、且排序不会破坏接触映射正确性的关键设计。性能特征边界根据 changelog 的表述该排序优化具有明确的适用边界更快shape 多或 world 多的场景如大规模多世界仿真、含大量软体对象的场景略慢小规模高密度接触场景shape 数量少、接触候选密集时排序开销占比上升。因此在评估性能回归时应优先在多 shape、多 world 的典型生产场景中测量而不是只跑单 shape 的小场景。加速原理二CUDA 上的 tile 弹性求解与软接触路径SolverVBD在 CUDA 上的弹性加速依赖particle_enable_tile_solve默认True参数。其实现要点见 solver_vbd.py构造时根据模型静态数据四面体/三角形/边刚度特化生成tile 弹性内核make_solve_elasticity_tile(include_triangles, include_tets, two_particles_per_warp)。纯四面体模型会编译出不包含三角形/边代码路径的专用内核特化依赖元素材料因此修改三角形或边刚度后需要重建求解器否则特化不会更新__init__注释明确提示tile 求解仅 CUDA 生效在 CPU 设备上会自动禁用源码在 debug 级别输出 Tiled solve requires model.devicecuda.刚-软接触侧的加速还包括粒子-刚体接触邻接表的构建build_particle_body_contact_adjacency_active、每刚体接触列表的累积accumulate_body_particle_contacts_per_body等内核它们共同支撑“刚性 body 与软性 particle/edge/face 全表面接触”的高效求解。使用前提与注意事项SolverVBD需要模型调用ModelBuilder.color()或set_coloring()为粒子和刚体生成颜色分组否则构造时会抛出ValueError刚-软接触缓冲rigid_body_particle_contact_buffer_size等是固定预分配的运行期不会动态扩容设置过小会导致接触溢出源码中通过check_contact_overflow等机制检测设置过大会增加显存占用若开启rigid_contact_history进行接触数值 warm start需要CollisionPipeline使用contact_matchinglatest或sticky且在图捕获graph capture场景中应按CollisionPipeline→SolverVBD的顺序构造。数值变化与兼容性浮点容差内结果可偏移本次变更最重要的用户可见影响是结果数值的微小变化浮点容差偏移由于接触内核的计算顺序、归约方式原子累加、tile 分块发生变化接触力、弹性形变等数值可能相对旧版本在浮点容差内偏移接触容量溢出的子集差异当接触数超过预分配容量时保留的接触子集可能不同排序改变了截断的选择顺序这属于文档明确声明的可接受行为确定性SolverVBD支持通过deterministic参数或继承wp.config.deterministic对原子发射内核启用确定性模式需要逐 bit 可复现结果的场景应显式配置并注意确定性模式下可能需要更大的记录缓冲deterministic_max_records。验证与回归建议对依赖旧版精确数值的测试建议以物理量级与行为一致性如接触力误差在容差内、稳定站立/堆叠行为不变而非逐 bit 相等作为回归标准对比新旧性能时构造“多 shape、多 world、多接触”的压力场景同时记录小规模密集接触场景作为对照以覆盖两条性能边界。相关源码与测试定位若需深入验证或继续调研本次变更可在仓库中定位以下文件solver_vbd.pySolverVBD求解器主体含全部参数解析、tile 弹性内核特化与刚-软接触状态初始化rigid_vbd_kernels.py刚体 AVBD 与 body-particle 软接触内核接触列表构建、力/Hessian 累积、双变量更新particle_vbd_kernels.py粒子 VBD 弹性求解、自接触与粒子-刚体接触力/Hessian 内核collide.py_world_compatible_pairs实现软接触候选对按 shape 稳定排序np.argsort(..., kindstable)vbd_coupling_kernels.py粒子-刚体耦合接触力收割内核用于代理刚体proxy body耦合场景test_solver_vbd.py求解器行为与正确性的测试覆盖如软体-刚体接触、弹性形变相关用例CHANGELOG.md发行说明中 4141 条目的最终渲染位置由 Towncrier 在发布分支合并生成。小结SolverVBD的本次 CUDA 加速changelog 4141通过两条技术路线实现一是软接触候选对按 shape 稳定排序提升 CUDA warp 级数据复用与分支效率二是粒子弹性 tile 求解等内核层面的加速。它显著利好多 shape、多 world、多接触的规模化场景同时以浮点容差内的结果偏移和“小密集场景略慢”为代价。对于正在使用 NewtonSolverVBD做大规模软体-刚体耦合仿真的开发者建议在升级后以物理行为一致性为回归基准并在多接触压力场景中重新测量性能。【免费下载链接】newtonAn open-source, GPU-accelerated physics simulation engine built upon NVIDIA Warp, specifically targeting roboticists and simulation researchers.项目地址: https://gitcode.com/GitHub_Trending/newton9/newton创作声明：本文部分内容由AI辅助生成（AIGC），仅供参考