updated version of code for pseudo spec paper

Parallel-in-Time · Jan 3, 2025 · 5a1f55b · 5a1f55b
1 parent 183d283
commit 5a1f55b
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/scripts/pseudo-spectrum/figure_13-16.py b/scripts/pseudo-spectrum/figure_13-16.py
@@ -84,6 +84,7 @@
 
 
 Pmat, Bmat = para.get_parareal_matrix()
+print("|| E ||_2 = %5.2f" % np.linalg.norm(Pmat.todense(),2))
 
 defect_inf = np.zeros((1,maxiter))
 defect_l2  = np.zeros((1,maxiter))

diff --git a/scripts/pseudo-spectrum/figure_9-12.py b/scripts/pseudo-spectrum/figure_9-12.py
@@ -82,7 +82,7 @@
 
 
 Pmat, Bmat = para.get_parareal_matrix()
-
+print("|| E ||_2 = %5.2f" % np.linalg.norm(Pmat.todense(), 2))
 nreal = 20
 nimag = 20
 lambda_real = np.linspace(-3.0, 3.0,  nreal)
@@ -96,7 +96,6 @@
 Also, the eps-isolines are very much circles.
 QUESTION: can we work out the D matrix above and say something about how it looks like for diffusive/non-diffusive problems?
 '''
-print("2-Norm of E: %5.3f" % np.linalg.norm(Pmat.todense(),2))
 for i in range(0,nreal):
   for j in range(0,nimag):
     z = lambda_real[i] + 1j*lambda_imag[j]