Working rotation optimisation (angle-axis, not quaternion).
Good options are num_threads=4, robustify