align with paper

Author: pancetta

pySDC/implementations/problem_classes/AllenCahn_2D_FD.py

@@ -175,7 +175,7 @@ def u_exact(self, t):
         for i in range(self.params.nvars[0]):
             for j in range(self.params.nvars[1]):
                 r2 = self.xvalues[i] ** 2 + self.xvalues[j] ** 2
-                me.values[i, j] = np.tanh((self.params.radius - np.sqrt(r2)) / self.params.eps)
+                me.values[i, j] = np.tanh((self.params.radius - np.sqrt(r2)) / (np.sqrt(2) * self.params.eps))
 
         return me
 

pySDC/projects/TOMS/AllenCahn_contracting_circle.py

@@ -294,7 +294,7 @@ def show_results(fname, cwd=''):
     ax.plot(xcoords, init_width, color='k', linestyle='--', linewidth=1, label='exact')
 
     ax.yaxis.set_major_formatter(ticker.FormatStrFormatter('%1.2f'))
-    ax.set_ylabel('interface')
+    ax.set_ylabel(r'interface width ($\epsilon$)')
     ax.set_xlabel('time')
     ax.grid()
     ax.legend(loc=3)
@@ -318,19 +318,19 @@ def main(cwd=''):
         cwd (str): current working directory (need this for testing)
     """
 
-    # Loop over variants, exact and inexact solves
-    results = {}
-    for variant in ['multi-implicit', 'semi-implicit', 'fully-implicit', 'semi-implicit_v2', 'multi-implicit_v2']:
-
-        results[(variant, 'exact')] = run_SDC_variant(variant=variant, inexact=False)
-        results[(variant, 'inexact')] = run_SDC_variant(variant=variant, inexact=True)
-
-    # dump result
+    # # Loop over variants, exact and inexact solves
+    # results = {}
+    # for variant in ['multi-implicit', 'semi-implicit', 'fully-implicit', 'semi-implicit_v2', 'multi-implicit_v2']:
+    #
+    #     results[(variant, 'exact')] = run_SDC_variant(variant=variant, inexact=False)
+    #     results[(variant, 'inexact')] = run_SDC_variant(variant=variant, inexact=True)
+    #
+    # # dump result
     fname = 'data/results_SDC_variants_AllenCahn_1E-03'
-    file = open(cwd + fname + '.pkl', 'wb')
-    dill.dump(results, file)
-    file.close()
-    assert os.path.isfile(cwd + fname + '.pkl'), 'ERROR: dill did not create file'
+    # file = open(cwd + fname + '.pkl', 'wb')
+    # dill.dump(results, file)
+    # file.close()
+    # assert os.path.isfile(cwd + fname + '.pkl'), 'ERROR: dill did not create file'
 
     # visualize
     show_results(fname, cwd=cwd)