1 files changed, 83 insertions, 0 deletions

diff --git a/unit_trace/gedf_inversion_stat_printer.py b/unit_trace/gedf_inversion_stat_printer.py
new file mode 100644
index 0000000..c314389
--- /dev/null
+++ b/unit_trace/gedf_inversion_stat_printer.py
@@ -0,0 +1,83 @@
+###############################################################################
+# Description
+###############################################################################
+# Compute and print G-EDF inversion statistics
+
+
+###############################################################################
+# Public Functions
+###############################################################################
+
+def gedf_inversion_stat_printer(stream,num):
+
+    # State
+    min_inversion = -1
+    max_inversion = -1
+    sum_inversions = 0
+    num_inversions = 0
+    longest_inversions = []
+
+    # Iterate over records, updating state
+    for record in stream:
+        if record.type_name == 'inversion_end':
+            length = record.job.inversion_end - record.job.inversion_start
+            if length > 0:
+                num_inversions += 1
+                if length > max_inversion:
+                    max_inversion = length
+                if length < min_inversion or min_inversion == -1:
+                    min_inversion = length
+                sum_inversions += length
+                if len(longest_inversions) == num:
+                    si = longest_inversions[0]
+                    if length > (si.job.inversion_end -
+                        si.job.inversion_start):
+                        longest_inversions.append(record)
+                        longest_inversions = _sort_longest_inversions(
+                            longest_inversions)
+                        del longest_inversions[0]
+                else:
+                    longest_inversions.append(record)
+                    longest_inversions = _sort_longest_inversions(
+                        longest_inversions)
+
+    # We've seen all records.
+    # Further update state
+    if num_inversions > 0:
+        avg_inversion = int(sum_inversions / num_inversions)
+    else:
+        avg_inversion = 0
+
+    # Print out our information
+    print "Num inversions: %d" % (num_inversions)
+    print "Min inversion: %d" % (min_inversion)
+    print "Max inversion: %d" % (max_inversion)
+    print "Avg inversion: %d" % (avg_inversion)
+    for inv in longest_inversions:
+        print ""
+        print "Inversion record IDs: (%d, %d)" % (inv.inversion_start_id,
+            inv.id)
+        print("Triggering Event IDs: (%d, %d)" %
+            (inv.inversion_start_triggering_event_id,
+            inv.triggering_event_id))
+        print "Time: %d" % (inv.job.inversion_end)
+        print "Duration: %d" % (
+            inv.job.inversion_end - inv.job.inversion_start)
+        print "Job: %d.%d" % (inv.job.pid,inv.job.job)
+        print "Deadline: %d" % (inv.job.deadline)
+        print "Eligible: ",
+        for job in inv.eligible:
+            print str(job) + " ",
+        print
+        print "On CPU: ",
+        for job in inv.on_cpu:
+            print str(job) + " ",
+        print ""
+
+
+def _sort_longest_inversions(longest_inversions):
+    """ Sort longest inversions"""
+    def sortkey(x):
+        return x.job.inversion_end - x.job.inversion_start
+    longest_inversions.sort(key=sortkey)
+    return longest_inversions