An Efficient Multi-solution Solver for the Inverse Kinematics of 3-Section Constant-Curvature Robots