Core: review computation of pole figure
If I understood & remember correctly, there is a change of algorithm at alpha=10deg in order to avoid a spurious singularity at the pole alpha=0. Perhaps one can do better, and generate the entire pole figure with one single algorithm.