A parallel physics solver may solve an equation to determine the characteristics of motion after a collision. In some embodiments, the physics solver solves the equation AX=B, where A is a sparse constrained matrix and B is the right hand vector. The sparse constrained matrix may be formed of 6×K blocks, where K is a tuning parameter that divides into the width of a single instruction multiple data processor used to implement the physics solver, without residue.