diff --git a/docs/index.rst b/docs/index.rst
index e999b810..216607ff 100644
--- a/docs/index.rst
+++ b/docs/index.rst
@@ -23,6 +23,8 @@ Contents
 .. toctree::
 
     source/expression_tree/index
+    source/geometry/index
+    source/meshes/index
     source/models/index
     source/parameters/index
     source/discretisations/index
diff --git a/docs/source/geometry/geometry.rst b/docs/source/geometry/geometry.rst
new file mode 100644
index 00000000..d99e3948
--- /dev/null
+++ b/docs/source/geometry/geometry.rst
@@ -0,0 +1,14 @@
+Geometry
+==========
+
+.. autoclass:: pybamm.Geometry
+    :members:
+
+.. autoclass:: pybamm.Geometry1DMacro
+    :members:
+
+.. autoclass:: pybamm.Geometry1DMicro
+    :members:
+    
+.. autoclass:: pybamm.Geometry3DMacro
+    :members:
diff --git a/docs/source/geometry/index.rst b/docs/source/geometry/index.rst
new file mode 100644
index 00000000..d2dc988d
--- /dev/null
+++ b/docs/source/geometry/index.rst
@@ -0,0 +1,7 @@
+Geometry
+========
+
+.. toctree::
+
+  geometry
+ 
diff --git a/docs/source/meshes/index.rst b/docs/source/meshes/index.rst
new file mode 100644
index 00000000..7a518d02
--- /dev/null
+++ b/docs/source/meshes/index.rst
@@ -0,0 +1,7 @@
+Meshes
+======
+
+.. toctree::
+  
+  meshes
+  submeshes
diff --git a/docs/source/meshes/meshes.rst b/docs/source/meshes/meshes.rst
new file mode 100644
index 00000000..f6788bb7
--- /dev/null
+++ b/docs/source/meshes/meshes.rst
@@ -0,0 +1,5 @@
+Meshes
+==========
+
+.. autoclass:: pybamm.Mesh
+    :members:
diff --git a/docs/source/meshes/submeshes.rst b/docs/source/meshes/submeshes.rst
new file mode 100644
index 00000000..ea9d6b44
--- /dev/null
+++ b/docs/source/meshes/submeshes.rst
@@ -0,0 +1,8 @@
+Sub Meshes
+==========
+
+.. autoclass:: pybamm.SubMesh1D
+    :members:
+
+.. autoclass:: pybamm.Uniform1DSubMesh
+    :members:
diff --git a/pybamm/__init__.py b/pybamm/__init__.py
index 0a27774a..7bf5b0a3 100644
--- a/pybamm/__init__.py
+++ b/pybamm/__init__.py
@@ -125,8 +125,9 @@ from .discretisations.finite_volume_discretisations import (
     FiniteVolumeDiscretisation,
     NodeToEdge,
 )
-from .discretisations.meshes import KNOWN_DOMAINS
-from .discretisations.meshes import Mesh, SubMesh1D, Uniform1DSubMesh
+from .meshes.meshes import KNOWN_DOMAINS
+from .meshes.meshes import Mesh
+from .meshes.submeshes import SubMesh1D, Uniform1DSubMesh
 
 #
 # Simulation class
diff --git a/pybamm/discretisations/base_discretisation.py b/pybamm/discretisations/base_discretisation.py
index e5326af3..067cfdb7 100644
--- a/pybamm/discretisations/base_discretisation.py
+++ b/pybamm/discretisations/base_discretisation.py
@@ -23,25 +23,14 @@ class BaseDiscretisation(object):
         the types of submeshes being employed on each of the subdomains
     """
 
-    def __init__(self, mesh_type, submesh_pts, submesh_types):
-        self.mesh_type = mesh_type
-        self.submesh_pts = submesh_pts
-        self.submesh_types = submesh_types
-
-        self._mesh = {}
+    def __init__(self, mesh):
+        self._mesh = mesh
 
     @property
     def mesh(self):
         return self._mesh
 
-    @mesh.setter
-    def mesh(self, mesh_in):
-        self._mesh = mesh_in
-
-    def mesh_geometry(self, geometry):
-        self._mesh = self.mesh_type(geometry, self.submesh_types, self.submesh_pts)
-
-    def process_model(self, model, geometry):
+    def process_model(self, model):
         """Discretise a model.
         Currently inplace, could be changed to return a new model.
 
@@ -53,9 +42,6 @@ class BaseDiscretisation(object):
 
         """
 
-        # Set the mesh
-        self._mesh = self.mesh_type(geometry, self.submesh_types, self.submesh_pts)
-
         # Set the y split for variables
         variables = list(model.rhs.keys()) + list(model.algebraic.keys())
         y_slices = self.get_variable_slices(variables)
diff --git a/pybamm/discretisations/finite_volume_discretisations.py b/pybamm/discretisations/finite_volume_discretisations.py
index 86ac9133..7f50883d 100644
--- a/pybamm/discretisations/finite_volume_discretisations.py
+++ b/pybamm/discretisations/finite_volume_discretisations.py
@@ -20,8 +20,8 @@ class FiniteVolumeDiscretisation(pybamm.BaseDiscretisation):
     **Extends:** :class:`pybamm.BaseDiscretisation`
     """
 
