From 2e1c170af6a764ffc5e2557d8e7e9f2ba0fd681f Mon Sep 17 00:00:00 2001
From: Apostolof <apotwohd@gmail.com>
Date: Sat, 15 Sep 2018 12:13:11 +0300
Subject: [PATCH] Fix serial implementation Add extension to executable
 filename so git ignore can catch it Fix matrix transpose function

---
 .gitignore                                 |   3 +
 datasets/smallset                          |   5 +
 datasets/toronto_death_penalty_refined_adj |   2 +-
 serial/Makefile                            |   5 +-
 serial/pagerank                            | Bin 25344 -> 0 bytes
 serial/serial_gs_pagerank.c                |   4 +-
 serial/serial_gs_pagerank_functions.c      | 427 +++++++++++++--------
 serial/serial_gs_pagerank_functions.h      |  52 ++-
 8 files changed, 304 insertions(+), 194 deletions(-)
 create mode 100644 datasets/smallset
 delete mode 100644 serial/pagerank

diff --git a/.gitignore b/.gitignore
index c6127b3..25734a9 100644
--- a/.gitignore
+++ b/.gitignore
@@ -50,3 +50,6 @@ modules.order
 Module.symvers
 Mkfile.old
 dkms.conf
+
+#Output files
+*_output
\ No newline at end of file
diff --git a/datasets/smallset b/datasets/smallset
new file mode 100644
index 0000000..0fd3edd
--- /dev/null
+++ b/datasets/smallset
@@ -0,0 +1,5 @@
+0 0 0 0 0
+1 0 1 0 0
+1 1 0 1 1
+0 0 0 0 1
+0 0 0 1 0
\ No newline at end of file
diff --git a/datasets/toronto_death_penalty_refined_adj b/datasets/toronto_death_penalty_refined_adj
index 3b4fc82..e38c650 100644
--- a/datasets/toronto_death_penalty_refined_adj
+++ b/datasets/toronto_death_penalty_refined_adj
@@ -1797,4 +1797,4 @@
 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
-0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
\ No newline at end of file
diff --git a/serial/Makefile b/serial/Makefile
index d6b3229..78d2e2e 100644
--- a/serial/Makefile
+++ b/serial/Makefile
@@ -5,14 +5,15 @@ SHELL := /bin/bash
 
 CC = gcc
 RM = rm -f
-CFLAGS=-O0 -g -I.
+CFLAGS_DEBUG=-O0 -g -I.
+CFLAGS=-O3 -I.
 OBJ=serial_gs_pagerank.o serial_gs_pagerank_functions.o
 DEPS=serial_gs_pagerank_functions.h
 
 # ==========================================
 # TARGETS
 
-EXECUTABLES = pagerank
+EXECUTABLES = pagerank.out
 
 .PHONY: all clean
 
diff --git a/serial/pagerank b/serial/pagerank
deleted file mode 100644
index 94f5af26eacb5ef9679323f48b26c29b3d3bfa33..0000000000000000000000000000000000000000
GIT binary patch
literal 0
HcmV?d00001

