1 files changed, 0 insertions, 102 deletions
diff --git a/src/libcola/shortest_paths.cpp b/src/libcola/shortest_paths.cpp
deleted file mode 100644
index 514721fb5..000000000
--- a/src/libcola/shortest_paths.cpp
+++ /dev/null
@@ -1,102 +0,0 @@
-// vim: set cindent 
-// vim: ts=4 sw=4 et tw=0 wm=0
-#include "shortest_paths.h"
-#include <float.h>
-#include <cassert>
-#include <iostream>
-#include <libvpsc/pairingheap/PairingHeap.h>
-using namespace std;
-namespace shortest_paths {
-// O(n^3) time.  Slow, but fool proof.  Use for testing.
-void floyd_warshall(
-        unsigned n,
-        double** D, 
-        vector<Edge>& es,
-        double* eweights) 
-{
-    for(unsigned i=0;i<n;i++) {
-        for(unsigned j=0;j<n;j++) {
-            if(i==j) D[i][j]=0;
-            else D[i][j]=DBL_MAX;
-        }
-    }
-    for(unsigned i=0;i<es.size();i++) {
-        unsigned u=es[i].first, v=es[i].second;
-        assert(u<n&&v<n);
-        D[u][v]=D[v][u]=eweights[i];
-    }
-    for(unsigned k=0; k<n; k++) {
-        for(unsigned i=0; i<n; i++) {
-            for(unsigned j=0; j<n; j++) {
-                D[i][j]=min(D[i][j],D[i][k]+D[k][j]);
-            }
-        }
-    }
-}
-static void dijkstra_init(Node* vs, vector<Edge>& es, double* eweights) {
-    for(unsigned i=0;i<es.size();i++) {
-        unsigned u=es[i].first, v=es[i].second;
-        vs[u].neighbours.push_back(&vs[v]);
-        vs[u].nweights.push_back(eweights[i]);
-        vs[v].neighbours.push_back(&vs[u]);
-        vs[v].nweights.push_back(eweights[i]);
-    }
-}
-static void dijkstra(
-        unsigned s,
-        unsigned n,
-        Node* vs,
-        double* d)
-{
-    assert(s<n);
-    for(unsigned i=0;i<n;i++) {
-        vs[i].id=i;
-        vs[i].d=DBL_MAX;
-        vs[i].p=NULL;
-    }
-    vs[s].d=0;
-    PairingHeap<Node*> Q(&compareNodes);
-    for(unsigned i=0;i<n;i++) {
-        vs[i].qnode = Q.insert(&vs[i]);
-    }
-    while(!Q.isEmpty()) {
-        Node *u=Q.extractMin();
-        d[u->id]=u->d;
-        for(unsigned i=0;i<u->neighbours.size();i++) {
-            Node *v=u->neighbours[i];
-            double w=u->nweights[i];
-            if(v->d>u->d+w) {
-                v->p=u;
-                v->d=u->d+w;
-                Q.decreaseKey(v->qnode,v);
-            }
-        }
-    }
-}
-
-void dijkstra(
-        unsigned s,
-        unsigned n,
-        double* d,
-        vector<Edge>& es,
-        double* eweights)
-{
-    assert(s < n);
-    std::vector<Node> vs(n);
-    dijkstra_init(&vs[0], es, eweights);
-    dijkstra(s, n, &vs[0], d);
-}
-
-void johnsons(
-        unsigned n,
-        double** D, 
-        vector<Edge>& es,
-        double* eweights) 
-{
-    std::vector<Node> vs(n);
-    dijkstra_init(&vs[0], es, eweights);
-    for (unsigned k = 0; k < n; k++) {
-        dijkstra(k,n,&vs[0],D[k]);
-    }
-}
-}