diff --git a/docs/pyerrors/input/bdio.html b/docs/pyerrors/input/bdio.html
index 3873932a..b9efe50f 100644
--- a/docs/pyerrors/input/bdio.html
+++ b/docs/pyerrors/input/bdio.html
@@ -466,7 +466,7 @@
 </span><span id="L-384"><a href="#L-384"><span class="linenos">384</span></a>            <span class="n">d_buf</span> <span class="o">=</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_double</span> <span class="o">*</span> <span class="p">(</span><span class="mi">2</span> <span class="o">+</span> <span class="n">d0</span> <span class="o">*</span> <span class="n">d1</span> <span class="o">*</span> <span class="mi">2</span><span class="p">)</span>
 </span><span id="L-385"><a href="#L-385"><span class="linenos">385</span></a>            <span class="n">pd_buf</span> <span class="o">=</span> <span class="n">d_buf</span><span class="p">()</span>
 </span><span id="L-386"><a href="#L-386"><span class="linenos">386</span></a>            <span class="n">ppd_buf</span> <span class="o">=</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_void_p</span><span class="p">(</span><span class="n">ctypes</span><span class="o">.</span><span class="n">addressof</span><span class="p">(</span><span class="n">pd_buf</span><span class="p">))</span>
-</span><span id="L-387"><a href="#L-387"><span class="linenos">387</span></a>            <span class="n">iread</span> <span class="o">=</span> <span class="n">bdio_read_f64</span><span class="p">(</span><span class="n">ppd_buf</span><span class="p">,</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_size_t</span><span class="p">(</span><span class="n">rlen</span><span class="p">),</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_void_p</span><span class="p">(</span><span class="n">fbdio</span><span class="p">))</span>
+</span><span id="L-387"><a href="#L-387"><span class="linenos">387</span></a>            <span class="n">bdio_read_f64</span><span class="p">(</span><span class="n">ppd_buf</span><span class="p">,</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_size_t</span><span class="p">(</span><span class="n">rlen</span><span class="p">),</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_void_p</span><span class="p">(</span><span class="n">fbdio</span><span class="p">))</span>
 </span><span id="L-388"><a href="#L-388"><span class="linenos">388</span></a>            <span class="k">if</span> <span class="n">corr_type</span><span class="p">[</span><span class="n">corr_no</span><span class="p">]</span> <span class="o">==</span> <span class="s1">&#39;complex&#39;</span><span class="p">:</span>
 </span><span id="L-389"><a href="#L-389"><span class="linenos">389</span></a>                <span class="n">tmp_mean</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">mean</span><span class="p">(</span><span class="n">np</span><span class="o">.</span><span class="n">asarray</span><span class="p">(</span><span class="n">np</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="n">np</span><span class="o">.</span><span class="n">asarray</span><span class="p">(</span><span class="n">pd_buf</span><span class="p">[</span><span class="mi">2</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">d1</span><span class="p">:</span><span class="o">-</span><span class="mi">2</span> <span class="o">*</span> <span class="n">d1</span><span class="p">:</span><span class="mi">2</span><span class="p">]),</span> <span class="n">d0</span> <span class="o">-</span> <span class="mi">2</span><span class="p">)),</span> <span class="n">axis</span><span class="o">=</span><span class="mi">1</span><span class="p">)</span>
 </span><span id="L-390"><a href="#L-390"><span class="linenos">390</span></a>            <span class="k">else</span><span class="p">:</span>
@@ -675,7 +675,7 @@
 </span><span id="L-593"><a href="#L-593"><span class="linenos">593</span></a>            <span class="n">d_buf</span> <span class="o">=</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_double</span> <span class="o">*</span> <span class="p">(</span><span class="mi">2</span> <span class="o">+</span> <span class="n">d0</span><span class="p">)</span>
 </span><span id="L-594"><a href="#L-594"><span class="linenos">594</span></a>            <span class="n">pd_buf</span> <span class="o">=</span> <span class="n">d_buf</span><span class="p">()</span>
 </span><span id="L-595"><a href="#L-595"><span class="linenos">595</span></a>            <span class="n">ppd_buf</span> <span class="o">=</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_void_p</span><span class="p">(</span><span class="n">ctypes</span><span class="o">.</span><span class="n">addressof</span><span class="p">(</span><span class="n">pd_buf</span><span class="p">))</span>
