From 3cf78785a9b0d3e87788b304dcac81fe7c98698f Mon Sep 17 00:00:00 2001
From: Martin Reinecke <martin@mpa-garching.mpg.de>
Date: Wed, 24 May 2017 07:34:47 +0200
Subject: [PATCH] convert z to colatitudes

---
 pyHealpix.cc | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/pyHealpix.cc b/pyHealpix.cc
index 18673e5..a59c039 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;
   }
 
-- 
GitLab