Added timings to LUStability scripts

drreynolds · Feb 10, 2021 · 0f4de61 · 0f4de61
1 parent a65a4f8
commit 0f4de61
Show file tree

Hide file tree

Showing 2 changed files with 13 additions and 0 deletions.
diff --git a/LUFactorization/LUStability.m b/LUFactorization/LUStability.m
@@ -24,17 +24,23 @@
 
    % test LUPFactors_simple
    fprintf('  LUPFactors_simple:\n');
+   ts = tic;
    [L,U,P] = LUPFactors_simple(A);
+   fprintf('    runtime = %g\n', toc(ts))
    fprintf('    norm(A - P^T L U) = %9.2e\n', norm(A - P'*L*U));
 
    % test LUPFactors
    fprintf('  LUPFactors:\n');
+   ts = tic;
    [L,U,P] = LUPFactors(A);
+   fprintf('    runtime = %g\n', toc(ts))
    fprintf('    norm(A - P^T L U) = %9.2e\n', norm(A - P'*L*U));
 
    % test LUPPFactors
    fprintf('  LUPPFactors:\n');
+   ts = tic;
    [L,U,P1,P2] = LUPPFactors(A);
+   fprintf('    runtime = %g\n', toc(ts))
    fprintf('    norm(A - P1^T L U P2^T) = %9.2e\n', norm(A - P1'*L*U*P2'));
 
 end
diff --git a/LUFactorization/LUStability.py b/LUFactorization/LUStability.py
@@ -9,6 +9,7 @@
 
 # imports
 import numpy
+import time
 from LUPFactors_simple import LUPFactors_simple
 from LUPFactors import LUPFactors
 from LUPPFactors import LUPPFactors
@@ -31,17 +32,23 @@
 
   # test LUPFactors_simple
   print("  LUPFactors_simple:")
+  ts = time.time()
   L, U, P = LUPFactors_simple(A)
+  print("    runtime = ", time.time()-ts)
   print("    norm(A - P^T L U) = ", numpy.linalg.norm(A - P.T@L@U))
 
   # test LUPFactors
   print("  LUPFactors:")
+  ts = time.time()
   L, U, P = LUPFactors(A)
+  print("    runtime = ", time.time()-ts)
   print("    norm(A - P^T L U) = ", numpy.linalg.norm(A - P.T@L@U))
 
   # test LUPPFactors
   print("  LUPPFactors:")
+  ts = time.time()
   L, U, P1, P2 = LUPPFactors(A)
+  print("    runtime = ", time.time()-ts)
   print("    norm(A - P1^T L U P2^T) = ", numpy.linalg.norm(A - P1.T@L@U@P2.T))
 
 # end of script