diff --git a/pyHealpix.cc b/pyHealpix.cc
index 18673e54c3a9b2651033e049ac0bb65dd40efd81..a59c039e8a35cdcb8a4f4047c7e4ca7e86d3c47a 100644
--- a/pyHealpix.cc
+++ b/pyHealpix.cc
@@ -491,6 +491,8 @@ a_d_c GL_thetas(int64_t nlat)
   auto rr=res.mutable_unchecked<1>();
   vector<double> dummy_weights(nlat);
   sharp_legendre_roots(nlat, &rr[0], dummy_weights.data());
+  for (size_t i=0; i<rr.shape(0); ++i)
+    rr[i]=acos(-rr[i]);
   return res;
   }