FA-87159 / Physics integrator stability / Member archive

Quaternion derivative omits the one-half factor · case 04

Bodies spin at twice the commanded angular velocity.

Member previewVariant 4 · 3 implementations · 8 checks per implementation

Case contract

solve(q, omega, dt, steps): orientation quaternion [w,x,y,z], world-frame angular velocity omega. Each step q += 0.5*dt*(0,omega) (x) q (Hamilton product, omega on the left), then renormalize to unit length. Return q with w>=0 (negate all components otherwise), rounded to 6.

Why this case matters

Game and robotics physics loops depend on integrator update order, step control and stabilization terms; a wrong decision point turns a stable simulation into drifting or exploding motion.

One recorded failure

Sample boundary fixture

This sample comes from the broken implementation of a controlled reproducer.

Boundary fixture	Actual	Expected	Outcome
case 0	`[0.383734, 0.504892, 0.64242, -0.430267]`	`[0.378645, 0.510325, 0.556833, -0.534914]`	Failed

MEMBER ARCHIVE

The complete case is available to members.

This record includes three runnable implementations, regression fixtures, execution results, and source hashes.

Member access is invitation-based. Sign in with your invited account to inspect the sources.