AENC/switchchain Changeset - 28b33414825e · Centrum Wiskunde & Informatica (CWI)

Changeset - 28b33414825e

Parent rev.

Child rev.

[Not reviewed]

0 2 0

Tom Bannink - 8 years ago 2017-05-15 17:13:00
tombannink@gmail.com

Add initial code for degree-structure-plots

2 files changed with 26 insertions and 1 deletions:

cpp/graph.hpp

cpp/switchchain.cpp

0 comments (0 inline, 0 general)

cpp/graph.hpp

➞

Show inline comments

@@ @@ -56,6 +56,8 @@ class Graph { @@
     const Edge &getEdge(unsigned int i) const { return edges[i].e; }
     const auto& getAdj() const { return adj; }
     // When the degree sequence is not graphics, the Graph can be
     // in any state, it is not neccesarily empty
     bool createFromDegreeSequence(const DegreeSequence &d) {

cpp/switchchain.cpp

➞

Show inline comments

@@ @@ -2,6 +2,7 @@ @@
 #include "graph.hpp"
 #include "powerlaw.hpp"
 #include <algorithm>
 #include <array>
 #include <fstream>
 #include <iostream>
 #include <numeric>
@@ @@ -63,6 +64,27 @@ class SwitchChain { @@
     std::bernoulli_distribution permutationDistribution;
 };
 void getTriangleDegrees(const Graph& g) {
     std::vector<std::array<std::size_t,3>> trids;
     const auto& adj = g.getAdj();
     int triangles = 0;
     for (auto& v : adj) {
         for (unsigned int i = 0; i < v.size(); ++i) {
             for (unsigned int j = i + 1; j < v.size(); ++j) {
                 if (g.hasEdge({v[i], v[j]})) {
                     ++triangles;
                     std::array<std::size_t, 3> ds = {v.size(), adj[v[i]].size(),
                                                      adj[v[j]].size()};
                     std::sort(ds.begin(), ds.end());
                     trids.push_back(ds);
+                }
+            }
+        }
+    }
     assert(triangles % 3 == 0);
     return;
+}
 //
 // Assumes degree sequence does NOT contain any zeroes!
 //
@@ @@ -201,7 +223,7 @@ int main() { @@
             // For a single n,tau take samples over several instances of
             // the degree distribution.
             // 500 samples seems to give reasonable results
-            for (int degreeSample = 0; degreeSample < 200; ++degreeSample) {
+            for (int degreeSample = 0; degreeSample < 1; ++degreeSample) {
                 // Generate a graph
                 // might require multiple tries
                 for (int i = 1; ; ++i) {
@@ @@ -277,6 +299,7 @@ int main() { @@
                 constexpr int measureSkip = 1;
                 int movesDone = 0;
                 int triangles[measurements];

0 comments (0 inline, 0 general)