diff --git a/LDDsimulation/LDDsimulation.py b/LDDsimulation/LDDsimulation.py
index ede0b8af4d58b4edd256882b06538b15f7ab7b0b..662d5c41476fa35927182fe924f307bf8553d2a7 100644
--- a/LDDsimulation/LDDsimulation.py
+++ b/LDDsimulation/LDDsimulation.py
@@ -1050,17 +1050,7 @@ class LDDsimulation(object):
                 subdomain.pressure_exact = dict()
                 for phase in subdomain.has_phases:
                     print(f"initialising exact solutions for subdomain{subdomain.subdomain_index} and phase {phase}")
-                    # if np.abs(self.t-self.starttime) < self.calc_tol:
-                    #     # here t = 0 and we have to initialise the sources.
-
-                    # p0 contains both pw_0 and pnw_0 for subdomain[subdomain_index]
-                    # print(f"subdomain{subdomain.subdomain_index} has phases = {subdomain.has_phases}")
                     pa_exact = self.exact_solution[sd_ind][phase]
-                    # self.exact_solution_expr.update({sd_ind: dict()})
-                    # note that the default interpolation degree is 2
-                    # self.exact_solution_expr[sd_ind].update(
-                    #     {phase: df.Expression(pa_exact, domain = mesh, degree=self.interpolation_degree, t=self.starttime) }
-                    #     )
                     subdomain.pressure_exact.update(
                         {phase: df.Expression(pa_exact, domain = mesh, degree=self.interpolation_degree, t=self.starttime) }
                         )