-</span><span id="L-596"><a href="#L-596"><span class="linenos">596</span></a>            <span class="n">iread</span> <span class="o">=</span> <span class="n">bdio_read_f64</span><span class="p">(</span><span class="n">ppd_buf</span><span class="p">,</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_size_t</span><span class="p">(</span><span class="n">rlen</span><span class="p">),</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_void_p</span><span class="p">(</span><span class="n">fbdio</span><span class="p">))</span>
+</span><span id="L-596"><a href="#L-596"><span class="linenos">596</span></a>            <span class="n">bdio_read_f64</span><span class="p">(</span><span class="n">ppd_buf</span><span class="p">,</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_size_t</span><span class="p">(</span><span class="n">rlen</span><span class="p">),</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_void_p</span><span class="p">(</span><span class="n">fbdio</span><span class="p">))</span>
 </span><span id="L-597"><a href="#L-597"><span class="linenos">597</span></a>            <span class="n">tmp_mean</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">mean</span><span class="p">(</span><span class="n">np</span><span class="o">.</span><span class="n">asarray</span><span class="p">(</span><span class="n">pd_buf</span><span class="p">[</span><span class="mi">2</span><span class="p">:]))</span>
 </span><span id="L-598"><a href="#L-598"><span class="linenos">598</span></a>
 </span><span id="L-599"><a href="#L-599"><span class="linenos">599</span></a>            <span class="n">data</span><span class="p">[</span><span class="n">corr_no</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">tmp_mean</span><span class="p">)</span>
@@ -1208,7 +1208,7 @@ adding the flag -fPIC to CC and changing the all target to</p>
 </span><span id="read_mesons-385"><a href="#read_mesons-385"><span class="linenos">385</span></a>            <span class="n">d_buf</span> <span class="o">=</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_double</span> <span class="o">*</span> <span class="p">(</span><span class="mi">2</span> <span class="o">+</span> <span class="n">d0</span> <span class="o">*</span> <span class="n">d1</span> <span class="o">*</span> <span class="mi">2</span><span class="p">)</span>
 </span><span id="read_mesons-386"><a href="#read_mesons-386"><span class="linenos">386</span></a>            <span class="n">pd_buf</span> <span class="o">=</span> <span class="n">d_buf</span><span class="p">()</span>
 </span><span id="read_mesons-387"><a href="#read_mesons-387"><span class="linenos">387</span></a>            <span class="n">ppd_buf</span> <span class="o">=</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_void_p</span><span class="p">(</span><span class="n">ctypes</span><span class="o">.</span><span class="n">addressof</span><span class="p">(</span><span class="n">pd_buf</span><span class="p">))</span>
-</span><span id="read_mesons-388"><a href="#read_mesons-388"><span class="linenos">388</span></a>            <span class="n">iread</span> <span class="o">=</span> <span class="n">bdio_read_f64</span><span class="p">(</span><span class="n">ppd_buf</span><span class="p">,</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_size_t</span><span class="p">(</span><span class="n">rlen</span><span class="p">),</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_void_p</span><span class="p">(</span><span class="n">fbdio</span><span class="p">))</span>
+</span><span id="read_mesons-388"><a href="#read_mesons-388"><span class="linenos">388</span></a>            <span class="n">bdio_read_f64</span><span class="p">(</span><span class="n">ppd_buf</span><span class="p">,</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_size_t</span><span class="p">(</span><span class="n">rlen</span><span class="p">),</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_void_p</span><span class="p">(</span><span class="n">fbdio</span><span class="p">))</span>
 </span><span id="read_mesons-389"><a href="#read_mesons-389"><span class="linenos">389</span></a>            <span class="k">if</span> <span class="n">corr_type</span><span class="p">[</span><span class="n">corr_no</span><span class="p">]</span> <span class="o">==</span> <span class="s1">&#39;complex&#39;</span><span class="p">:</span>
 </span><span id="read_mesons-390"><a href="#read_mesons-390"><span class="linenos">390</span></a>                <span class="n">tmp_mean</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">mean</span><span class="p">(</span><span class="n">np</span><span class="o">.</span><span class="n">asarray</span><span class="p">(</span><span class="n">np</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="n">np</span><span class="o">.</span><span class="n">asarray</span><span class="p">(</span><span class="n">pd_buf</span><span class="p">[</span><span class="mi">2</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">d1</span><span class="p">:</span><span class="o">-</span><span class="mi">2</span> <span class="o">*</span> <span class="n">d1</span><span class="p">:</span><span class="mi">2</span><span class="p">]),</span> <span class="n">d0</span> <span class="o">-</span> <span class="mi">2</span><span class="p">)),</span> <span class="n">axis</span><span class="o">=</span><span class="mi">1</span><span class="p">)</span>
 </span><span id="read_mesons-391"><a href="#read_mesons-391"><span class="linenos">391</span></a>            <span class="k">else</span><span class="p">:</span>
