The current code for computing spherical barycenters in ConvexArea2S.getBarycenter() seems to suffer from floating point accuracy issues. The ConvexArea2STest.checkBarycenterConsistency() method checks the consistency of the barycenter computation of a region by splitting the region into two sections, computing the area and barycenter of each section, and then computing the combined barycenter of the sections by adding the barycenter of each scaled by its corresponding area. It is expected that the combined barycenter computed in this way should equal the barycenter computed for the region as a whole. However, in practice, a large epsilon value is needed in the comparison in order for the tests to pass. We need to investigate why this is the case.