Use Huilier's theorem to compute the spherical excess in the solidangle() method.

The former trigonometric method led to different results on different platforms due to numerical imprecisions. We now use Huilier's theorem to be less sensitive to numerical problems. The polygon is decomposed in two triangles for this purpose.
Adapt reference values to new GWcs::solidangle method using Huilier's theorem.