literal 25344
zcmeHve|%KcweLPNlgSwZnF#}!h}Rzc3D*2TQ4ubf5DrW%fe1mdg2QBHl8huX>HG*_
zt5E~iX`~RZ*ZTTsFZH&)^yk&y*2@*!2C;(LcRy+Cv)tDEKB+}J5v?iyxT&J^zH9%O
zIWtMr_df5x$Id5b@3r<?Yp=ET+WVZn&z{p~u32T%G^UcBEoTJPPZOAWSyBI4SwL3L
z7Bd&F7qSbO6R?BhM18%;P%F~76s>Yy;_b+ZZ?R1viu3a%oywJll$Xf5$+G%+0;_Vy
z42~c^(sgk)2ay}Rpz?94M;$4~lfkcqCcYNQ*CP2;K1PC(RIcnJ+vvYl$}P@U5)+cj
ziI!5!E6EU*7oX3;3HgHp%~Y=VZU7(Y`AL&U60TS6P2rd0Qh$|8b6IUyr0wcOwOyg=
zu1GYsqk6~UtE;bGRFjC+T*2)o`9!;V?MB`bxmAkd^|JsIj_R8;9pf90?(TT0;`G1N
z{r2RI&)@l?hj`ULik~?B29(Am)N6s8ISGEw6!;}m;KNhk=T3pYd<y)PQ{WR*;C)lz
zADRNccMANQz#I54P7Odz=AT=qz;{i7r+%MI&WFGo_%BXnKuktIQ=6RrF3?|4EQDah
zW01n{Y8KuRNwW4>PdLhYQpp7K`xD7PaGO8axy|1mh;*@@cqE!^XYug%NHoL}$#^m*
zvo6-&5l#ld5=q2@*Iw&yPXq(ecGey44)*M1?ZK{CBFwr2U0tyt17$MO9gek!0z1JG
z+!64%N1}nQ$n9Yg8tm=?meiv5<4J#aAQEM(*EFqY^j}f4a8a(9FI`b{RStU%^Eb6N
z`$OS)xFeEChU2Zxja{*5xHZt$6^7uB?pRa`_A~$Wz5ex*uCXhSNQ7ZFu}{iO?y_O@
z(Vs)+0-`Zl!ZBiW+L;ClbOlrM#e;(nVVWsn2PJM+Br=0$x_^;y>fDHncz-J4b9lXt
zJqwtUns3!SSx*6gR4?aMg=Z44*BQKC;T7yD!f&Mkh*VjG4@$gIgm0Jl`9=7M#1GdR
zTprD-WWP!$EO@CRuNtx7;gDR)Sn$?3AG6?TE>vmUf}gFB$TvxRDavXdrRuE~eD2l-
ze5VDk?h!=aW5HYRqdP2k>wV-d3!d6m={^fy^$peSv)~a3x%99FFCQ&<-2n@JdXAd0
zgBE<b1^=i8@3!C%S@0@Wsd~tQudvV`x8S`Nys!Ut-8WG3i%SgV+dY)DWruzJPwPi>
z6xpk0qAGjwEL`26sz-_7PO2Xr$s%1mjbI9y(P0kT2&OI@9pdnN7XzlQ8hw<*rwOJm
z8a=?_*9fMr8QsU>mkFjW8NG|czazMm;2sYDnqcaJ(XAYQj$m^AXbXp*BA8r0TF>Dp
z2&T{-t>f^I2qxE$8XW#E!4$%y9u7Z3a5=$@!z8bLA7GvqzR{Ha)vdnt%f9}TV=b*~
z_LSTP_xtwDcniJb+p}ysJdvFp1<%{=Pp#yre*@s#Q$pEB--2;p`W0XQczN$^xP3WQ
z!R@Xs*^zcmzvU=-xxL2ysa?E1h}_BfuZj!)4~hTChjw3j%y%Snolkql_p1++l@L}T
zg-sX21b?x-yOu2g8B1O2>tFV#i%{AKdCBR%fn~1&%FKYt86PfNo+)7<)ZUt?k9Gl+
z?!S(u%0_PmrO+Ggry+oh86}a~;Vs<%R7OPImA+%YFP+*K@(nmHo)61>X>)wxkZ;e%
zVc$UWz5x@p|DHaUxfky9?P(rH1u9dAe0$7CeZcplgMUU2fj_97%*hY3SwUnTq8QGm
z#<HnQHZ_t>oyew!ed!m8Z4~{C*6Y&d_uY@1nO#H4IRobR)6YQ4^XTG3Br&rSqLl2v
zA=&ql?19-;WY0j!1(Yp#%9kE)cfUDw+ULG(=(O82qYV9h+I{PhcK7__pgMUWq9k9x
z=?HKm0{1NO96#;eeB`vd!UQq?)8>5;V&3Q5b5$IQK%Z}rJ~WInqi3>Nt{FTt9`ao|
zojWX9k{Lq9vqyzTU-$UZcZ@?8;)%?8FlWwpLG^Pks19~0>hv9F(kakx92+ptEH%ed
zbHF=BKKefTP|!2b_nT)lP~@6vN#F5Hblea)4;6J@{~f>7?sVch?M^p@4fzJptH&~%
zArIi8W8_>o5T$yDD#GZJZ(zW{6ZM#P%tM5DSRy=xc&5IR9W!4b14f_3G=h%Ugn%0`
zU+_Kb+cW(aQqqu5gD24-@8@^l&t`WIxp#k2v<Ou^1I@n#4-cyG53|{{`GOQN%OV27
z2j6czA;nw=G3`G0ov_x#czkoH-Cbjj+&(W?gxB^gtAeXD4v-43J%^IC#2n#vjnW+p
zdLNwkae6P8dfzjt-t|K7pV4GCrSi=u(x)+Y$n%UeNt+p8{~Z~Y`V^EpLF18hz6TRr
zrsX}Q^BkdbCjAs?o^ElBHq*cxnEiqAaVif$2)X*eq$=-02z??~zwp<*tH+do$8Imn
z6_Z{XAY-`9=z+u06GwUrOnM$;J97ubh;C?vUaA|wRZ;Ib^#vHT=wl6f0wQ(~rKS-x
zf=@eT4w%2}&uB}{7gF;kcRLSMU%JIORRu0Y8xQ<oZs7eZMD`<w!M)V{rThNH@X9_=
z(r6abA4c=(aW$GZU;<l=;f!ek1C0x;H!w3Wn!V7kn;<UVuhgBRxSwV2h5%vvY7iVV
zPoPKGDNp|$C$q_UxK2ws`ppv>1ZGIeQy@>9CsDytzdxmsK11%u>np{ic!JL&--28*
zDTY8z6ayPiq{&<^!mp4FWAVrc?m@>g-y<B4BdG5=X1+vPj1J%eE*hc(8(;DTf9*T+
zvfa1mJQ6hocM-dd*t<Bp{g^pIe5Y1UWYUJYkP#uo=F=!<GjNX=%0h(H_rItYlNJ5u
z2w99amYO4dkHOA<^Q9c-CHLNTF+B-;z_btdhsC#S;nOO0#C&uF=7ZGmfa%Fw`3q1A
zD^Uk4xwfO{@<D@9MpjA%xs?UHY0e81Pm+np%nW6xst3#wn$dF8(XdZin#ox@Jis*c
zn3+K&d6t@)zUPmbuM_tuJ&*97G+}H}`%=v7%In(jz(x;^tJChwo}%$&z;U>B9mlPU
zZb|G^!zU3UC735k#D8K=vWn<8Ul%&Q?%o@Mhp8L5&Li%<zr!>}UY;+#3<qb9nd97!
zui<*kJj2U(<I201yg8oNo_qI94ki<t4Mo@t!e@O){$LmPgMnrcDj0-N@@Cw6JsS;N
zN-pJEZ+;4z!J1#gAHy&L0ymK#Y~tPrN1t);O>iT|g%RWKy*B~gmyZ%YT+^wi(;J7Z
zcg4)Yp#k&heAs;X9X@hL<lr2=Uzj+8QKjs`m_WGZS#IGhh!#GOPTc@Qt>m(Eb3=-c
z20qN1fm4=~%bL-If`uo@!k1vtOH0k6FZ>hD4Tznk=F{%`J84ktK|rR><1{g3tP@q!
z+o;SB5{!|0>OzbW(9^vWk81<wu$)V#V|wbt<9g=hx4<?)&v!I=@Y}Ro+;$|XdzYLg
z5axXL=>eM<&98D^b9kWMfvMRGkq_qTo<rTKbE!qta=rt(x*tkD-#`mR2zu93@1cTk
z0FN8#7x}{)Yy=*2)1GBtCQjU}SthAkv`lTogyU1281y_$Pc4BwG2$#S$<+}Lfb#k9
zFhpBM2oEknh4|__jAF(&pu<5U^!Q3O*&ZX?Q<?4YWP9{%9lxjlHWMb7;R*Jm2!EUT
z7H&O)3SQv`lW229JQu8{+mC#rL&v34C&gH<FvpLXZ`0U2b=iRVcEN0v!E7X=#?w45
z2Z%Y|Z@!I~d3&jO@(ce?J{%V#LWTsqL{Huuk7K?WXeKpI2nC*jdUE<3r42mQ{Fs^y
zZzNhi&QGdwQ8*c5c#}Fhd98TT6l&g#A?#@$Q>u<xX3a++{dlo-QXiAT&Gd*-!TQaU
zVyK^V@1;czECAE62-~BWz(7PdL&B0MQSn_~;n^~Thx-$Vt)oN*pYl2SR7<Ye{8UW$
zF0M*5aF)@DIY3R(5XN1u=I2eApJ?Q&3AXQGv3olZBdUMJ;N(3a?^2MQpItr}qmVW;
z1;~s*j$w&3A5Q}`f`5u9_KX3%Ahh14cw}Maz#lokMy!00H;=)z6BNtYR9!YTKbtbL
zsmg51BPm}%5K@TBI5^cZV4f`KN<<uaXPD1j2m_0E@bt@Rb5H+cnC;;vT`sEN9o2o@
zI}no*J*1;m^+b>8rDubU<8<dmlu`s7=M%^nwkO2etNTsrYEN?|RgUwA&T|ILS@ghr
z8l8<I9xL2?U*+YIXOHe{xT)dBru1(cHa1`{;Nll&(i(qHHRT%`(jPRX-(9l+>t77;
zBOll^d;bg{?0zF@Ed0G%->*r(w<i7W%JiEJ*;z07`j2Y9YyUI#8XTqHvZY~j!<L3y
z8~jK2<@W>LI*OYl$&<S@n%o_mR!zxRA5NsYk_m%DE|+UsYb+V)G9uVd4J0G6Xu?=#
zTpYq}uqDs|RCJp$e{~?0NK|hKM?&E)qdT1JjD;3lZY09-2r+i|bcMUa(In?)7blGE
z*bFwXFC5%v&~7lvYo`SY(p|f18mfb=x|>x8SamNWG&pXl4jRE&v^N~@2uFkAEkJY|
z-GLqcyjFw@7=f;y&HxpAw-|rh5fAiqVuQIWyc}i}iU7iu><n9_8o@Y;gNa04ET0(A
zVx&9OZA4SuZQ-~PYqxqO88dpq@%C7}o7kl)s+mxryC)LuFxmsbWGqeuy+$}n`_cxC
zYJ;0%sU)`EQNx;|z1R^C6{r@7ccjRNxyroZ#)LZVggiR_X=cXbG3bl^eXayNVz4!l
zD4Cy%Cdm-;k`YOuXN_qL7x-O;Lm)&<+?hsuJk~8}c`c|IQ*o1ua_QS6$xfpu5Dy>_
z!f{I+t<8C@2a%Ilg0=*I^GCRKi4noegik^zr9q#@e5B<26wS3$G_Ne^358uvp~4M|
zLT_Lb)0ABs=uHH9z~XXYe3Kdyuw;pW*iQzmu`IM6yY5P!PF<(XwJ)W2khG^d{#rIW
zLiwp|R>u-Re?6Oh6!79VvRMQ0moUuhksm_78Tnt3e*yWX)7k8P<bS|(#Srq<c+@+I
zy!72{)`6!GCl(6xksqZs5%MZ5>u4wSAF;H(4|yeSqX&`ShkOh<CAF8X-M*e_J3QLl
z>81KUt<*#K^YGsS9Y0NkGEGj;*hQ!(dJk;lF96CttI8{{cb9F~?_$@TyY#XvFS>w4
z5dTX2XFx}ruqqMnX8e2Li}kRHHI#eqwKdLg+M6|~#43UM0{&OPM>WWE^(YeW!}xDS
z-L=SbbxGnB_4HiWjrtBu)wz0mo1G{Opibcb$G91#t!<nAc6%Rhp3|HGO%?i}y_jaR
zpdmdB)ZIe;gL1AN+hMMV0aQKe-$#AcT90;Y@8oq2f#yX#H2>7vj{OjCgxc5zd;)!b
zqg96O1}=c=zlr+qp#Dn*{JV4fPXq5oj3f){W4xaF<|Oc)h@CLWUs3MaYg<uX`DMFV
zZtQWaD4&0SiLbowUZ<~oaerxZd40TmaYK1sL;3s_<p$_LzoJ~{@4qfbjO~LQ8EPc$
zY&z?KvmQ9>fwLYs>w&W#IO~D49ysfPvmW^Wzys>snmVVZOr|9WC6!MXAoIw>Q0i>h
z)8~qcDp{}2#mVDmOsyp)S?+ust%oS7`gcCe#t42`E=o!WKsxRaxjM_mclCjiTT|?E
zsUTi`@kHxUO4c)H1k>`8Qhlw#bG7#2C!By$>u2ol^K_2P4@!i>_ee#Q-iu}X(k*;F
zOswZienp3^Ty7u8ISLz%95$pptUOgc=fski!*5GRV7bO&rH2&%|14AIz>5X0l7`(V
z^LCl{%KToLKP2;S%lvVfKO^%OW&Wzn-<7#jcH|tHU!pp;SdmW%t8<m=9Hlxpsm@8N
zbCIhX8<!aKH@0D?EoEF)b5%`U^};Js99+0SAdD;O7GATkZebm}9F2>B5#8A%qm6Hf
zOrn^ePk=?xKP5-c9yMx;cw(Rwajt89sGEK<F0Ud~CTqy1>7PTj%~txqD07t35f7VQ
z>Z43^CU{-R5Fopw1P6psSMnV}LZ`4CH1cc@-ioVZEfpW4Vu^kql{M$>$efpe&STeq
zjB@GMNYJz{6j|vuDwl0XQD;t*GRLC-po&rweDzHzveGlSIIf|x17b;rPBb>d9)OuP
zk7EvqZ4SpaLh6pVKz386m(UgoY-vR1{0#`^m&~~h<<bN$CFh6+PUOmt(i$pP{*@Xj
z{ZX#$;iIIaidZX4_Y(fxdq88rLfcNX<e+{jnL8Ubu6ej<I*Bg!wzhvZXs6wsBd`zN
z-#6_y9Hbii;9Y%XEhuQNxg0#uSN2(|opvb)5B8O9hFEP{H3uK<D<hM%vd@9WHdtTw
z7lLL{%l5(jePz!P>v?30eeh5pe0L9uCHg!tmvH?`bP`me6UlTs*sqoB#nokRsWXrR
z>BU7M>lKpAbrO@yzEvXG`xt)?M!TCb<c?#8q>w9)RXPh{`|oJic49H8?IO_ZA_(mw
z2<;*W?IQZ@HXeNr5q)+UeSHvNzdr&Iw!1N>+4oQp(kLcNbf^CJ1d8}8{T>Fx3CBC<
zfkpR2FUPxWsAc*WfpEM}HtTvl;m=Unqkk0w9UoA+Qs0UR)$t*f4gCY6$x?Z~>jB8b
z5pV+PT+fmqjRT8a?L=kcK)tJu$Q&H-xh^AFB^+pR9i&E_9N6Ug1Ci;PA9A+39wceg
zv@o(x*Ibe|U28)-J+93}=HWonbtl!#)Ot{}!xblBj<yQGU9K4<`W&qlz&)-<NLr<~
z4OzeIMiOwI_IUvJx$YsJi?z>#Y@h2V#8a*90r0TvNdjs(aKLqxfLab5bUjBv9S0tD
z{fdAqv}*7ia?#Io99MGPhFpK5nya+C(azJ%_Rpv)DcujUX$OevtEGQJ)0K_`WUu~t
zqWuOHnO;WKkC5fM{t}hH-3GLWIrSRgoI1%Uqpl?jrr(PKVQF}9?4ZJWf&B`A89WRO
zK2_MoRACoWg<VV)K+;r!sG_L?@k3L^j5dfZ<@1Ge&ZSUQr}@Hp4xyP&^Mxoobeb=m
zm5)J!uG4%W${s#VJFAGbQm6UCc`lD2L?oXtY((Ik`)lx(p9YcUQp1Grd81O>9|qoi
z5qzMzh}SmQ$vh|EDs2v_F(z*T*P={G+y5BR`w(G*z&6<B?#_{eTgaLtWNl^T15|nb
z0(lGJ`{i8DU1-w%8nt^qHDn)rvd_Ja(^8Fn@X<aupLWi_mxBlUDy||676`$<3Lmjv
zMy$5Ml`~(K^6vpmX(?(e=%GuiBd^&8TiDF6mW#SOM8^C?^J5ZoErphCa1--3yG7;K
zB~8WS)aVK!V=J7gQ~L-Lpvy0$y4HOjLCd%LQPlR=zzjE^U|Xq3)1t7YRoL=qU&W2&
zz>W2&MY9z%NswRI(pvEfsVyO&*IyyTe_Ld1R)lFOLAu8K6<K+Us=SS9J4Iy+P)>aj
zT5#%1_?7TZoo1)<F(9;ERq_`R?T8#o+WzTK%gyfvyQpXz^a(q5lN`+b!i?Qyy={<r
zO)20CDS&iK%jUf|$V$3zQmR-&Ht!MDtuybI)yqXKJB1czLH)gg0DbXuNpVn8Q0)Vv
zmhMlmWhNv#^-6xFE}DHQ+4c24x}P8o@hjQ&AV`E=Ul(?55|(^jS@K26@)OBImV8sR
z)#?2^S^3{oB^3FVsH6vs@-I>M{op$JEY4rW0Ils`fcv@o2(P48c^rHnYHB|qX8Yie
z``i^2eiB4jEGO_WG7Ayw?UfS0soI57KN9Wwya#0^!Ia25KM|E{0G6Leu6SzsY+lLu
zyB`#_{a+)ccyK;N!HMt_j(JKrrp3+g{YQlkKH;XL6u0Pp?-KgG9uhr)Q6e|}PBaQn
z%(0QJf25vxUlQ_nN+8ts)2lN#RcU{u4x@If$l5;&_STACQuriMU>5PZ&*2TQr6S`k
zy(&mUu<onW(i=BOqA%qF_%&*U2k;vzfW2r<dqcEY?`|f>Ka+L&Vfkkn1r@iF3hzl3
zX6}^Ie<)RW9<s^czetU}yCuN`O2<U1n`x&7|5oppWp#(DMi0F$csBv%4a!AdK-EL`
z2cjN*aL!vcDC6Kay5CDm{)P?I+Wx=Cb5JF>)1gU2D|yE`G|_YQm3;PeXwtT-zYt#u
z7l+uYIz!T(nzXg*HtLE}P4-L`_n9uit=_w<go;YjbD^nMio)H0?l^JV2O-jBSHBFK
zdks0$A&}~4VRN4*$`XN8zX_ZBc|tk`QvD!o?u(>g>Hpv|yVSB=tI#ek_m<hq=ayG!
zGfO@2u<asNR<;}{^r$SNpB!_Ty?ja8(y}Iu-4ZuQsNnQ!_VTha!fOCaXXJ<$V`-#c
z!-&XLqg_`em@lnr@=`nXl@G0`B93YJb)FP!L!Xo_DgdYF!RxB1o=v8&RWqeAly~Hn
z=lVihrFLnd+IpedE(w-2dC7`tN<P%nv?ed9J$)W5t3~CEOSp;URuc;hbQc=vu@V;f
zq9SiC**dcf#?CCI=Deba(E<`(^a7F0wn`Jtent3oj<rU$BkRt|)s;->k4mX7Y^MgQ
zq)HMzS2i!fxp{?U_wXehJ4qd-y<a{9SB+*8{A;zTM7%Z<4R)nM;adE<+JE&T`BUqT
zXsWg?f*-biJV|YGXHWR!2opOKc>m7NH~Bjfemc)WC!uPBH0HQ~nod#)ziCYclYTno
zgfl(BVQ3J566Umt0{rc%XpkSW0tHVln_4PXogCHI=-c&KwrX#`zD4&UPduB@FR)!U
z%QgpoSOk$a_FrI|+%x!HZEZ&|IMFj#DW`n`_k02p>A{aSk9#VS428n&!d1#^@DLP@
z#Dud}HLWrKPn@NmmCrqFcj?J*==1syA4t{dHLqYBKy$V`JO1<!KhV9pcX#j2wuN4O
zm#x;T-+cIhUMrx5UQU;E4gxn%Q2=Li>YI?+^%|;r3$BNpw_iD6-+aCOk|*>|eUAPQ
zsl$(?>N-12z4Ow;`sTr(>CZiByL{Flaog%ExZvxZ&k<<9-H@Pt`pLt)^;!CXCv6L6
z^#W-pdb@3bqThWuVl!t6CR+oJ>cNcb33S@GKOrEe-TTC@`wj!J^NZ~wvceASDbxQ*
z*CusQ-2Vmy%FZLp^_q<SnZveAs1ezEvuqt0o3l?++Feh4QU7PX_wbW?-CO#uALw=Z
z1GXhq`m=eytCiM|+UmTvOT7Bv@*g~quf9Z8-+S{yeYw8uu%&LjwaMMou>E%11=3z8
z+SYH@m+BYhg{{d6)0ZuOt5^SwzNY_YwhL$3YAlB2h<57_e8zU6_h#^G&hPW<(*7rx
z*=!gqbO4Wjqu!D`v&O>F&~_T9_|Z6J{^TUVjKR)-q%H`F>~{Rt92|*MB8<a|p#<2}
z%y46v&i%3Wu2>++Qqe@DBO1mDzs^9MKIdRGE?ii&5olC*CsWbzGMsnBG!q0EhzC2D
ziEl<2)g4Clx;mq}o&QX~x+fk>@&>9g;l!dHXjUD+V~NgKJZa>bq+PoFK{l4W-5*Nu
z-{e!V8{ef!h7;Svi6rdc$}YTOF@C&{U*VI1n8TRm7-pk7FUjyL`woL&3K~jcQcOf>
zK5Gsn<B=UK+3P3gFmZr#T{|5iO)wl_gfH<80Gb0BKWMlL#}Vly9&r#uz5y~oo~B&I
z4^v7D37-gmj`tTy?1kDpawqgKp?Ai*!?l5)SRxtginZ7JR<2xFds8gFEzuJQh7+}o
zwbcJzT{yv-3?#N0i)tHE$xel9h=x|;lLGNb8ypncuoGVw=uR}oQt?FC+!2lkBZ=^$
zS}~uqP$Z7d2#0P8xA9|zaAhjGEgIV%_4|b(p;(H4kzhqE)`i+wdpknK-?R?j9tec|
zJ;^x3mlw$C<P@sno!FMb2N+`86RcMRlcgt`qBuLsC4k7^7C<zw-PmldZ)$8<<6qIR
z!OVG_PJJeFVafQRMJaGqB%VmF#raiy<{=Pjq}b$#XQMHW;I^#7`A_N_8X+w)_&f|E
zr8sab8s1KwN_RsU^N6EB2Yp*%ou8B=B7ZOy_XlHLsqQHA$GSq~fkZgTc)Xy8!ZKW<
z;T_4Ac&shZ7U_y4cjkHk0mb+sS7CRz2(FL3o#F!(PymBy8#LzUYd0pKC&zB5I2LYC
zkduXG%Fi3<;}jTSh|BJ^6p|z;5$FvIKO+!&8sc$$&qAHUZHk7%JD8u&H^YHs{1wdK
z9Sen}pI1iO+wqACoMUD;n(J4fht_RuZQ0oBZ&<$?owl~sq7?TQ+9P2yM(CDk8N$K^
zD;t_yn%1uNuWD#)UANvM3Eo8G+-yG3_$gupE=Eh^y0vJ3wYj#@v@rB0=~S}nkw_H7
zSA3tNJFqR>6px2PsbGYnx4B^xzL#OHZ)j~=w{}Bbf^viG{wM;<VlfR$?jOQ4KlQK}
z0aq{L(<KH}gFNdjVulY19??0!bEn{;E3^=-$lzYet`ylKkzK`to#EiNrf4wEPoFak
zLd-oN=2L1Q8KzL@riEII=c%4}xHl3@C5l6tVi?~=N%e@nr`dg7Dj5oQ#KROLf=Sqy
zo5n1?$;}YYbyHmF>miHr2vnTxc(PO9L$cP<H%O9`(mQhRHw)=Y#IKdCNR;i&nS1|f
zt+KOHrrrw{;&Emc3-aP*$5*b@drNDb9dA+9`^`dpS?)b=A$~gc&aD_xwzC<mZj$$Q
zc2=JI{+5-?&fHABH!j3`*syw0Y^kv0bhi4Alm!j{%w)soTQQ<+XWrcRS**2oyyaEj
zXDP(bX1;2}isEHEn}hFVO@==w_a52GWT&^k>II9X!p^E#-K6+Qj0K;asb#uF2*;mS
zh>(z-{mm5kNu0ywd};zcE!Q0EFKW*%S3eyk4R6zn{X-iWHuP^#QNNSsgnzUMUnB9w
z=L}YI`Fa149+WoVUz?=A==~(g&z)N!Vx)WVuT9eL-vWLzyPlW=Kb0T;gUg@E{&pS;
zmh>l-C-V=HOxEsOoYy-c4jpJ%n<U<5aEDBcXM^J>j*Ba~oQd(-1pG|+ZAgZknx1Z!
z_~I>W`dtIj7jFUI54^!`D7?H3fyN~9e~`;5q$TJm$Cppo-=w#za=W^K>6Iz)e+Axy
zb|;K4=rfs|4P4H|`DfM?^mS9<S5AT7#PQ{<_+`dzQ_$~{^u_Oo{%#8TZ%g{(P4dSj
zzIdPd$P{u$fG58d&sT2%UxE3mUd~tYoeX<-3OUoTzd>=L*7!8eDa{7n;HARLMV!76
zN6=#6E3oyX&Tl9Ra|-%RL=Owp{Q^6mJkfhG(x*knd067NbNVW#UX?4m=$-0haq<wS
zuViW$7F(e_(FEy%e>#ii^H;fl^7p$>OZx9|IUc5VUX`9tNcn0*SIPgi#H;OFg&zf;
z+AWUrcYr5*)rOQhlJcI!7r!rbK>?yK-U*)xJoWolX_snuE~l5qSh5=64GXM_E4Z9W
zW)#W!>=bgkh#vi3Y}bw{<kM?>(!cnP#SekMKq+C#Kgazr$@5Gwo=hZD?d{kx@cXZC
zT<>4gw4oJi!Tg8n{K;-Veh^7ZGrvC+^LKP%(bMG*Vb37p52SWrNkYG(#G<n1nyatI
zckU)5qEG2X`~g~vV3n1O<Kug@3*Zl>y1VIXdltw~&-DeQ6BkdkzcFEl)3QNH#;jW3
z&}{n6wJTwbFin*U8o-Kg<;`mwnwuIyWvSu&CqQ8f#BchfG~deg%)ffgI=(PkwQ7Ue
z>ThjWvBsnp_@@Dbi4-^WlRU^5FIqnO1Azr!Fs!dvA0E7B!UCB-l9;cjZxCW{0G~uG
z#QDRaKr$ddfml>k^nt~E?IQWD#XN?@iw`dr;{u_=MY;9!iiNaAYa#y0#zKnXwKsK>
zT=`A-Bx4~bf)Qql4=d((&fpMy&asfvy6!Ig8Y6dKS0OdzCSv~1Ks1Earudv>o&X<%
z<m<sCew<2sN(ixhbdc`W%{lS)$viW8bkfyt(Wc2%OWdiJ?f$kzLi&O4nPD5wx}2ru
zB5jQ1_o?_NHw*Qe?7Ngxt?ucwokiS!eDhKb$SJ>O`BD2Sg<s_?l+J1rJG-$zfE?k@
z^G;RbdqF*{CW@`dnuZlk)yY7Igz)yIrY#l0hi@YxRzpx{AkoQcLOY|7Cvxn0@fvK#
zCTM%d3iweM4|fHKL6&;DlB|Yza}A0$9Wj7O+DoXRNf1OaKJ;qB`29DZ1Uo|zq9DO4
z=0Cxq%C`kUhGqiYkst~*d~49DtOgS>|NV7MDBPCn0C^xvTbB~VZ<J%I9B6C9Mu39w
zso$_F!&R2s#s77Y);)N%#M(q8wO^p};*}fGTGNF9aG)T^ul5yGZdi#?#_g^kEd*rc
zPfA5peo+DSME;ei$4^t~Bq~o2%Yw?4{&*M7lX^aoAV0`}_o6xd4@rKNdn8QT6J)#M
z#}g<|biM(%>Kwm1_pNeu&Rb9m8LCYAt@x)i)JlJ~Kce!>C9%?9(W!hJ>gi0jg4I5X
z%IOpkrDFXBmN8nFQlBV(wcn!h#j?Xy`^EO}lKd;AK(#NUa&=Cf_~|%AG5?o=qbAh$
zk2<HXa(RRT6e5wX)|3Wd<yZSWD(96J2KD@`@&}9g<@#Uk3n_h-2nDO0J{V<fzh3gI
z{2tYi6i@`Re=Oox`%5Y>-nJ^X|2vXj^`BZ-sGNW4RJKi-(p#Y(LxJ>H{A#~S<>Vi0
zDsKPhsI>B{eJz#uD8W)lar@5#kHZgq3=xU?9ftZ&d@;R(RsKK0Lr3;#3{q12Wc75S
z4@uDpC=V<AA5cO1D}J?K<og(YCGS-fD0UUU+ILgmt5@HpSMjadRdMhzc&+WL-+9#6
ziWY`Mmb)r`mA{D^D}QcXTqjs$Hi2Iu-a&z6DgD)YrcN?0R)jKFdMosMC{X*w?I?cf
zMn%qP2_POyb`(@uV8nHc6u8xjM7fxD<`n!tmHctq2(+YP{`05cuU#gX))mEXF+a`k
z#W>=;^E$!5zG%P{)6)D(xMKeeHwylk3LF9E+EQ{TFH(_H)?Y6e#^eMt58xvFtGH43
lFRE63FF3@t4#&DEJ@+yCSAb?R{`WR<c1FjXt*HoP{~O6ZL2v*7

diff --git a/serial/serial_gs_pagerank.c b/serial/serial_gs_pagerank.c
index 8f15730..5b056b9 100644
--- a/serial/serial_gs_pagerank.c
+++ b/serial/serial_gs_pagerank.c
@@ -14,7 +14,7 @@ int main(int argc, char **argv) {
 
 	initialize(&directedWebGraph, &transitionMatrix, &pagerankVector, &parameters);
 
-	//Starts wall-clock timer
+	// Starts wall-clock timer
 	gettimeofday (&startwtime, NULL);
 
 	int iterations = pagerank(&transitionMatrix, &pagerankVector, parameters);
@@ -23,7 +23,7 @@ int main(int argc, char **argv) {
 			\nTotal iterations = %d\n", iterations);
 	}
 
-	//Stops wall-clock timer
+	// Stops wall-clock timer
 	gettimeofday (&endwtime, NULL);
 	double seq_time = (double)((endwtime.tv_usec - startwtime.tv_usec)/1.0e6 +
 		endwtime.tv_sec - startwtime.tv_sec);
diff --git a/serial/serial_gs_pagerank_functions.c b/serial/serial_gs_pagerank_functions.c
index 380a3a1..46cbd7f 100644
--- a/serial/serial_gs_pagerank_functions.c
+++ b/serial/serial_gs_pagerank_functions.c
@@ -1,32 +1,210 @@
 #include "serial_gs_pagerank_functions.h"
 
-const char *CONVERGENCE_ARGUMENT = "-c";
-const char *MAX_ITERATIONS_ARGUMENT = "-m";
-const char *DAMPING_FACTOR_ARGUMENT = "-a";
-const char *VERBAL_OUTPUT_ARGUMENT = "-v";
+const char *ARGUMENT_CONVERGENCE_TOLERANCE = "-c";
+const char *ARGUMENT_MAX_ITERATIONS = "-m";
+const char *ARGUMENT_DAMPING_FACTOR = "-a";
+const char *ARGUMENT_VERBAL_OUTPUT = "-v";
+const char *ARGUMENT_OUTPUT_HISTORY = "-h";
+const char *ARGUMENT_OUTPUT_FILENAME = "-o";
+
 const int NUMERICAL_BASE = 10;
+char *DEFAULT_OUTPUT_FILENAME = "pagerank_output";
 
-void validUsage(char *programName) {
-	printf("%s [-c convergence] [-m max_iterations] [-a alpha] [-v] <graph_file>\
-		\n-c convergence\
-		\n\tthe convergence criterion\
-		\n-m max_iterations\
-		\n\tmaximum number of iterations to perform\
-		\n-a alpha\
-		\n\tthe damping factor\
-		\n-v enable verbal output\
-		\n", programName);
-	exit(EXIT_FAILURE);
+// ==================== PAGERANK ====================
+
+int pagerank(double ***transitionMatrix, double **pagerankVector, Parameters parameters) {
+	int iterations = 0;
+	double delta,
+	*vectorDifference = (double *) malloc(parameters.numberOfPages * sizeof(double)),
+	*previousPagerankVector = (double *) malloc(parameters.numberOfPages * sizeof(double));
+
+	if (parameters.verbose) {
+		printf("\n----- Starting iterations -----\n");
+	}
+
+	do {
+		memcpy(previousPagerankVector, *pagerankVector, parameters.numberOfPages * sizeof(double));
+
+		matrixVectorMultiplication(transitionMatrix, previousPagerankVector,
+			pagerankVector, parameters.numberOfPages, parameters.dampingFactor);
+
+		if (parameters.history) {
+			savePagerankToFile(parameters.outputFilename, iterations != 0,
+				*pagerankVector, parameters.numberOfPages);
+		}
+
+		for (int i=0; i<parameters.numberOfPages; ++i) {
+			vectorDifference[i] = (*pagerankVector)[i] - previousPagerankVector[i];
+		}
+		delta = vectorNorm(vectorDifference, parameters.numberOfPages);
+
+		++iterations;
+		printf("Iteration %d: delta = %f\n", iterations, delta);
+	} while (delta > parameters.convergenceCriterion &&
+		(parameters.maxIterations == 0 || iterations < parameters.maxIterations));
+
+	if (!parameters.history) {
+		savePagerankToFile(parameters.outputFilename, false, *pagerankVector,
+			parameters.numberOfPages);
+	}
+
+	return iterations;
 }
 
-int checkIncrement(int previousIndex, int maxIndex, char *programName) {
-	if (previousIndex == maxIndex) {
-		validUsage(programName);
-		exit(EXIT_FAILURE);
+// ==================== INITIALIZATION ====================
+
+/*
+ * initialize allocates required memory for arrays, reads the web graph from the
+ * from the file and creates the initial transition probability distribution
+ * matrix.
+*/
+void initialize(int ***directedWebGraph, double ***transitionMatrix,
+	double **pagerankVector, Parameters *parameters) {
+
+	// Reads web graph from file
+	if ((*parameters).verbose) {
+		printf("----- Reading graph from file -----\n");
 	}
-	return ++previousIndex;
+	readGraphFromFile(directedWebGraph, parameters);
+
+	// Outputs the algorithm parameters to the console
+	if ((*parameters).verbose) {
+		printf("\n----- Running with parameters -----\
+			\nNumber of pages: %d", (*parameters).numberOfPages);
+		if (!(*parameters).maxIterations) {
+			printf("\nMaximum number of iterations: inf");
+		} else {
+			printf("\nMaximum number of iterations: %d", (*parameters).maxIterations);
+		}
+		printf("\nConvergence criterion: %f\
+			\nDamping factor: %f\
+			\nGraph filename: %s\n", (*parameters).convergenceCriterion,
+			(*parameters).dampingFactor, (*parameters).graphFilename);
+	}
+
+	// Allocates memory for the pagerank vector
+	(*pagerankVector) = (double *) malloc((*parameters).numberOfPages * sizeof(double));
+	double webUniformProbability = 1. / (*parameters).numberOfPages;
+	for (int i=0; i<(*parameters).numberOfPages; ++i) {
+		(*pagerankVector)[i] = webUniformProbability;
+	}
+
+	// Generates the initial transition matrix (matrix P).
+	generateNormalizedTransitionMatrix(transitionMatrix, *directedWebGraph, *parameters);
+	// Transposes the transition matrix (P^T).
+	transposeMatrix(transitionMatrix, (*parameters).numberOfPages, (*parameters).numberOfPages);
+}
+
+/*
+ * generateNormalizedTransitionMatrix generates the normalized transition matrix
+ * from the graph data (matrix P').
+*/
+void generateNormalizedTransitionMatrix(double ***transitionMatrix,
+	int **directedWebGraph, Parameters parameters) {
+	// Allocates memory for the transitionMatrix rows
+	(*transitionMatrix) = (double **) malloc(parameters.numberOfPages * sizeof(double *));
+
+	for (int i=0; i<parameters.numberOfPages; ++i) {
+		// Allocates memory for this row's columns
+		(*transitionMatrix)[i] = (double *) malloc(parameters.numberOfPages * sizeof(double));
+
+		// Calculates the outdegree of this page
+		int pageOutdegree = 0;
+		for (int j=0; j<parameters.numberOfPages; ++j) {
+			pageOutdegree += directedWebGraph[i][j];
+		}
+
+		// Populates this row of the transition matrix
+		if (pageOutdegree != 0) {
+			// Calculates the uniform probability once.
+			double pageUniformProbability = 1. / pageOutdegree;
+
+			for (int j=0; j<parameters.numberOfPages; ++j) {
+				if (directedWebGraph[i][j] == 1){
+					(*transitionMatrix)[i][j] = pageUniformProbability;
+				} else {
+					(*transitionMatrix)[i][j] = 0;
+				}
+			}
+		} else {
+			for (int j=0; j<parameters.numberOfPages; ++j) {
+				(*transitionMatrix)[i][j] = 0;
+			}
+		}
+	}
+}
+
+// ==================== MATH UTILS ====================
+
+/*
+ * matrixVectorMultiplication calculates the product of the multiplication
+ * between a matrix and the a vector in a cheap way.
+*/
+void matrixVectorMultiplication(double ***matrix, double *vector,
+	double **product, int vectorSize, double dampingFactor) {
+	double webUniformProbability = 1. / vectorSize;
+
+	for (int i=0; i<vectorSize; ++i) {
+		double sum = 0;
+
+		for (int j=0; j<vectorSize; ++j) {
+			sum += (*matrix)[i][j] * vector[j];
+		}
+		(*product)[i] = dampingFactor * sum;
+	}
+
+	double normDifference = vectorNorm(vector, vectorSize) -
+	vectorNorm((*product), vectorSize);
+
+	for (int i=0; i<vectorSize; ++i) {
+		(*product)[i] += normDifference * webUniformProbability;
+	}
+}
+
+/*
+ * vectorNorm calculates the first norm of a vector.
+*/
+double vectorNorm(double *vector, int vectorSize) {
+	double norm = 0.;
+
+	for (int i=0; i<vectorSize; ++i) {
+		norm += fabs(vector[i]);
+	}
+
+	return norm;
+}
+
+/*
+ * transposeMatrix transposes the matrix passed (by reference) in the arguments.
+*/
+void transposeMatrix(double ***matrix, int rows, int columns) {
+	// Transposes the matrix
+	// Rows become columns and vice versa
+
+	double **tempArray = (double **) malloc(columns * sizeof(double *));
+	for (int i=0; i<columns; ++i) {
+		tempArray[i] = malloc(rows * sizeof(double));
+
+		for (int j=0; j<rows; ++j) {
+			tempArray[i][j] = (*matrix)[j][i];
+		}
+	}
+
+	// TODO free memory
+
+	//double **pointerToFreeMemoryLater = *matrix;
+	*matrix = tempArray;
+	/*for (int i=0; i<rows; ++i) {
+		free(pointerToFreeMemoryLater[i]);
+	}
+	free(pointerToFreeMemoryLater);*/
 }
 
+// ==================== PROGRAM INPUT AND OUTPUT UTILS ====================
+
+/*
+ * parseArguments parses the command line arguments given by the user.
+*/
 void parseArguments(int argumentCount, char **argumentVector, Parameters *parameters) {
 	if (argumentCount < 2 || argumentCount > 10) {
 		validUsage(argumentVector[0]);
@@ -37,12 +215,14 @@ void parseArguments(int argumentCount, char **argumentVector, Parameters *parame
 	(*parameters).convergenceCriterion = 1;
 	(*parameters).dampingFactor = 0.85;
 	(*parameters).verbose = false;
+	(*parameters).history = false;
+	(*parameters).outputFilename = DEFAULT_OUTPUT_FILENAME;
 
 	char *endPointer;
 	int argumentIndex = 1;
 
 	while (argumentIndex < argumentCount) {
-		if (!strcmp(argumentVector[argumentIndex], CONVERGENCE_ARGUMENT)) {
+		if (!strcmp(argumentVector[argumentIndex], ARGUMENT_CONVERGENCE_TOLERANCE)) {
 			argumentIndex = checkIncrement(argumentIndex, argumentCount, argumentVector[0]);
 
 			double convergenceInput = strtod(argumentVector[argumentIndex], &endPointer);
@@ -51,7 +231,7 @@ void parseArguments(int argumentCount, char **argumentVector, Parameters *parame
 				exit(EXIT_FAILURE);
 			}
 			(*parameters).convergenceCriterion = convergenceInput;
-		} else if (!strcmp(argumentVector[argumentIndex], MAX_ITERATIONS_ARGUMENT)) {
+		} else if (!strcmp(argumentVector[argumentIndex], ARGUMENT_MAX_ITERATIONS)) {
 			argumentIndex = checkIncrement(argumentIndex, argumentCount, argumentVector[0]);
 
 			size_t iterationsInput = strtol(argumentVector[argumentIndex], &endPointer, NUMERICAL_BASE);
@@ -60,7 +240,7 @@ void parseArguments(int argumentCount, char **argumentVector, Parameters *parame
 				exit(EXIT_FAILURE);
 			}
 			(*parameters).maxIterations = iterationsInput;
-		} else if (!strcmp(argumentVector[argumentIndex], DAMPING_FACTOR_ARGUMENT)) {
+		} else if (!strcmp(argumentVector[argumentIndex], ARGUMENT_DAMPING_FACTOR)) {
 			argumentIndex = checkIncrement(argumentIndex, argumentCount, argumentVector[0]);
 
 			double alphaInput = strtod(argumentVector[argumentIndex], &endPointer);
@@ -69,8 +249,18 @@ void parseArguments(int argumentCount, char **argumentVector, Parameters *parame
 				exit(EXIT_FAILURE);
 			}
 			(*parameters).dampingFactor = alphaInput;
-		} else if (!strcmp(argumentVector[argumentIndex], VERBAL_OUTPUT_ARGUMENT)) {
+		} else if (!strcmp(argumentVector[argumentIndex], ARGUMENT_VERBAL_OUTPUT)) {
 			(*parameters).verbose = true;
+		} else if (!strcmp(argumentVector[argumentIndex], ARGUMENT_OUTPUT_HISTORY)) {
+			(*parameters).history = true;
+		} else if (!strcmp(argumentVector[argumentIndex], ARGUMENT_OUTPUT_FILENAME)) {
+			argumentIndex = checkIncrement(argumentIndex, argumentCount, argumentVector[0]);
+
+			if (fopen(argumentVector[argumentIndex], "w") == NULL) {
+				printf("Invalid output filename. Reverting to default.\n");
+				continue;
+			}
+			(*parameters).outputFilename = argumentVector[argumentIndex];
 		} else if (argumentIndex == argumentCount - 1) {
 			(*parameters).graphFilename = argumentVector[argumentIndex];
 		} else {
@@ -81,6 +271,10 @@ void parseArguments(int argumentCount, char **argumentVector, Parameters *parame
 	}
 }
 
+/*
+ * readGraphFromFile loads the file supplied in the command line arguments to an
+ * array (directedWebGraph) that represents the graph.
+*/
 void readGraphFromFile(int ***directedWebGraph, Parameters *parameters) {
 	FILE *graphFile;
 
@@ -103,11 +297,11 @@ void readGraphFromFile(int ***directedWebGraph, Parameters *parameters) {
 	}
 
 	if ((*parameters).verbose) {
-		printf("Line count of file is %d \n", numberOfLines);
+		printf("Line count of file is %d \n", numberOfLines + 1);
 	}
 
 	// Each line of the file represents one page of the graph
-	(*parameters).numberOfPages = numberOfLines;
+	(*parameters).numberOfPages = numberOfLines + 1;
 	rewind(graphFile);
 
 	// Allocates memory and loads values into directedWebGraph (matrix A)
@@ -122,165 +316,62 @@ void readGraphFromFile(int ***directedWebGraph, Parameters *parameters) {
 			if (!fscanf(graphFile, "%d ", &(*directedWebGraph)[i][j])) {
 				break;
 			}
-			//printf("directedWebGraph[%d][%d] = %d", i , j, (*directedWebGraph)[i][j]);
 		}
 	}
 
 	fclose(graphFile);
 }
 
-void generateNormalizedTransitionMatrix(double ***transitionMatrix,
-	int **directedWebGraph, Parameters parameters) {
-	// Allocates memory for the transitionMatrix rows
-	(*transitionMatrix) = (double **) malloc(parameters.numberOfPages * sizeof(double *));
-
-	for (int i=0; i<parameters.numberOfPages; ++i) {
-		// Allocates memory for this row's columns
-		(*transitionMatrix)[i] = (double *) malloc(parameters.numberOfPages * sizeof(double));
-
-		int pageOutdegree = 0;
-		//Calculates the outdegree of this page
-		for (int j=0; j<parameters.numberOfPages; ++j) {
-			pageOutdegree += directedWebGraph[i][j];
-		}
-		for (int j=0; j<parameters.numberOfPages; ++j) {
-			if (pageOutdegree == 0) {
-				// Introduces random jumps from dangling nodes (P' = P + D)
-				// This makes sure that there are no pages with zero outdegree.
-				(*transitionMatrix)[i][j] = 1. / parameters.numberOfPages;
-			} else {
-				(*transitionMatrix)[i][j] = 1. / pageOutdegree;
-			}
-		}
-	}
+/*
+ * validUsage outputs a message to the console that informs the user of the
+ * correct (valid) way to use the program.
+*/
+void validUsage(char *programName) {
+	printf("%s [-c convergence_criterion] [-m max_iterations] [-a alpha] [-v] [-h] [-o output_filename] <graph_file>\
+		\n-c convergence_criterion\
+		\n\tthe convergence tolerance criterion\
+		\n-m max_iterations\
+		\n\tmaximum number of iterations to perform\
+		\n-a alpha\
+		\n\tthe damping factor\
+		\n-v enable verbal output\
+		\n-h enable history output to file\
+		\n-o output_filename\
+		\n\tfilename and path for the output\
+		\n", programName);
+	exit(EXIT_FAILURE);
 }
 
-void makeIrreducible(double ***transitionMatrix, Parameters parameters) {
-	// Manipulates the values of transitionMatrix to make it irreducible. A
-	// uniform probability (1/number_of_pages) and no personalization are used
-	// here.
-
-	// Introduces teleportation (P'' = cP' + (1 - c)E)
-	for (int i=0; i<parameters.numberOfPages; ++i) {
-		for (int j=0; j<parameters.numberOfPages; ++j) {
-			(*transitionMatrix)[i][j] =
-			parameters.dampingFactor *(*transitionMatrix)[i][j] +
-			(1 - parameters.dampingFactor) / parameters.numberOfPages;
-		}
+/*
+ * checkIncrement is a helper function for parseArguments function.
+*/
+int checkIncrement(int previousIndex, int maxIndex, char *programName) {
+	if (previousIndex == maxIndex) {
+		validUsage(programName);
+		exit(EXIT_FAILURE);
 	}
+	return ++previousIndex;
 }
 
-void transposeMatrix(double ***matrix, int rows, int columns) {
-	// Transposes the matrix
-	// Rows become columns and vice versa
-
-	double **tempArray = (double **) malloc(rows * sizeof(double *));
-	for (int i=0; i<rows; ++i) {
-		tempArray[i] = malloc(columns * sizeof(double));
+void savePagerankToFile(char *filename, bool append, double *pagerankVector,
+	int vectorSize) {
+	FILE *outputFile;
 
-		for (int j=0; j<columns; ++j) {
-			tempArray[i][j] = (*matrix)[j][i];
-		}
+	if (append) {
+		outputFile = fopen(filename, "a");
+	} else {
+		outputFile = fopen(filename, "w");
 	}
 
-	//double **pointerToFreeMemoryLater = *matrix;
-	matrix = &tempArray;
-	/*for (int i=0; i<rows; ++i) {
-		free(pointerToFreeMemoryLater[i]);
+	if (outputFile == NULL) {
+		printf("Error while opening the output file.\n");
+		return;
 	}
-	free(pointerToFreeMemoryLater);*/
-}
-
-void initialize(int ***directedWebGraph, double ***transitionMatrix,
-	double **pagerankVector, Parameters *parameters) {
-
-	if ((*parameters).verbose) {
-		printf("----- Reading graph from file -----\n");
-	}
-	readGraphFromFile(directedWebGraph, parameters);
-
-	if ((*parameters).verbose) {
-		printf("\n----- Running with parameters -----\
-			\nNumber of pages: %d", (*parameters).numberOfPages);
-		if (!(*parameters).maxIterations) {
-			printf("\nMaximum number of iterations: inf");
-		} else {
-			printf("\nMaximum number of iterations: %d", (*parameters).maxIterations);
-		}
-		printf("\nConvergence criterion: %f\
-			\nDamping factor: %f\
-			\nGraph filename: %s\n", (*parameters).convergenceCriterion,
-			(*parameters).dampingFactor, (*parameters).graphFilename);
-	}
-
-	// Allocates memory for the pagerank vector
-	(*pagerankVector) = (double *) malloc((*parameters).numberOfPages * sizeof(double));
-	for (int i=0; i<(*parameters).numberOfPages; ++i) {
-		(*pagerankVector)[i] = 1. / (*parameters).numberOfPages;
-	}
-
-	generateNormalizedTransitionMatrix(transitionMatrix, *directedWebGraph, *parameters);
-	makeIrreducible(transitionMatrix, *parameters);
-	transposeMatrix(transitionMatrix, (*parameters).numberOfPages, (*parameters).numberOfPages);
-}
-
-double vectorFirstNorm(double *vector, int vectorSize) {
-	double norm = 0;
 
 	for (int i=0; i<vectorSize; ++i) {
-		norm += vector[i];
-	}
-
-	return norm;
-}
-
-void nextProbabilityDistribution(double ***transitionMatrix, double *previousPagerankVector,
-	double **newPagerankVector, Parameters parameters) {
-
-	transposeMatrix(transitionMatrix, parameters.numberOfPages, parameters.numberOfPages);
-	for (int i=0; i<parameters.numberOfPages; ++i) {
-		double sum = 0;
-
-		for (int j=0; j<parameters.numberOfPages; ++j) {
-			sum += (*transitionMatrix)[i][j] * previousPagerankVector[j];
-		}
-		(*newPagerankVector)[i] = parameters.dampingFactor * sum;
-	}
-
-	double normDifference = vectorFirstNorm(previousPagerankVector, parameters.numberOfPages) -
-	vectorFirstNorm((*newPagerankVector), parameters.numberOfPages);
-
-	for (int i=0; i<parameters.numberOfPages; ++i) {
-		(*newPagerankVector)[i] += normDifference / parameters.numberOfPages;
+		fprintf(outputFile, "%f ", pagerankVector[i]);
 	}
+	fprintf(outputFile, "\n");
 
-	transposeMatrix(transitionMatrix, parameters.numberOfPages, parameters.numberOfPages);
-}
-
-int pagerank(double ***transitionMatrix, double **pagerankVector, Parameters parameters) {
-	int iterations = 0;
-	double delta,
-	*vectorDifference = (double *) malloc(parameters.numberOfPages * sizeof(double)),
-	*previousPagerankVector = (double *) malloc(parameters.numberOfPages * sizeof(double));
-
-	if (parameters.verbose) {
-		printf("\n----- Starting iterations -----\n");
-	}
-
-	do {
-		memcpy(previousPagerankVector, *pagerankVector, parameters.numberOfPages * sizeof(double));
-
-		nextProbabilityDistribution(transitionMatrix, previousPagerankVector, pagerankVector, parameters);
-
-		for (int i=0; i<parameters.numberOfPages; ++i) {
-			vectorDifference[i] = (*pagerankVector)[i] - previousPagerankVector[i];
-		}
-		delta = vectorFirstNorm(vectorDifference, parameters.numberOfPages);
-
-		++iterations;
-		printf("Iteration %d: delta = %f\n", iterations, delta);
-	} while (delta > parameters.convergenceCriterion &&
-		(parameters.maxIterations != 0 || iterations < parameters.maxIterations));
-
-	return iterations;
+	fclose(outputFile);
 }
\ No newline at end of file
diff --git a/serial/serial_gs_pagerank_functions.h b/serial/serial_gs_pagerank_functions.h
index 9eea0e1..0166f33 100644
--- a/serial/serial_gs_pagerank_functions.h
+++ b/serial/serial_gs_pagerank_functions.h
@@ -5,25 +5,30 @@
 #include <stdio.h>
 #include <stdlib.h>
 #include <string.h>
+#include <math.h>
 
 /*
  * Constant strings that store the command line options available.
 */
-extern const char *CONVERGENCE_ARGUMENT;
-extern const char *MAX_ITERATIONS_ARGUMENT;
-extern const char *DAMPING_FACTOR_ARGUMENT;
-extern const char *VERBAL_OUTPUT_ARGUMENT;
+extern const char *ARGUMENT_CONVERGENCE_TOLERANCE;
+extern const char *ARGUMENT_MAX_ITERATIONS;
+extern const char *ARGUMENT_DAMPING_FACTOR;
+extern const char *ARGUMENT_VERBAL_OUTPUT;
+extern const char *ARGUMENT_OUTPUT_HISTORY;
+extern const char *ARGUMENT_OUTPUT_FILENAME;
 
 // This is the numerical base used when parsing the numerical command line
 // arguments.
 extern const int NUMERICAL_BASE;
+// Default filename used for the output.
+extern char *DEFAULT_OUTPUT_FILENAME;
 
-// Declares a data structure to conveniently hold the algorithm's parameters
+// Declares a data structure to conveniently hold the algorithm's parameters.
 typedef struct parameters {
 	int numberOfPages, maxIterations;
 	double convergenceCriterion, dampingFactor;
-	bool verbose;
-	char* graphFilename;
+	bool verbose, history;
+	char *outputFilename, *graphFilename;
 } Parameters;
 
 // Function validUsage outputs the correct way to use the program with command
@@ -37,19 +42,21 @@ int checkIncrement(int previousIndex, int maxIndex, char *programName);
 // Function parseArguments parses command line arguments.
 void parseArguments(int argumentCount, char **argumentVector, Parameters *parameters);
 
-// Function readGraphFromFile loads the graph stored in the file provided in the
-// command line arguments to the array directedWebGraph.
+// Function readGraphFromFile loads adjacency matrix, that represents the web
+// graph, stored in the file provided in the command line arguments to the array
+// directedWebGraph.
 void readGraphFromFile(int ***directedWebGraph, Parameters *parameters);
 
-// Function generateNormalizedTransitionMatrix generates the normalized transition
-// matrix from the graph data.
+// Function savePagerankToFile appends or overwrites the pagerank vector
+// "pagerankVector" to the file with the filename supplied in the arguments 
+void savePagerankToFile(char *filename, bool append, double *pagerankVector,
+	int vectorSize);
+
+// Function generateNormalizedTransitionMatrix generates the normalized
+// transition matrix from the web graph data.
 void generateNormalizedTransitionMatrix(double ***transitionMatrix,
 	int **directedWebGraph, Parameters parameters);
 
-// Function makeIrreducible introduces teleportation to the transition matrix,
-// making it irreducible.
-void makeIrreducible(double ***transitionMatrix, Parameters parameters);
-
 // Function transposeMatrix transposes a matrix.
 void transposeMatrix(double ***matrix, int rows, int columns);
 
@@ -62,14 +69,17 @@ void initialize(
 	Parameters *parameters
 	);
 
-// Function vectorFirstNorm calculates the first norm of a vector.
-double vectorFirstNorm(double *vector, int vectorSize);
+// Function vectorNorm calculates the first norm of a vector.
+double vectorNorm(double *vector, int vectorSize);
 
-// Function nextProbabilityDistribution calculates the product of the transition
-// matrix and the pagerank vector.
-void nextProbabilityDistribution(double ***transitionMatrix, double *previousPagerankVector,
-	double **newPagerankVector, Parameters parameters);
+// Function matrixVectorMultiplication calculates the product of the
+// multiplication between a matrix and the a vector.
+void matrixVectorMultiplication(double ***matrix, double *vector,
+	double **product, int vectorSize, double dampingFactor);
 
+// Function pagerank iteratively calculates the pagerank of each page until
+// either the convergence criterion is met or the maximum number of iterations
+// is reached.
 int pagerank(double ***transitionMatrix, double **pagerankVector, Parameters parameters);
 
 #endif	// SERIAL_GS_PAGERANK_FUNCTIONS_H
\ No newline at end of file