include GL in readme

Author: CosmoMatt

README.md

@@ -57,8 +57,9 @@ isolattitude sampling scheme. A number of sampling schemes are currently
 supported.
 
 The equiangular sampling schemes of [McEwen & Wiaux
-(2012)](https://arxiv.org/abs/1110.6298) and [Driscoll & Healy
-(1995)](https://www.sciencedirect.com/science/article/pii/S0196885884710086)
+(2012)](https://arxiv.org/abs/1110.6298), [Driscoll & Healy
+(1995)](https://www.sciencedirect.com/science/article/pii/S0196885884710086) 
+and [Gauss-Legendre (1986)](https://link.springer.com/article/10.1007/BF02519350)
 are supported, which exhibit associated sampling theorems and so
 harmonic transforms can be computed to machine precision. Note that the
 McEwen & Wiaux sampling theorem reduces the Nyquist rate on the sphere

docs/index.rst

@@ -21,7 +21,8 @@ Algorithms |:zap:|
 ``S2FFT`` leverages new algorithmic structures that can he highly parallelised and
 distributed, and so map very well onto the architecture of hardware accelerators (i.e.
 GPUs and TPUs).  In particular, these algorithms are based on new Wigner-d recursions
-that are stable to high angular resolution :math:`L`.  The diagram below illustrates the recursions (for further details see Price & McEwen 2023).
+that are stable to high angular resolution :math:`L`.  The diagram below illustrates the 
+recursions (for further details see Price & McEwen 2023).
 
 .. image:: ./assets/figures/Wigner_recursion_github_docs.png
 
@@ -46,7 +47,8 @@ Sampling |:earth_africa:|
 
 The structure of the algorithms implemented in ``S2FFT`` can support any isolattitude sampling scheme.  A number of sampling schemes are currently supported.
 
-The equiangular sampling schemes of `McEwen & Wiaux (2012) <https://arxiv.org/abs/1110.6298>`_ and `Driscoll & Healy (1995) <https://www.sciencedirect.com/science/article/pii/S0196885884710086>`_ are supported, which exhibit associated sampling theorems and so harmonic transforms can be computed to machine precision.  Note that the McEwen & Wiaux sampling theorem reduces the Nyquist rate on the sphere by a factor of two compared to the Driscoll & Healy approach, halving the number of spherical samples required. 
+The equiangular sampling schemes of `McEwen & Wiaux (2012) <https://arxiv.org/abs/1110.6298>`_, 
+`Driscoll & Healy (1995) <https://www.sciencedirect.com/science/article/pii/S0196885884710086>`_, and `Gauss-Legendre (1986) <https://link.springer.com/article/10.1007/BF02519350>`_ are supported, which exhibit associated sampling theorems and so harmonic transforms can be computed to machine precision.  Note that the McEwen & Wiaux sampling theorem reduces the Nyquist rate on the sphere by a factor of two compared to the Driscoll & Healy approach, halving the number of spherical samples required. 
 
 The popular `HEALPix <https://healpix.jpl.nasa.gov>`_ sampling scheme (`Gorski et al. 2005 <https://arxiv.org/abs/astro-ph/0409513>`_) is also supported.  The HEALPix sampling does not exhibit a sampling theorem and so the corresponding harmonic transforms do not achieve machine precision but exhibit some error.  However, the HEALPix sampling provides pixels of equal areas, which has many practical advantages.
 

requirements/requirements-core.txt

@@ -3,4 +3,7 @@ numpy>=1.20
 colorlog
 pyyaml
 jax>=0.3.13
-jaxlib
+jaxlib
+
+# Remove when subpackage functionality is fixed.
+torch