Physics & Rendering Expert: Rope Dynamics & Constraint Solving
Expert in computational physics for real-time rope/cable dynamics, constraint solving, and physically-based simulations.
When to Use This Skill
Use for:
- Real-time rope/cable/chain simulation (leashes, climbing ropes)
- Position-Based Dynamics (PBD) implementation
- Constraint solvers (Gauss-Seidel, Jacobi)
- Quaternion/dual-quaternion rotation math
- Verlet integration for particle systems
- Tangle detection (multi-rope collisions)
Do NOT use for:
- Fluid dynamics → specialized SPH/MPM solvers
- Fracture simulation → requires FEM or MPM
- Offline cinematic physics → different constraints
- Unity/Unreal physics → use built-in systems
Expert vs Novice Shibboleths
| Topic |
Novice |
Expert |
| Constraint approach |
Uses spring forces (F=ma) |
Uses PBD (directly manipulates positions) |
| Why PBD |
"Springs work fine" |
Springs require tiny timesteps; PBD is unconditionally stable |
| Solver choice |
"Just iterate until done" |
Gauss-Seidel for chains, Jacobi for GPU |
| Iterations |
20+ iterations |
5-10 is optimal; diminishing returns after |
| Rotation |
Uses Euler angles |
Uses quaternions (no gimbal lock) |
| Integration |
Forward Euler |
Verlet (symplectic, energy-conserving) |
Common Anti-Patterns
Force-Based Springs for Stiff Constraints
| What it looks like |
Why it's wrong |
force = k * (distance - rest_length) with high k |
High k requires tiny dt for stability; low k gives squishy ropes |
| Instead: Use PBD - directly move particles to satisfy constraints |
|
Euler Angles for Rotation
| What it looks like |
Why it's wrong |
rotation = vec3(pitch, yaw, roll) |
Gimbal lock at 90° pitch; unstable composition |
| Instead: Use quaternions - 4 numbers, no gimbal lock, stable SLERP |
|
Over-Iteration
| What it looks like |
Why it's wrong |
solver_iterations = 50 |
Diminishing returns after 5-10; wastes cycles |
| Instead: Use 5-10 iterations; if more needed, use XPBD compliance |
|
Single-Threaded Gauss-Seidel for Large Systems
| What it looks like |
Why it's wrong |
| Gauss-Seidel on 1000+ constraints |
Gauss-Seidel is inherently sequential |
| Instead: Use Jacobi solver for GPU parallelization |
|
Quick Reference
Why PBD Beats Force-Based Physics
- Unconditionally stable (large timesteps OK)
- Direct control over constraint satisfaction
- No spring constants to tune
- Predictable behavior
Solver Choice
| Solver |
Parallelizable |
Convergence |
Use Case |
| Gauss-Seidel |
No |
Fast |
Chains, ropes |
| Jacobi |
Yes (GPU) |
Slower |
Large meshes, cloth |
Rotation Representation
- 3D rotation → Quaternion (never Euler)
- Rotation + translation → Dual quaternion
- Skinning/blending → Dual quaternion (no candy-wrapper artifact)
Performance Targets
| System |
Budget |
Notes |
| Single rope (100 particles) |
<0.5ms |
5 iterations sufficient |
| Three-dog leash (60 particles) |
<0.7ms |
Include tangle detection |
| Cloth (1000 particles) |
<2ms |
Use Jacobi on GPU |
Evolution Timeline
| Era |
Key Development |
| Pre-2006 |
Mass-spring systems, stability issues |
| 2006-2015 |
PBD introduced (Müller et al.), unconditional stability |
| 2016-2020 |
XPBD adds compliance for soft constraints |
| 2021-2024 |
ALEM (2024 SIGGRAPH), BDEM, neural physics |
| 2025+ |
XPBD standard, hybrid CPU/GPU, learned corrections |
Decision Trees
Choosing constraint solver:
- Sequential structure (rope/chain)? → Gauss-Seidel
- Large parallel system (cloth/hair)? → Jacobi (GPU)
- Need soft constraints? → XPBD with compliance
Choosing integration:
- Position-only needed? → Basic Verlet
- Need velocity for forces? → Velocity Verlet
- High accuracy required? → RK4 (but PBD usually sufficient)
Integrates With
- metal-shader-expert - GPU compute shaders for Jacobi solver
- native-app-designer - Visualization and debugging UI
Reference Files
| File |
Contents |
references/core-algorithms.md |
PBD loop, Verlet, quaternions, solver implementations |
references/tangle-physics.md |
Multi-rope collision, Capstan friction, TangleConstraint |
Remember: Real-time physics is about stability and visual plausibility, not physical accuracy. PBD with 5-10 iterations at 60fps looks great and runs fast.
1---2name: physics-rendering-expert-23description: Real-time rope/cable physics using Position-Based Dynamics (PBD), Verlet integration, and constraint solvers. Expert in quaternion math, Gauss-Seidel/Jacobi solvers, and tangling detection. Activate on 'rope simulation', 'PBD', 'Position-Based Dynamics', 'Verlet', 'constraint solver', 'quaternion', 'cable dynamics', 'cloth simulation', 'leash physics'. NOT for fluid dynamics (SPH/MPM), fracture simulation (FEM), offline cinematic physics, molecular dynamics, or general game physics engines (use Unity/Unreal built-ins).4---56# Physics & Rendering Expert: Rope Dynamics & Constraint Solving78Expert in computational physics for real-time rope/cable dynamics, constraint solving, and physically-based simulations.910## When to Use This Skill1112**Use for:**13- Real-time rope/cable/chain simulation (leashes, climbing ropes)14- Position-Based Dynamics (PBD) implementation15- Constraint solvers (Gauss-Seidel, Jacobi)16- Quaternion/dual-quaternion rotation math17- Verlet integration for particle systems18- Tangle detection (multi-rope collisions)1920**Do NOT use for:**21- Fluid dynamics → specialized SPH/MPM solvers22- Fracture simulation → requires FEM or MPM23- Offline cinematic physics → different constraints24- Unity/Unreal physics → use built-in systems2526## Expert vs Novice Shibboleths2728| Topic | Novice | Expert |29|-------|--------|--------|30| **Constraint approach** | Uses spring forces (F=ma) | Uses PBD (directly manipulates positions) |31| **Why PBD** | "Springs work fine" | Springs require tiny timesteps; PBD is unconditionally stable |32| **Solver choice** | "Just iterate until done" | Gauss-Seidel for chains, Jacobi for GPU |33| **Iterations** | 20+ iterations | 5-10 is optimal; diminishing returns after |34| **Rotation** | Uses Euler angles | Uses quaternions (no gimbal lock) |35| **Integration** | Forward Euler | Verlet (symplectic, energy-conserving) |3637## Common Anti-Patterns3839### Force-Based Springs for Stiff Constraints40| What it looks like | Why it's wrong |41|--------------------|----------------|42| `force = k * (distance - rest_length)` with high k | High k requires tiny dt for stability; low k gives squishy ropes |43| **Instead**: Use PBD - directly move particles to satisfy constraints |4445### Euler Angles for Rotation46| What it looks like | Why it's wrong |47|--------------------|----------------|48| `rotation = vec3(pitch, yaw, roll)` | Gimbal lock at 90° pitch; unstable composition |49| **Instead**: Use quaternions - 4 numbers, no gimbal lock, stable SLERP |5051### Over-Iteration52| What it looks like | Why it's wrong |53|--------------------|----------------|54| `solver_iterations = 50` | Diminishing returns after 5-10; wastes cycles |55| **Instead**: Use 5-10 iterations; if more needed, use XPBD compliance |5657### Single-Threaded Gauss-Seidel for Large Systems58| What it looks like | Why it's wrong |59|--------------------|----------------|60| Gauss-Seidel on 1000+ constraints | Gauss-Seidel is inherently sequential |61| **Instead**: Use Jacobi solver for GPU parallelization |6263## Quick Reference6465### Why PBD Beats Force-Based Physics6667- Unconditionally stable (large timesteps OK)68- Direct control over constraint satisfaction69- No spring constants to tune70- Predictable behavior7172### Solver Choice7374| Solver | Parallelizable | Convergence | Use Case |75|--------|---------------|-------------|----------|76| **Gauss-Seidel** | No | Fast | Chains, ropes |77| **Jacobi** | Yes (GPU) | Slower | Large meshes, cloth |7879### Rotation Representation8081- 3D rotation → Quaternion (never Euler)82- Rotation + translation → Dual quaternion83- Skinning/blending → Dual quaternion (no candy-wrapper artifact)8485### Performance Targets8687| System | Budget | Notes |88|--------|--------|-------|89| Single rope (100 particles) | <0.5ms | 5 iterations sufficient |90| Three-dog leash (60 particles) | <0.7ms | Include tangle detection |91| Cloth (1000 particles) | <2ms | Use Jacobi on GPU |9293## Evolution Timeline9495| Era | Key Development |96|-----|-----------------|97| Pre-2006 | Mass-spring systems, stability issues |98| 2006-2015 | PBD introduced (Müller et al.), unconditional stability |99| 2016-2020 | XPBD adds compliance for soft constraints |100| 2021-2024 | ALEM (2024 SIGGRAPH), BDEM, neural physics |101| 2025+ | XPBD standard, hybrid CPU/GPU, learned corrections |102103## Decision Trees104105**Choosing constraint solver:**106- Sequential structure (rope/chain)? → Gauss-Seidel107- Large parallel system (cloth/hair)? → Jacobi (GPU)108- Need soft constraints? → XPBD with compliance109110**Choosing integration:**111- Position-only needed? → Basic Verlet112- Need velocity for forces? → Velocity Verlet113- High accuracy required? → RK4 (but PBD usually sufficient)114115## Integrates With116117- **metal-shader-expert** - GPU compute shaders for Jacobi solver118- **native-app-designer** - Visualization and debugging UI119120## Reference Files121122| File | Contents |123|------|----------|124| `references/core-algorithms.md` | PBD loop, Verlet, quaternions, solver implementations |125| `references/tangle-physics.md` | Multi-rope collision, Capstan friction, TangleConstraint |126127---128129**Remember**: Real-time physics is about stability and visual plausibility, not physical accuracy. PBD with 5-10 iterations at 60fps looks great and runs fast.
curiositech