diff --git a/README.rst b/README.rst
index 0997e9bc97ff11be477955163c44f5099358c543..0940247b080c8f37832b70ba875d6b9422119a22 100644
--- a/README.rst
+++ b/README.rst
@@ -1,5 +1,7 @@
-.. image:: https://gitlab.mpcdf.mpg.de/mjr/bfps/badges/develop/build.svg
-   :target: https://gitlab.mpcdf.mpg.de/mjr/bfps/commits/develop
+.. only:: html
+
+    .. image:: https://gitlab.mpcdf.mpg.de/mjr/bfps/badges/develop/build.svg
+        :target: https://gitlab.mpcdf.mpg.de/mjr/bfps/commits/develop
 
 =========================================
 Turbulence Tools: Lagrangian and Eulerian
diff --git a/cpp/kspace.cpp b/cpp/kspace.cpp
index 213fbcbb48e229c4217a1c61d109121ae47e843c..d759219872b7fd9efb1bf6ca98909aa766881649 100644
--- a/cpp/kspace.cpp
+++ b/cpp/kspace.cpp
@@ -352,25 +352,33 @@ void kspace<be, dt>::Gauss_filter(
  *  a Gaussian filter and a sharp real space spherical filter.
  *
  *  Filter expressions in real space are as follows:
- *  \f{eqnarray*}{
- *       \phi^b_\ell(r) &=&
- *           \frac{1}{\ell^3}\frac{6}{\pi} H(\ell/2 - r) \\
- *       \phi^g_\ell(r) &=&
+ *  \f[
+ *       \phi^b_\ell(r) =
+ *          \frac{1}{\ell^3}\frac{6}{\pi} H(\ell/2 - r)
+ *  \f]
+ *  \f[
+ *       \phi^g_\ell(r) =
  *           \frac{1}{\sigma_\ell^3}\frac{1}{(2\pi)^{3/2}}
- *           \exp\left(-\frac{1}{2}\left(\frac{r}{\sigma_\ell}\right)^2\right) \\
- *       \phi^s_\ell(r) &=&
+ *           \exp\left(-\frac{1}{2}\left(\frac{r}{\sigma_\ell}\right)^2\right)
+ *  \f]
+ *  \f[
+ *       \phi^s_\ell(r) =
  *           \frac{1}{2 \pi^2 r^3}
  *           \left(\sin k_\ell r - k_\ell r \cos k_\ell r\right)
- *  \f}
+ *  \f]
  *  and the corresponding expressions in Fourier space are:
- *  \f{eqnarray*}{
- *       \hat{\phi^b_\ell}(k) &=&
+ *  \f[
+ *       \hat{\phi^b_\ell}(k) =
  *       \frac{3}{2(k\ell/2)^3}
- *       \left(2\sin (k \ell/2) - k \ell \cos (k \ell/2)\right) \\
- *       \hat{\phi^g_\ell}(k) &=&
- *       \exp\left(-\frac{1}{2}k^2 \sigma_\ell^2\right) \\
- *       \hat{\phi^s_\ell}(k) &=& H(k_\ell - k)
- *  \f}
+ *       \left(2\sin (k \ell/2) - k \ell \cos (k \ell/2)\right)
+ *  \f]
+ *  \f[
+ *       \hat{\phi^g_\ell}(k) =
+ *       \exp\left(-\frac{1}{2}k^2 \sigma_\ell^2\right)
+ *  \f]
+ *  \f[
+ *       \hat{\phi^s_\ell}(k) = H(k_\ell - k)
+ *  \f]
  *
  *  \f$ k_\ell \f$ is given as a parameter, and then we use
  *  \f[
@@ -425,7 +433,7 @@ int kspace<be, dt>::filter(
  *  a Gaussian filter and a sharp real space spherical filter.
  *
  *  Filter expressions in real space are as follows:
- *  \f{eqnarray*}{
+ *  \begin{eqnarray*}
  *      \phi^b_\ell(r) &=&
  *          \frac{1}{\ell^3}\frac{6}{\pi} H(\ell/2 - r) \\
  *      \phi^g_\ell(r) &=&
@@ -434,16 +442,16 @@ int kspace<be, dt>::filter(
  *      \phi^s_\ell(r) &=&
  *          \frac{1}{2 \pi^2 r^3}
  *          \left(\sin k_\ell r - k_\ell r \cos k_\ell r\right)
- *  \f}
+ *  \end{eqnarray*}
  *  and the corresponding expressions in Fourier space are:
- *  \f{eqnarray*}{
+ *  \begin{eqnarray*}
  *      \hat{\phi^b_\ell}(k) &=&
  *      \frac{3}{2(k\ell/2)^3}
  *      \left(2\sin (k \ell/2) - k \ell \cos (k \ell/2)\right) \\
  *      \hat{\phi^g_\ell}(k) &=&
  *      \exp\left(-\frac{1}{2}k^2 \sigma_\ell^2\right) \\
  *      \hat{\phi^s_\ell}(k) &=& H(k_\ell - k)
- *  \f}
+ *  \end{eqnarray*}
  *
  *  \f$\sigma_\ell\f$ and \f$k_\ell\f$ are calibrated such that the energy of
  *  the large scales is approximately the same (within the inertial range)
diff --git a/documentation/conf.py b/documentation/conf.py
index e37f8e523e2039897c993372686bc0c15ef0ccdf..be67c310772ab154210d7bb5c348b0e594d11e8b 100644
--- a/documentation/conf.py
+++ b/documentation/conf.py
@@ -225,7 +225,7 @@ latex_elements = {
 #'pointsize': '10pt',
 
 # Additional stuff for the LaTeX preamble.
-#'preamble': '',
+'preamble': '\usepackage{amssymb,amsmath,amsfonts}',
 
 # Latex figure (float) alignment
 #'figure_align': 'htbp',