-    def __init__(self, mesh_type, submesh_pts, submesh_types):
-        super().__init__(mesh_type, submesh_pts, submesh_types)
+    def __init__(self, mesh):
+        super().__init__(mesh)
 
     def compute_diffusivity(self, symbol):
         """
diff --git a/pybamm/discretisations/meshes.py b/pybamm/meshes/meshes.py
similarity index 69%
rename from pybamm/discretisations/meshes.py
rename to pybamm/meshes/meshes.py
index cbc267d9..f31ab1fe 100644
--- a/pybamm/discretisations/meshes.py
+++ b/pybamm/meshes/meshes.py
@@ -74,7 +74,7 @@ class Mesh(dict):
                 [self[submeshnames[0]].edges]
                 + [self[submeshname].edges[1:] for submeshname in submeshnames[1:]]
             )
-            return SubMesh1D(combined_submesh_edges)
+            return pybamm.SubMesh1D(combined_submesh_edges)
         else:
             raise pybamm.DomainError("submesh edges are not aligned")
 
@@ -93,47 +93,8 @@ class Mesh(dict):
 
             # left ghost cell: two edges, one node, to the left of existing submesh
             lgs_edges = np.array([2 * edges[0] - edges[1], edges[0]])
-            self[submeshname + "_left ghost cell"] = SubMesh1D(lgs_edges)
+            self[submeshname + "_left ghost cell"] = pybamm.SubMesh1D(lgs_edges)
 
             # right ghost cell: two edges, one node, to the right of existing submesh
             rgs_edges = np.array([edges[-1], 2 * edges[-1] - edges[-2]])
-            self[submeshname + "_right ghost cell"] = SubMesh1D(rgs_edges)
-
-
-class SubMesh1D:
-    """ Submesh class.
-        Contains the position of the nodes and the number of mesh points.
-
-        Parameters
-        ----------
-        domain : dict
-            A dictionary that contains the limits of the spatial variables
-        npts : dict
-            A dictionary that contains the number of points to be used on each
-            spatial variable
-        """
-
-    def __init__(self, edges):
-        self.edges = edges
-        self.nodes = (self.edges[1:] + self.edges[:-1]) / 2
-        self.d_edges = np.diff(self.edges)
-        self.d_nodes = np.diff(self.nodes)
-        self.npts = self.nodes.size
-
-
-class Uniform1DSubMesh(SubMesh1D):
-    """
-    A class to generate a uniform submesh on a 1D domain
-    """
-
-    def __init__(self, lims, npts):
-
-        # currently accept lims and npts as dicts. This may get changed at a future
-        # date depending on the form of mesh we desire for 2D/3D
-
-        spatial_lims = list(lims.values())[0]
-        npts = list(npts.values())[0]
-
-        edges = np.linspace(spatial_lims["min"], spatial_lims["max"], npts + 1)
-
-        super().__init__(edges)
+            self[submeshname + "_right ghost cell"] = pybamm.SubMesh1D(rgs_edges)
diff --git a/pybamm/meshes/submeshes.py b/pybamm/meshes/submeshes.py
new file mode 100644
index 00000000..9914b305
--- /dev/null
+++ b/pybamm/meshes/submeshes.py
@@ -0,0 +1,46 @@
+#
+# Native PyBaMM Meshes
+#
+from __future__ import absolute_import, division
+from __future__ import print_function, unicode_literals
+
+import numpy as np
+
+
+class SubMesh1D:
+    """ Submesh class.
+        Contains the position of the nodes and the number of mesh points.
+
+        Parameters
+        ----------
+        domain : dict
+            A dictionary that contains the limits of the spatial variables
+        npts : dict
+            A dictionary that contains the number of points to be used on each
+            spatial variable
+        """
+
+    def __init__(self, edges):
+        self.edges = edges
+        self.nodes = (self.edges[1:] + self.edges[:-1]) / 2
+        self.d_edges = np.diff(self.edges)
+        self.d_nodes = np.diff(self.nodes)
+        self.npts = self.nodes.size
+
+
+class Uniform1DSubMesh(SubMesh1D):
+    """
+    A class to generate a uniform submesh on a 1D domain
+    """
+
+    def __init__(self, lims, npts):
+
+        # currently accept lims and npts as dicts. This may get changed at a future
+        # date depending on the form of mesh we desire for 2D/3D
+
+        spatial_lims = list(lims.values())[0]
+        npts = list(npts.values())[0]
+
+        edges = np.linspace(spatial_lims["min"], spatial_lims["max"], npts + 1)
+
+        super().__init__(edges)
diff --git a/pybamm/models/base_model.py b/pybamm/models/base_model.py
index d79435a2..78184343 100644
--- a/pybamm/models/base_model.py
+++ b/pybamm/models/base_model.py
@@ -67,11 +67,8 @@ class BaseModel(object):
             "positive electrode": pybamm.Uniform1DSubMesh,
         }
 
-        mesh_type = pybamm.Mesh
-
-        self.default_discretisation = pybamm.FiniteVolumeDiscretisation(
-            mesh_type, submesh_pts, submesh_types
-        )
+        self.mesh = pybamm.Mesh(self.default_geometry, submesh_types, submesh_pts)
+        self.default_discretisation = pybamm.FiniteVolumeDiscretisation(self.mesh)
         self.default_solver = pybamm.ScipySolver(method="RK45")
 
     def _set_dict(self, dict, name):

Refactor Meshing to Occur Outside Discretisation
Take mesh creation out of the discretisation class. The mesh should be created by taking in a geometry, submesh_types, and submesh_pts. 

The discretisation should now take in a mesh either in the init function or when we process model. 



diff --git a/tests/shared.py b/tests/shared.py
index 2ab35152..20fc12ec 100644
--- a/tests/shared.py
+++ b/tests/shared.py
@@ -7,7 +7,7 @@ import numpy as np
 
 
 class TestDefaults1DMacro:
-    def __init__(self):
+    def __init__(self, npts=0):
         self.param = pybamm.ParameterValues(
             base_parameters={"Ln": 0.3, "Ls": 0.3, "Lp": 0.3}
         )
@@ -20,21 +20,22 @@ class TestDefaults1DMacro:
             "separator": {"x": 25},
             "positive electrode": {"x": 35},
         }
+
         self.submesh_types = {
             "negative electrode": pybamm.Uniform1DSubMesh,
             "separator": pybamm.Uniform1DSubMesh,
             "positive electrode": pybamm.Uniform1DSubMesh,
         }
 
-        self.mesh_type = pybamm.Mesh
+        if npts != 0:
+            n = 3 * round(npts / 3)
+            self.submesh_pts = {
+                "negative electrode": {"x": n},
+                "separator": {"x": n},
+                "positive electrode": {"x": n},
+            }
 
-    def set_equal_pts(self, npts):
-        n = 3 * round(npts / 3)
-        self.submesh_pts = {
-            "negative electrode": {"x": n},
-            "separator": {"x": n},
-            "positive electrode": {"x": n},
-        }
+        self.mesh = pybamm.Mesh(self.geometry, self.submesh_types, self.submesh_pts)
 
 
 class TestDefaults1DParticle:
@@ -46,19 +47,17 @@ class TestDefaults1DParticle:
         }
         self.param = pybamm.ParameterValues(base_parameters={})
         self.param.process_geometry(self.geometry)
-        self.mesh_type = pybamm.Mesh
         self.submesh_pts = {"negative particle": {"r": n}}
         self.submesh_types = {"negative particle": pybamm.Uniform1DSubMesh}
 
+        self.mesh = pybamm.Mesh(self.geometry, self.submesh_types, self.submesh_pts)
+
 
 class DiscretisationForTesting(pybamm.BaseDiscretisation):
     """Identity operators, no boundary conditions."""
 
-    def __init__(self, mesh_type, submesh_pts, submesh_types):
-        super().__init__(mesh_type, submesh_pts, submesh_types)
-
-    def set_mesh(self, geometry):
-        self._mesh = self.mesh_type(geometry, self.submesh_types, self.submesh_pts)
+    def __init__(self, mesh):
+        super().__init__(mesh)
 
     def gradient(self, symbol, y_slices, boundary_conditions):
         discretised_symbol = self.process_symbol(symbol, y_slices, boundary_conditions)
diff --git a/tests/test_discretisations/test_base_discretisations.py b/tests/test_discretisations/test_base_discretisations.py
index e23d5bc3..8f0379df 100644
--- a/tests/test_discretisations/test_base_discretisations.py
+++ b/tests/test_discretisations/test_base_discretisations.py
@@ -22,10 +22,7 @@ class TestDiscretise(unittest.TestCase):
 
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = shared.DiscretisationForTesting(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = shared.DiscretisationForTesting(defaults.mesh)
 
         result = disc._concatenate_init(initial_conditions, y_slices)
 
@@ -41,10 +38,7 @@ class TestDiscretise(unittest.TestCase):
     def test_discretise_slicing(self):
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = shared.DiscretisationForTesting(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = shared.DiscretisationForTesting(defaults.mesh)
         mesh = disc.mesh
 
         whole_cell = ["negative electrode", "separator", "positive electrode"]
@@ -78,10 +72,7 @@ class TestDiscretise(unittest.TestCase):
     def test_process_symbol_base(self):
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = shared.DiscretisationForTesting(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = shared.DiscretisationForTesting(defaults.mesh)
 
         # variable
         var = pybamm.Variable("var")
@@ -140,10 +131,7 @@ class TestDiscretise(unittest.TestCase):
 
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = shared.DiscretisationForTesting(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = shared.DiscretisationForTesting(defaults.mesh)
 
         y_slices = {var1.id: slice(53), var2.id: slice(53, 59)}
         exp_disc = disc.process_symbol(expression, y_slices)
@@ -180,10 +168,7 @@ class TestDiscretise(unittest.TestCase):
     def test_discretise_spatial_operator(self):
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = shared.DiscretisationForTesting(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = shared.DiscretisationForTesting(defaults.mesh)
         mesh = disc.mesh
 
         whole_cell = ["negative electrode", "separator", "positive electrode"]
@@ -225,7 +210,7 @@ class TestDiscretise(unittest.TestCase):
 
     def test_core_NotImplementedErrors(self):
         # create discretisation
-        disc = pybamm.BaseDiscretisation(None, None, None)
+        disc = pybamm.BaseDiscretisation(None)
 
         with self.assertRaises(NotImplementedError):
             disc.gradient(None, None, {})
@@ -244,10 +229,7 @@ class TestDiscretise(unittest.TestCase):
 
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = shared.DiscretisationForTesting(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = shared.DiscretisationForTesting(defaults.mesh)
         mesh = disc.mesh
 
         combined_submesh = mesh.combine_submeshes(*whole_cell)
@@ -310,13 +292,11 @@ class TestDiscretise(unittest.TestCase):
 
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = shared.DiscretisationForTesting(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.process_model(model, defaults.geometry)
+        disc = shared.DiscretisationForTesting(defaults.mesh)
         mesh = disc.mesh
 
         combined_submesh = mesh.combine_submeshes(*whole_cell)
+        disc.process_model(model)
 
         y0 = model.concatenated_initial_conditions
         np.testing.assert_array_equal(y0, 3 * np.ones_like(combined_submesh.nodes))
@@ -344,7 +324,7 @@ class TestDiscretise(unittest.TestCase):
         }
         model.variables = {"ST": S * T}
 
-        disc.process_model(model, defaults.geometry)
+        disc.process_model(model)
         y0 = model.concatenated_initial_conditions
         np.testing.assert_array_equal(
             y0,
@@ -370,7 +350,7 @@ class TestDiscretise(unittest.TestCase):
         model.boundary_conditions = {}
         model.variables = {"ST": S * T}
         with self.assertRaises(pybamm.ModelError):
-            disc.process_model(model, defaults.geometry)
+            disc.process_model(model)
 
     def test_process_model_dae(self):
         # one rhs equation and one algebraic
@@ -391,12 +371,10 @@ class TestDiscretise(unittest.TestCase):
 
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = shared.DiscretisationForTesting(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.process_model(model, defaults.geometry)
+        disc = shared.DiscretisationForTesting(defaults.mesh)
         mesh = disc.mesh
 
+        disc.process_model(model)
         combined_submesh = mesh.combine_submeshes(*whole_cell)
 
         y0 = model.concatenated_initial_conditions
@@ -441,17 +419,14 @@ class TestDiscretise(unittest.TestCase):
         model.variables = {"c": c, "N": N, "d": d}
 
         with self.assertRaises(pybamm.ModelError):
-            disc.process_model(model, defaults.geometry)
+            disc.process_model(model)
 
     def test_broadcast(self):
         whole_cell = ["negative electrode", "separator", "positive electrode"]
 
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = shared.DiscretisationForTesting(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = shared.DiscretisationForTesting(defaults.mesh)
         mesh = disc.mesh
 
         combined_submesh = mesh.combine_submeshes(*whole_cell)
@@ -495,10 +470,7 @@ class TestDiscretise(unittest.TestCase):
 
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = shared.DiscretisationForTesting(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = shared.DiscretisationForTesting(defaults.mesh)
 
         conc = disc.concatenate(a, b, c)
         self.assertIsInstance(conc, pybamm.NumpyModelConcatenation)
@@ -510,10 +482,7 @@ class TestDiscretise(unittest.TestCase):
 
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = shared.DiscretisationForTesting(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = shared.DiscretisationForTesting(defaults.mesh)
         mesh = disc.mesh
 
         var = pybamm.Variable("var", domain=whole_cell)
@@ -533,10 +502,7 @@ class TestDiscretise(unittest.TestCase):
     def test_discretise_space(self):
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = shared.DiscretisationForTesting(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = shared.DiscretisationForTesting(defaults.mesh)
 
         # space
         x1 = pybamm.Space(["negative electrode"])
diff --git a/tests/test_discretisations/test_finite_volume_discretisations.py b/tests/test_discretisations/test_finite_volume_discretisations.py
index 37164dff..4e9cade2 100644
--- a/tests/test_discretisations/test_finite_volume_discretisations.py
+++ b/tests/test_discretisations/test_finite_volume_discretisations.py
@@ -38,10 +38,7 @@ class TestFiniteVolumeDiscretisation(unittest.TestCase):
 
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = pybamm.FiniteVolumeDiscretisation(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = pybamm.FiniteVolumeDiscretisation(defaults.mesh)
         mesh = disc.mesh
 
         combined_submesh = mesh.combine_submeshes(*whole_cell)
@@ -96,10 +93,7 @@ class TestFiniteVolumeDiscretisation(unittest.TestCase):
 
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = pybamm.FiniteVolumeDiscretisation(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = pybamm.FiniteVolumeDiscretisation(defaults.mesh)
         mesh = disc.mesh
 
         # Add ghost nodes
@@ -142,10 +136,7 @@ class TestFiniteVolumeDiscretisation(unittest.TestCase):
         whole_cell = ["negative electrode", "separator", "positive electrode"]
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = pybamm.FiniteVolumeDiscretisation(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = pybamm.FiniteVolumeDiscretisation(defaults.mesh)
         mesh = disc.mesh
 
         combined_submesh = mesh.combine_submeshes(*whole_cell)
@@ -192,10 +183,7 @@ class TestFiniteVolumeDiscretisation(unittest.TestCase):
         """
         # create discretisation
         defaults = shared.TestDefaults1DParticle(10)
