diff --git a/.gitignore b/.gitignore
index 9d73845..254131b 100644
--- a/.gitignore
+++ b/.gitignore
@@ -4,3 +4,4 @@ make_graphs/bin
 make_graphs/test_graphs
 make_planar/bin
 test_area
+dpt
diff --git a/deferred_planarity_test/main.cpp b/deferred_planarity_test/main.cpp
index 9ff7ed2..c713148 100644
--- a/deferred_planarity_test/main.cpp
+++ b/deferred_planarity_test/main.cpp
@@ -2,15 +2,32 @@
 // A simple code that test the MPS algorighm.
 //-----------------------------------------------------------------------------------
 
+#include "mps.h"
 #include <iostream>
 #include <cstdlib>
 #include <climits>
 #include <string>
-#include "mps.h"
+#include <algorithm>
+#include <iterator>
+#include <random>
+#include <vector>
+ 
+#include <ogdf/fileformats/GraphIO.h>
+
 using namespace std;
 
-int find_mps(string input_file);
+int compute_removed_edge_size(string input_file, vector<int> post_order);
 
+int get_graph_size(string input_file) {
+    ogdf::Graph G;
+
+    // utilize OGDF readGML
+    if (!ogdf::GraphIO::read(G, input_file, ogdf::GraphIO::readGML)) {
+        std::cerr << "Could not read " << input_file << ".gml" << std::endl;
+    }
+
+    return G.numberOfNodes();
+}
 
 //-----------------------------------------------------------------------------------
 // Main function.
@@ -19,7 +36,25 @@ int find_mps(string input_file);
 
 int main(int argc, char* argv[]) {
     string input_file = argv[1];
-	std::cout << "Number of removed edges: " << find_mps(input_file) << std::endl;
+    int num_runs = stoi(argv[2]);
+    // find the size of the graph here
+    int node_size = get_graph_size(input_file);
+    // generate order here
+    vector<int> post_order;
+    for (int i=0; i < node_size; ++i) {
+        post_order.push_back(i);
+    }
+
+    for (int i=0; i< num_runs; ++i) {
+        std::random_device rd;
+        std::mt19937 g(rd());
+    
+        std::shuffle(post_order.begin(), post_order.end(), g);
+        // print order
+        std::copy(post_order.begin(), post_order.end(), std::ostream_iterator<int>(std::cout, " "));
+        std::cout << std::endl;
+        std::cout << "Number of removed edges: " << compute_removed_edge_size(input_file, post_order) << std::endl;
+    }
 
 	return 0;
 }
diff --git a/deferred_planarity_test/mps.cpp b/deferred_planarity_test/mps.cpp
index e5cbb20..0dc746c 100644
--- a/deferred_planarity_test/mps.cpp
+++ b/deferred_planarity_test/mps.cpp
@@ -31,6 +31,14 @@ maximal_planar_subgraph_finder::postOrderTraversal() {
 	}
 }
 
+//Set the post-order-list via given list
+void
+maximal_planar_subgraph_finder::set_post_order(vector<int> post_order) {
+    for (int i = 0; i < _node_list.size(); ++i) {
+        _node_list[i]->set_post_order_index(post_order[i]);
+    }
+}
+
 //Sort the adj-list of every node increasingly according to post-order-index.
 void
 maximal_planar_subgraph_finder::sort_adj_list() {
diff --git a/deferred_planarity_test/mps.h b/deferred_planarity_test/mps.h
index e43b745..8a9dd6c 100644
--- a/deferred_planarity_test/mps.h
+++ b/deferred_planarity_test/mps.h
@@ -142,10 +142,12 @@ public:
 	maximal_planar_subgraph_finder();
 	~maximal_planar_subgraph_finder();
 	int find_mps(string input_file);
+	int compute_removed_edge_size(string input_file, vector<int> post_order);
 	node* get_new_node(node_type t);
     void read_from_gml(string input_file);
 	int output_removed_edge_size();
 	void postOrderTraversal();
+    void set_post_order(vector<int> post_order);
 	void sort_adj_list();
 	void determine_edges();
 	void back_edge_traversal();
diff --git a/deferred_planarity_test/mps_test.cpp b/deferred_planarity_test/mps_test.cpp
index 819954a..795ca0c 100644
--- a/deferred_planarity_test/mps_test.cpp
+++ b/deferred_planarity_test/mps_test.cpp
@@ -13,6 +13,11 @@ int find_mps(string input_file) {
 	return m.find_mps(input_file);
 }
 
+int compute_removed_edge_size(string input_file, vector<int> post_order) {
+    maximal_planar_subgraph_finder m;
+    return m.compute_removed_edge_size(input_file, post_order);
+}
+
 int maximal_planar_subgraph_finder::find_mps(string input_file) {
 	read_from_gml(input_file);
 	postOrderTraversal();
@@ -22,6 +27,16 @@ int maximal_planar_subgraph_finder::find_mps(string input_file) {
 	return output_removed_edge_size();
 }
 
+
+int maximal_planar_subgraph_finder::compute_removed_edge_size(string input_file, vector<int> post_order) {
+	read_from_gml(input_file);
+    set_post_order(post_order);
+	sort_adj_list();
+	determine_edges();
+	back_edge_traversal();
+	return output_removed_edge_size();
+}
+
 //-----------------------------------------------------------------------------------
 // Imput, output
 //-----------------------------------------------------------------------------------