From a486e61f9b940c43d27bfc24972ea8e2f8f02c05 Mon Sep 17 00:00:00 2001
From: Manos Katsomallos <manos@delkappa.com>
Date: Fri, 7 Jan 2022 06:17:37 +0100
Subject: [PATCH] code: Testing

---
 code/lib/lmdk_lib.py        | 17 +++++++++++++++++
 code/plots/laplace-graph.py | 22 ++++++++++++++++------
 2 files changed, 33 insertions(+), 6 deletions(-)

diff --git a/code/lib/lmdk_lib.py b/code/lib/lmdk_lib.py
index 80af98a..1072461 100644
--- a/code/lib/lmdk_lib.py
+++ b/code/lib/lmdk_lib.py
@@ -1072,3 +1072,20 @@ def is_landmark(p, lmdks):
   if len(lmdks) > 0 and any(np.equal(lmdks, p).all(1)):
     return True
   return False
+
+
+if __name__ == '__main__':
+  try:
+    seq = [1, 2, 3, 4, 5, 6, 7, 8]
+    lmdk = [1, 3, 5, 8]
+    print(get_hist(seq, lmdk))
+    start_time = time.process_time()
+    # main(args)
+    end_time = time.process_time()
+    print('##############################\n'
+          'Time       : %.4fs\n'
+          '##############################\n'
+          %(end_time - start_time), flush=True)
+  except KeyboardInterrupt:
+    print('Interrupted by user.')
+    exit()
diff --git a/code/plots/laplace-graph.py b/code/plots/laplace-graph.py
index fe60491..cb8a683 100644
--- a/code/plots/laplace-graph.py
+++ b/code/plots/laplace-graph.py
@@ -10,18 +10,28 @@ from scipy.stats import laplace
 # Initialize plot
 lmdk_lib.plot_init()
 # The x axis
-x = np.arange(0.0, 6.0, 0.01)
-plt.xlabel('Count')
+x = np.arange(-3, 3.1, .1)
+plt.xlabel('Noise')
 # Set y axis label
 plt.ylabel('Likelihood')
 # Draw the first laplace distribution
 plt.plot(x,
-         laplace.pdf(x, loc=2.0, scale=1.0),
-         label='Laplace(2, 1)',
-         linewidth=lmdk_lib.line_width)
+         laplace.pdf(x, loc=0, scale=1/1),
+         label='Lap(0, 1)',
+         linewidth=lmdk_lib.line_width*2)
+
+plt.plot(x,
+         laplace.pdf(x, loc=0, scale=1/.5),
+         label='Lap(0, 2)',
+         linewidth=lmdk_lib.line_width*2)
+
+plt.plot(x,
+         laplace.pdf(x, loc=0, scale=1/.25),
+         label='Lap(0, 4)',
+         linewidth=lmdk_lib.line_width*2)
 
 # Configure the plot
-plt.axis([0.0, 6.0, 0.0, 0.6])  # Set plot box.
+plt.axis([-3, 3, 0.0, 0.6])  # Set plot box.
 # Plot legend
 lmdk_lib.plot_legend()
 # Show the plot in a new window