-        disc = pybamm.FiniteVolumeDiscretisation(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = pybamm.FiniteVolumeDiscretisation(defaults.mesh)
         mesh = disc.mesh
 
         combined_submesh = mesh.combine_submeshes("negative particle")
@@ -248,10 +236,7 @@ class TestFiniteVolumeDiscretisation(unittest.TestCase):
 
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = pybamm.FiniteVolumeDiscretisation(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = pybamm.FiniteVolumeDiscretisation(defaults.mesh)
         mesh = disc.mesh
 
         combined_submesh = mesh.combine_submeshes(*whole_cell)
@@ -291,10 +276,7 @@ class TestFiniteVolumeDiscretisation(unittest.TestCase):
 
         # create discretisation
         defaults = shared.TestDefaults1DParticle(10)
-        disc = pybamm.FiniteVolumeDiscretisation(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = pybamm.FiniteVolumeDiscretisation(defaults.mesh)
         mesh = disc.mesh
 
         combined_submesh = mesh.combine_submeshes("negative particle")
@@ -338,10 +320,7 @@ class TestFiniteVolumeDiscretisation(unittest.TestCase):
         """
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = pybamm.FiniteVolumeDiscretisation(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = pybamm.FiniteVolumeDiscretisation(defaults.mesh)
         mesh = disc.mesh
 
         mesh.add_ghost_meshes()
@@ -396,12 +375,8 @@ class TestFiniteVolumeDiscretisation(unittest.TestCase):
             # Set up discretisation
             n = 3 * round(n / 3)
             # create discretisation
-            defaults = shared.TestDefaults1DMacro()
-            defaults.set_equal_pts(n)
-            disc = pybamm.FiniteVolumeDiscretisation(
-                defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-            )
-            disc.mesh_geometry(defaults.geometry)
+            defaults = shared.TestDefaults1DMacro(n)
+            disc = pybamm.FiniteVolumeDiscretisation(defaults.mesh)
             mesh = disc.mesh
 
             combined_submesh = mesh.combine_submeshes(*whole_cell)
@@ -440,12 +415,8 @@ class TestFiniteVolumeDiscretisation(unittest.TestCase):
         # Function for convergence testing
         def get_l2_error(n):
             # create discretisation
-            defaults = shared.TestDefaults1DMacro()
-            defaults.set_equal_pts(n)
-            disc = pybamm.FiniteVolumeDiscretisation(
-                defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-            )
-            disc.mesh_geometry(defaults.geometry)
+            defaults = shared.TestDefaults1DMacro(n)
+            disc = pybamm.FiniteVolumeDiscretisation(defaults.mesh)
             mesh = disc.mesh
 
             # Set up discretisation
@@ -490,12 +461,8 @@ class TestFiniteVolumeDiscretisation(unittest.TestCase):
         def get_l2_error(n):
 
             # create discretisation
-            defaults = shared.TestDefaults1DMacro()
-            defaults.set_equal_pts(n)
-            disc = pybamm.FiniteVolumeDiscretisation(
-                defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-            )
-            disc.mesh_geometry(defaults.geometry)
+            defaults = shared.TestDefaults1DMacro(n)
+            disc = pybamm.FiniteVolumeDiscretisation(defaults.mesh)
             mesh = disc.mesh
 
             whole_cell = ["negative electrode", "separator", "positive electrode"]
@@ -538,12 +505,8 @@ class TestFiniteVolumeDiscretisation(unittest.TestCase):
         def get_l2_error(n):
             whole_cell = ["negative electrode", "separator", "positive electrode"]
             # create discretisation
-            defaults = shared.TestDefaults1DMacro()
-            defaults.set_equal_pts(n)
-            disc = pybamm.FiniteVolumeDiscretisation(
-                defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-            )
-            disc.mesh_geometry(defaults.geometry)
+            defaults = shared.TestDefaults1DMacro(n)
+            disc = pybamm.FiniteVolumeDiscretisation(defaults.mesh)
             mesh = disc.mesh
 
             combined_submesh = mesh.combine_submeshes(*whole_cell)
@@ -585,10 +548,7 @@ class TestFiniteVolumeDiscretisation(unittest.TestCase):
         def get_l2_error(n):
             defaults = shared.TestDefaults1DParticle(n)
 
-            disc = pybamm.FiniteVolumeDiscretisation(
-                defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-            )
-            disc.mesh_geometry(defaults.geometry)
+            disc = pybamm.FiniteVolumeDiscretisation(defaults.mesh)
             mesh = disc.mesh["negative particle"]
             r = mesh.nodes
 
diff --git a/tests/test_expression_tree/test_concatenations.py b/tests/test_expression_tree/test_concatenations.py
index f6f14871..6043a316 100644
--- a/tests/test_expression_tree/test_concatenations.py
+++ b/tests/test_expression_tree/test_concatenations.py
@@ -80,10 +80,7 @@ class TestConcatenations(unittest.TestCase):
     def test_numpy_domain_concatenation(self):
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = shared.DiscretisationForTesting(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = shared.DiscretisationForTesting(defaults.mesh)
         mesh = disc.mesh
 
         a_dom = ["negative electrode"]
diff --git a/tests/test_expression_tree/test_unary_operators.py b/tests/test_expression_tree/test_unary_operators.py
index 4c96162b..f3bb3929 100644
--- a/tests/test_expression_tree/test_unary_operators.py
+++ b/tests/test_expression_tree/test_unary_operators.py
@@ -76,10 +76,7 @@ class TestUnaryOperators(unittest.TestCase):
     def test_numpy_broadcast(self):
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = shared.DiscretisationForTesting(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.mesh_geometry(defaults.geometry)
+        disc = shared.DiscretisationForTesting(defaults.mesh)
         mesh = disc.mesh
 
         whole_cell = ["negative electrode", "separator", "positive electrode"]
diff --git a/tests/test_models/standard_model_tests.py b/tests/test_models/standard_model_tests.py
index 4236cd80..2c303aa2 100644
--- a/tests/test_models/standard_model_tests.py
+++ b/tests/test_models/standard_model_tests.py
@@ -27,7 +27,7 @@ class StandardModelTest(object):
         # Overwrite discretisation if given
         if disc is not None:
             self.disc = disc
-        self.disc.process_model(self.model, self.geometry)
+        self.disc.process_model(self.model)
         # Model should still be well-posed after processing
         self.model.check_well_posedness()
 
diff --git a/tests/test_models/test_interface.py b/tests/test_models/test_interface.py
index 9d6dd7b0..5c3a1db2 100644
--- a/tests/test_models/test_interface.py
+++ b/tests/test_models/test_interface.py
@@ -4,6 +4,7 @@
 from __future__ import absolute_import, division
 from __future__ import print_function, unicode_literals
 import pybamm
+import tests.shared as shared
 
 import unittest
 
@@ -49,32 +50,16 @@ class TestHomogeneousReaction(unittest.TestCase):
         self.assertEqual(processed_rxn.children[2].evaluate() * lp.evaluate(), -1)
 
     def test_discretisation(self):
-        whole_cell = ["negative electrode", "separator", "positive electrode"]
-        param = pybamm.ParameterValues(
-            "input/parameters/lithium-ion/parameters/LCO.csv"
-        )
-        geometry = pybamm.Geometry1DMacro()
-        param.process_geometry(geometry)
+        defaults = shared.TestDefaults1DMacro()
+        disc = pybamm.FiniteVolumeDiscretisation(defaults.mesh)
 
-        mesh_type = pybamm.Mesh
-        submesh_pts = {
-            "negative electrode": {"x": 40},
-            "separator": {"x": 25},
-            "positive electrode": {"x": 35},
-        }
-        submesh_types = {
-            "negative electrode": pybamm.Uniform1DSubMesh,
-            "separator": pybamm.Uniform1DSubMesh,
-            "positive electrode": pybamm.Uniform1DSubMesh,
-        }
-
-        disc = pybamm.FiniteVolumeDiscretisation(mesh_type, submesh_pts, submesh_types)
-        disc._mesh = mesh_type(geometry, submesh_types, submesh_pts)
         rxn = pybamm.interface.homogeneous_reaction()
 
-        param_rxn = param.process_symbol(rxn)
+        param_rxn = defaults.param.process_symbol(rxn)
         processed_rxn = disc.process_symbol(param_rxn)
 
+        whole_cell = ["negative electrode", "separator", "positive electrode"]
+
         combined_submeshes = disc.mesh.combine_submeshes(*whole_cell)
         # processed_rxn should be a vector with the right shape
         self.assertIsInstance(processed_rxn, pybamm.Vector)
diff --git a/tests/test_models/test_lead_acid.py b/tests/test_models/test_lead_acid.py
index 2ccb0e60..9677f20d 100644
--- a/tests/test_models/test_lead_acid.py
+++ b/tests/test_models/test_lead_acid.py
@@ -23,7 +23,7 @@ class TestLeadAcidLOQS(unittest.TestCase):
         # process parameter values, discretise and solve
         model.default_parameter_values.process_model(model)
         disc = model.default_discretisation
-        disc.process_model(model, model.default_geometry)
+        disc.process_model(model)
         t_eval = np.linspace(0, 1, 100)
         solver = model.default_solver
         solver.solve(model, t_eval)
diff --git a/tests/test_models/test_simple_ode_model.py b/tests/test_models/test_simple_ode_model.py
index 904a1dd9..b6523c80 100644
--- a/tests/test_models/test_simple_ode_model.py
+++ b/tests/test_models/test_simple_ode_model.py
@@ -23,9 +23,8 @@ class TestSimpleODEModel(unittest.TestCase):
 
         # discretise and solve
         disc = model.default_discretisation
-        disc.mesh_geometry(model.default_geometry)
         combined_submesh = disc.mesh.combine_submeshes(*whole_cell)
-        disc.process_model(model, model.default_geometry)
+        disc.process_model(model)
         t_eval = np.linspace(0, 1, 100)
         solver = model.default_solver
         solver.solve(model, t_eval)
diff --git a/tests/test_solvers/test_scikits_solvers.py b/tests/test_solvers/test_scikits_solvers.py
index 59ab7d47..6d0e37b8 100644
--- a/tests/test_solvers/test_scikits_solvers.py
+++ b/tests/test_solvers/test_scikits_solvers.py
@@ -74,7 +74,7 @@ class TestScikitsSolver(unittest.TestCase):
         model.initial_conditions = {var: 1}
         model.initial_conditions_ydot = {var: 0.1}
         disc = model.default_discretisation
-        disc.process_model(model, model.default_geometry)
+        disc.process_model(model)
 
         # Solve
         solver = pybamm.ScikitsOdeSolver(tol=1e-9)
@@ -94,7 +94,7 @@ class TestScikitsSolver(unittest.TestCase):
         model.initial_conditions = {var1: 1, var2: 2}
         model.initial_conditions_ydot = {var1: 0.1, var2: 0.2}
         disc = model.default_discretisation
-        disc.process_model(model, model.default_geometry)
+        disc.process_model(model)
 
         # Solve
         solver = pybamm.ScikitsDaeSolver(tol=1e-8)
diff --git a/tests/test_solvers/test_scipy_solver.py b/tests/test_solvers/test_scipy_solver.py
index b7d7960f..599f3c84 100644
--- a/tests/test_solvers/test_scipy_solver.py
+++ b/tests/test_solvers/test_scipy_solver.py
@@ -46,10 +46,8 @@ class TestScipySolver(unittest.TestCase):
 
         # create discretisation
         defaults = shared.TestDefaults1DMacro()
-        disc = shared.DiscretisationForTesting(
-            defaults.mesh_type, defaults.submesh_pts, defaults.submesh_types
-        )
-        disc.process_model(model, defaults.geometry)
+        disc = shared.DiscretisationForTesting(defaults.mesh)
+        disc.process_model(model)
         # Solve
         solver = pybamm.ScipySolver(tol=1e-8, method="RK45")
         t_eval = np.linspace(0, 1, 100)

pip install -e .[dev,docs]
parse_log_pytest
requirements.txt
pytest
apt-get update apt-get install -y gcc gfortran libopenblas-dev
3.6
requirements.txt
pytest --no-header -rA --tb=line --color=no -p no:cacheprovider -W ignore::DeprecationWarning
alabaster==0.7.13
anyio==3.6.2
anytree==2.8.0
argon2-cffi==21.3.0
argon2-cffi-bindings==21.2.0
async-generator==1.10
attrs==22.2.0
Babel==2.11.0
backcall==0.2.0
bleach==4.1.0
certifi==2021.5.30
cffi==1.15.1
charset-normalizer==2.0.12
comm==0.1.4
contextvars==2.4
cycler==0.11.0
dataclasses==0.8
decorator==5.1.1
defusedxml==0.7.1
docutils==0.17.1
entrypoints==0.4
flake8==5.0.4
guzzle-sphinx-theme==0.7.11
idna==3.10
imagesize==1.4.1
immutables==0.19
importlib-metadata==4.2.0
iniconfig==1.1.1
ipykernel==5.5.6
ipython==7.16.3
ipython-genutils==0.2.0
ipywidgets==7.8.5
jedi==0.17.2
Jinja2==3.0.3
json5==0.9.16
jsonschema==3.2.0
jupyter==1.1.1
jupyter-client==7.1.2
jupyter-console==6.4.3
jupyter-core==4.9.2
jupyter-server==1.13.1
jupyterlab==3.2.9
jupyterlab-pygments==0.1.2
jupyterlab-server==2.10.3
jupyterlab_widgets==1.1.11
kiwisolver==1.3.1
MarkupSafe==2.0.1
matplotlib==3.3.4
mccabe==0.7.0
mistune==0.8.4
nbclassic==0.3.5
nbclient==0.5.9
nbconvert==6.0.7
nbformat==5.1.3
nest-asyncio==1.6.0
notebook==6.4.10
numpy==1.19.5
packaging==21.3
pandas==1.1.5
pandocfilters==1.5.1
parso==0.7.1
pexpect==4.9.0
pickleshare==0.7.5
Pillow==8.4.0
pluggy==1.0.0
prometheus-client==0.17.1
prompt-toolkit==3.0.36
ptyprocess==0.7.0
py==1.11.0
-e git+https://github.com/tinosulzer/PyBaMM.git@fd0bc6357f40c988ae992a35ca12f0d1ce5e3d9a#egg=pybamm
pycodestyle==2.9.1
pycparser==2.21
pyflakes==2.5.0
Pygments==2.14.0
pyparsing==3.1.4
pyrsistent==0.18.0
pytest==7.0.1
python-dateutil==2.9.0.post0
pytz==2025.2
pyzmq==25.1.2
requests==2.27.1
scipy==1.5.4
Send2Trash==1.8.3
six==1.17.0
sniffio==1.2.0
snowballstemmer==2.2.0
Sphinx==4.3.2
sphinxcontrib-applehelp==1.0.2
sphinxcontrib-devhelp==1.0.2
sphinxcontrib-htmlhelp==2.0.0
sphinxcontrib-jsmath==1.0.1
sphinxcontrib-qthelp==1.0.3
sphinxcontrib-serializinghtml==1.1.5
terminado==0.12.1
testpath==0.6.0
tomli==1.2.3
tornado==6.1
traitlets==4.3.3
typing_extensions==4.1.1
urllib3==1.26.20
wcwidth==0.2.13
webencodings==0.5.1
websocket-client==1.3.1
widgetsnbextension==3.6.10
zipp==3.6.0

name: PyBaMM
channels:
  - defaults
  - https://repo.anaconda.com/pkgs/main
  - https://repo.anaconda.com/pkgs/r
  - conda-forge
dependencies:
  - _libgcc_mutex=0.1=main
  - _openmp_mutex=5.1=1_gnu
  - ca-certificates=2025.2.25=h06a4308_0
  - certifi=2021.5.30=py36h06a4308_0
  - ld_impl_linux-64=2.40=h12ee557_0
  - libffi=3.3=he6710b0_2
  - libgcc-ng=11.2.0=h1234567_1
  - libgomp=11.2.0=h1234567_1
  - libstdcxx-ng=11.2.0=h1234567_1
  - ncurses=6.4=h6a678d5_0
  - openssl=1.1.1w=h7f8727e_0
  - pip=21.2.2=py36h06a4308_0
  - python=3.6.13=h12debd9_1
  - readline=8.2=h5eee18b_0
  - setuptools=58.0.4=py36h06a4308_0
  - sqlite=3.45.3=h5eee18b_0
  - tk=8.6.14=h39e8969_0
  - wheel=0.37.1=pyhd3eb1b0_0
  - xz=5.6.4=h5eee18b_1
  - zlib=1.2.13=h5eee18b_1
  - pip:
      - alabaster==0.7.13
      - anyio==3.6.2
      - anytree==2.8.0
      - argon2-cffi==21.3.0
      - argon2-cffi-bindings==21.2.0
      - async-generator==1.10
      - attrs==22.2.0
      - babel==2.11.0
      - backcall==0.2.0
      - bleach==4.1.0
      - cffi==1.15.1
      - charset-normalizer==2.0.12
      - comm==0.1.4
      - contextvars==2.4
      - cycler==0.11.0
      - dataclasses==0.8
      - decorator==5.1.1
      - defusedxml==0.7.1
      - docutils==0.17.1
      - entrypoints==0.4
      - flake8==5.0.4
      - guzzle-sphinx-theme==0.7.11
      - idna==3.10
      - imagesize==1.4.1
      - immutables==0.19
      - importlib-metadata==4.2.0
      - iniconfig==1.1.1
      - ipykernel==5.5.6
      - ipython==7.16.3
      - ipython-genutils==0.2.0
      - ipywidgets==7.8.5
      - jedi==0.17.2
      - jinja2==3.0.3
      - json5==0.9.16
      - jsonschema==3.2.0
      - jupyter==1.1.1
      - jupyter-client==7.1.2
      - jupyter-console==6.4.3
      - jupyter-core==4.9.2
      - jupyter-server==1.13.1
      - jupyterlab==3.2.9
      - jupyterlab-pygments==0.1.2
      - jupyterlab-server==2.10.3
      - jupyterlab-widgets==1.1.11
      - kiwisolver==1.3.1
      - markupsafe==2.0.1
      - matplotlib==3.3.4
      - mccabe==0.7.0
      - mistune==0.8.4
      - nbclassic==0.3.5
      - nbclient==0.5.9
      - nbconvert==6.0.7
      - nbformat==5.1.3
      - nest-asyncio==1.6.0
      - notebook==6.4.10
      - numpy==1.19.5
      - packaging==21.3
      - pandas==1.1.5
      - pandocfilters==1.5.1
      - parso==0.7.1
      - pexpect==4.9.0
      - pickleshare==0.7.5
      - pillow==8.4.0
      - pluggy==1.0.0
      - prometheus-client==0.17.1
      - prompt-toolkit==3.0.36
      - ptyprocess==0.7.0
      - py==1.11.0
      - pycodestyle==2.9.1
      - pycparser==2.21
      - pyflakes==2.5.0
      - pygments==2.14.0
      - pyparsing==3.1.4
      - pyrsistent==0.18.0
      - pytest==7.0.1
      - python-dateutil==2.9.0.post0
      - pytz==2025.2
      - pyzmq==25.1.2
      - requests==2.27.1
      - scipy==1.5.4
      - send2trash==1.8.3
      - six==1.17.0
      - sniffio==1.2.0
      - snowballstemmer==2.2.0
      - sphinx==4.3.2
      - sphinxcontrib-applehelp==1.0.2
      - sphinxcontrib-devhelp==1.0.2
      - sphinxcontrib-htmlhelp==2.0.0
      - sphinxcontrib-jsmath==1.0.1
      - sphinxcontrib-qthelp==1.0.3
      - sphinxcontrib-serializinghtml==1.1.5
      - terminado==0.12.1
      - testpath==0.6.0
      - tomli==1.2.3
      - tornado==6.1
      - traitlets==4.3.3
      - typing-extensions==4.1.1
      - urllib3==1.26.20
      - wcwidth==0.2.13
      - webencodings==0.5.1
      - websocket-client==1.3.1
      - widgetsnbextension==3.6.10
      - zipp==3.6.0
prefix: /opt/conda/envs/PyBaMM

tests/test_discretisations/test_base_discretisations.py::TestDiscretise::test_broadcast tests/test_discretisations/test_base_discretisations.py::TestDiscretise::test_concatenate_init tests/test_discretisations/test_base_discretisations.py::TestDiscretise::test_concatenation tests/test_discretisations/test_base_discretisations.py::TestDiscretise::test_concatenation_of_scalars tests/test_discretisations/test_base_discretisations.py::TestDiscretise::test_core_NotImplementedErrors tests/test_discretisations/test_base_discretisations.py::TestDiscretise::test_discretise_slicing tests/test_discretisations/test_base_discretisations.py::TestDiscretise::test_discretise_space tests/test_discretisations/test_base_discretisations.py::TestDiscretise::test_discretise_spatial_operator tests/test_discretisations/test_base_discretisations.py::TestDiscretise::test_process_complex_expression tests/test_discretisations/test_base_discretisations.py::TestDiscretise::test_process_dict tests/test_discretisations/test_base_discretisations.py::TestDiscretise::test_process_model_dae tests/test_discretisations/test_base_discretisations.py::TestDiscretise::test_process_model_ode tests/test_discretisations/test_base_discretisations.py::TestDiscretise::test_process_symbol_base tests/test_discretisations/test_finite_volume_discretisations.py::TestFiniteVolumeDiscretisation::test_add_ghost_nodes tests/test_discretisations/test_finite_volume_discretisations.py::TestFiniteVolumeDiscretisation::test_discretise_diffusivity_times_spatial_operator tests/test_discretisations/test_finite_volume_discretisations.py::TestFiniteVolumeDiscretisation::test_div_convergence tests/test_discretisations/test_finite_volume_discretisations.py::TestFiniteVolumeDiscretisation::test_div_convergence_internal tests/test_discretisations/test_finite_volume_discretisations.py::TestFiniteVolumeDiscretisation::test_grad_convergence_with_bcs tests/test_discretisations/test_finite_volume_discretisations.py::TestFiniteVolumeDiscretisation::test_grad_convergence_without_bcs tests/test_discretisations/test_finite_volume_discretisations.py::TestFiniteVolumeDiscretisation::test_grad_div_shapes_Dirichlet_bcs tests/test_discretisations/test_finite_volume_discretisations.py::TestFiniteVolumeDiscretisation::test_grad_div_shapes_Neumann_bcs tests/test_discretisations/test_finite_volume_discretisations.py::TestFiniteVolumeDiscretisation::test_grad_div_shapes_mixed_domain tests/test_discretisations/test_finite_volume_discretisations.py::TestFiniteVolumeDiscretisation::test_spherical_grad_div_shapes_Dirichlet_bcs tests/test_discretisations/test_finite_volume_discretisations.py::TestFiniteVolumeDiscretisation::test_spherical_grad_div_shapes_Neumann_bcs tests/test_discretisations/test_finite_volume_discretisations.py::TestFiniteVolumeDiscretisation::test_spherical_operators tests/test_expression_tree/test_concatenations.py::TestConcatenations::test_numpy_domain_concatenation tests/test_expression_tree/test_unary_operators.py::TestUnaryOperators::test_numpy_broadcast tests/test_models/test_interface.py::TestHomogeneousReaction::test_discretisation tests/test_models/test_lead_acid.py::TestLeadAcidLOQS::test_basic_processing tests/test_models/test_lead_acid.py::TestLeadAcidLOQS::test_solution tests/test_models/test_simple_ode_model.py::TestSimpleODEModel::test_basic_processing tests/test_models/test_simple_ode_model.py::TestSimpleODEModel::test_solution tests/test_solvers/test_scipy_solver.py::TestScipySolver::test_model_solver
tests/test_discretisations/test_finite_volume_discretisations.py::TestFiniteVolumeDiscretisation::test_node_to_edge tests/test_expression_tree/test_concatenations.py::TestConcatenations::test_base_concatenation tests/test_expression_tree/test_concatenations.py::TestConcatenations::test_concatenation_domains tests/test_expression_tree/test_concatenations.py::TestConcatenations::test_concatenation_orphans tests/test_expression_tree/test_concatenations.py::TestConcatenations::test_numpy_model_concatenation_vector_scalar tests/test_expression_tree/test_concatenations.py::TestConcatenations::test_numpy_model_concatenation_vectors tests/test_expression_tree/test_unary_operators.py::TestUnaryOperators::test_absolute tests/test_expression_tree/test_unary_operators.py::TestUnaryOperators::test_broadcast tests/test_expression_tree/test_unary_operators.py::TestUnaryOperators::test_gradient tests/test_expression_tree/test_unary_operators.py::TestUnaryOperators::test_id tests/test_expression_tree/test_unary_operators.py::TestUnaryOperators::test_negation tests/test_expression_tree/test_unary_operators.py::TestUnaryOperators::test_printing tests/test_expression_tree/test_unary_operators.py::TestUnaryOperators::test_unary_operator tests/test_models/test_interface.py::TestHomogeneousReaction::test_set_parameters tests/test_solvers/test_scipy_solver.py::TestScipySolver::test_integrate
pybamm/__init__.py docs/source/geometry/geometry.rst pybamm/discretisations/meshes.py docs/source/meshes/meshes.rst pybamm/discretisations/finite_volume_discretisations.py pybamm/meshes/submeshes.py docs/source/meshes/submeshes.rst docs/source/geometry/index.rst docs/source/meshes/index.rst pybamm/models/base_model.py docs/index.rst pybamm/discretisations/base_discretisation.py
pybamm/__init__.py docs/source/geometry/geometry.rst docs/source/meshes/meshes.rst pybamm/discretisations/finite_volume_discretisations.py pybamm/meshes/meshes.py pybamm/meshes/submeshes.py docs/source/meshes/submeshes.rst docs/source/geometry/index.rst docs/source/meshes/index.rst pybamm/models/base_model.py docs/index.rst pybamm/discretisations/base_discretisation.py