Commit 1d1b856e authored by Stefan Mauerberger's avatar Stefan Mauerberger
Browse files

equi_sph() got a new interface

parent 14b3423c
......@@ -5,7 +5,7 @@
proxys of these) for quantities of interest, such as the power spectrum,
the dipole intensity, the dipole location etc.
Copyright (C) 2019 Maximilian Schanner, GFZ Potsdam
Copyright (C) 2019 Maximilian Schanner, GFZ Potsdam
Copyright (C) 2019 Stefan Mauerberger, University of Potsdam
Cite as:
......@@ -149,7 +149,8 @@ def location(posterior, mu_dips, cov_dips, n_points=10000,
The dipole location pdf, evaluated at the gridpoints
"""
# set up grid to evaluate on
theta, phi = zip(*equi_sph(n_points, bounds=bounds))
(ti, tf), (pi, pf) = bounds
theta, phi = zip(*equi_sph(n_points, ti, tf, pi, pf))
phi = np.asarray(phi)
theta = np.pi/2. - np.asarray(theta)
# allocate array for the density
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment