Do not use underscores for folders

* This may fixes the build warning
* Also this may happen in the future:
readthedocs/readthedocs.org#1776 (comment)

Author: wangenau

docs/conf.py

@@ -24,11 +24,11 @@
     'sphinx_rtd_theme'
 ]
 exclude_patterns = []
-templates_path = ['_templates']
+templates_path = ['templates']
 pygments_style = 'sphinx'
 
 
 # -- Options for HTML output -------------------------------------------------
 
 html_theme = 'sphinx_rtd_theme'
-html_static_path = ['_static']
+html_static_path = ['static']

docs/_images/calculation_time.png renamed to docs/images/calculation_time.png

@@ -3,6 +3,7 @@
 Introduction
 ************
 
+
 What is this?
 =============
 
@@ -28,33 +29,33 @@ Minimization algorithm
 | The default setting for the variational mesh function is strictly bound to PySCF's grid levels (0-9) that are just lists of ten radial and angular grids for every group of elements.
 | To show how the algorithm works, a molecule with two atomic species will be used (:math:`\mathrm{H}_2\mathrm{O}`). Assuming that more grid points mean a better representation of the density (for an atomic species) and therefore results in a smaller error, we use this knowledge to generate more accurate grids. Our optimal grid for a specific threshold would be somewhere in the combination space shown in the figure below.
 
-.. figure:: _images/comb0.png
+.. figure:: images/comb0.png
    :align: center
    :figwidth: 50%
 
 First, we will do a coarse grid search (seen the figure below, colored in red) where the grid levels for every atomic species are the same. Eventually, we find a match (colored in blue), i.e., a combination whose error is below our threshold. If there is no match the best possible grid will be used.
 
-.. figure:: _images/comb1.png
+.. figure:: images/comb1.png
    :align: center
    :figwidth: 50%
 
 | With the found match we can eliminate further grid level pairs, shown below. It is safe to remove the upper right regime, since raising either one of the grid levels will only result in a larger grid than the found one.
 | We can also remove the bottom left part. With no match in the previous steps of the coarse grid search, it is unlikely to find a smaller grid below the threshold while reducing grid levels.
 
-.. figure:: _images/comb2.png
+.. figure:: images/comb2.png
    :align: center
    :figwidth: 50%
 
 At last, we can remove every combination that would have a larger number of grid points than our match since we are only interested in the minimum. For :math:`\mathrm{H}_2\mathrm{O}` this step can be seen, again, in the figure below. Afterward, we go through the list of possible candidates sorted by growing numbers of grid points and calculate the resulting error until we eventually get a new match (colored in green). This will be referred to as the fine grid search.
 
-.. figure:: _images/comb3.png
+.. figure:: images/comb3.png
    :align: center
    :figwidth: 50%
 
 Program flow
 ------------
 
-.. figure:: _images/uml_diagram.png
+.. figure:: images/uml_diagram.png
    :align: center
    :figwidth: 50%
 

docs/_images/comb0.png renamed to docs/images/comb0.png

@@ -14,4 +14,3 @@ Visualization
 =============
 .. automodule:: var_mesh.helpers
    :members:
-

docs/_images/comb1.png renamed to docs/images/comb1.png

@@ -12,6 +12,7 @@ Tutorial
 
    >>> import var_mesh
 
+
 Minimal example
 ===============
 
@@ -117,7 +118,7 @@ The package comes with functions to visualize meshes. The function :func:`plot_m
    >>> mesh.build()
    >>> plot_mesh_3d(mesh=mesh, weight=True)
 
-.. figure:: _images/plot_mesh_3d.png
+.. figure:: images/plot_mesh_3d.png
    :align: center
    :figwidth: 75%
 
@@ -128,7 +129,7 @@ The grid can be projected to a given plane, too. Also, the grid points can be sc
    >>> from var_mesh import plot_mesh_2d
    >>> plot_mesh_2d(mesh=mesh, weight=5, plane='xz')
 
-.. figure:: _images/plot_mesh_2d.png
+.. figure:: images/plot_mesh_2d.png
    :align: center
    :figwidth: 75%
 
@@ -167,6 +168,7 @@ This can be compared to the output, when the parameter ``precise`` is set to ``T
 
 The script for this example can be downloaded :download:`here </../examples/04_precise_option/04_precise_option.py>`.
 
+
 Mesh time
 =========
 
@@ -235,7 +237,7 @@ These result can be plotted afterwards.
    >>> plt.legend()
    >>> plt.show()
 
-.. figure:: _images/mesh_time.png
+.. figure:: images/mesh_time.png
    :align: center
    :figwidth: 75%
 
@@ -303,7 +305,7 @@ These result can be plotted as well.
   >>> plt.legend()
   >>> plt.show()
 
-.. figure:: _images/calculation_time.png
+.. figure:: images/calculation_time.png
    :align: center
    :figwidth: 75%