second-order IMEX, static mesh

growth_eqns.py

@@ -24,9 +24,6 @@ class Growth(object):
         self.fminus1 = df.Function(self.function_space)           # f at t = -1
         self.f0 = df.Function(self.function_space)                 # f at t =0
 
-        # df.plot(self.mesh)
-        # plt.show()
-
     def advection(self,c,v,tc):
         return (df.inner((v*c).dx(0),tc))
 
@@ -123,9 +120,6 @@ class Growth(object):
             c_array[i] = c.compute_vertex_values(self.mesh)
             self.f0.assign(self.f)
             self.fminus1.assign(self.f0)
-            # df.ALE.move(self.mesh,u)
-        # df.plot(self.mesh)
-        # plt.show()
         return (u_array, v_array, rho_array, c_array)