@@ -1460,7 +1460,7 @@ Manually overwrite ensemble name</li>
 </span><span id="read_dSdm-594"><a href="#read_dSdm-594"><span class="linenos">594</span></a>            <span class="n">d_buf</span> <span class="o">=</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_double</span> <span class="o">*</span> <span class="p">(</span><span class="mi">2</span> <span class="o">+</span> <span class="n">d0</span><span class="p">)</span>
 </span><span id="read_dSdm-595"><a href="#read_dSdm-595"><span class="linenos">595</span></a>            <span class="n">pd_buf</span> <span class="o">=</span> <span class="n">d_buf</span><span class="p">()</span>
 </span><span id="read_dSdm-596"><a href="#read_dSdm-596"><span class="linenos">596</span></a>            <span class="n">ppd_buf</span> <span class="o">=</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_void_p</span><span class="p">(</span><span class="n">ctypes</span><span class="o">.</span><span class="n">addressof</span><span class="p">(</span><span class="n">pd_buf</span><span class="p">))</span>
-</span><span id="read_dSdm-597"><a href="#read_dSdm-597"><span class="linenos">597</span></a>            <span class="n">iread</span> <span class="o">=</span> <span class="n">bdio_read_f64</span><span class="p">(</span><span class="n">ppd_buf</span><span class="p">,</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_size_t</span><span class="p">(</span><span class="n">rlen</span><span class="p">),</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_void_p</span><span class="p">(</span><span class="n">fbdio</span><span class="p">))</span>
+</span><span id="read_dSdm-597"><a href="#read_dSdm-597"><span class="linenos">597</span></a>            <span class="n">bdio_read_f64</span><span class="p">(</span><span class="n">ppd_buf</span><span class="p">,</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_size_t</span><span class="p">(</span><span class="n">rlen</span><span class="p">),</span> <span class="n">ctypes</span><span class="o">.</span><span class="n">c_void_p</span><span class="p">(</span><span class="n">fbdio</span><span class="p">))</span>
 </span><span id="read_dSdm-598"><a href="#read_dSdm-598"><span class="linenos">598</span></a>            <span class="n">tmp_mean</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">mean</span><span class="p">(</span><span class="n">np</span><span class="o">.</span><span class="n">asarray</span><span class="p">(</span><span class="n">pd_buf</span><span class="p">[</span><span class="mi">2</span><span class="p">:]))</span>
 </span><span id="read_dSdm-599"><a href="#read_dSdm-599"><span class="linenos">599</span></a>
 </span><span id="read_dSdm-600"><a href="#read_dSdm-600"><span class="linenos">600</span></a>            <span class="n">data</span><span class="p">[</span><span class="n">corr_no</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">tmp_mean</span><span class="p">)</span>