summaryrefslogtreecommitdiffstats
path: root/src/libcola/tests/scale_free.cpp
diff options
context:
space:
mode:
authorMarc Jeanmougin <marc.jeanmougin@telecom-paristech.fr>2018-04-29 14:25:32 +0000
committerMarc Jeanmougin <marc.jeanmougin@telecom-paristech.fr>2018-04-29 14:25:32 +0000
commitab5f8ff5869021958f4ae8b838c3d707a2e85eaa (patch)
tree4907675828a5401d013b7587538cc8541edd2764 /src/libcola/tests/scale_free.cpp
parentmoved libcroco, libuemf, libdepixelize to 3rdparty folder (diff)
downloadinkscape-ab5f8ff5869021958f4ae8b838c3d707a2e85eaa.tar.gz
inkscape-ab5f8ff5869021958f4ae8b838c3d707a2e85eaa.zip
Put adaptagrams into its own folder
Diffstat (limited to 'src/libcola/tests/scale_free.cpp')
-rw-r--r--src/libcola/tests/scale_free.cpp135
1 files changed, 0 insertions, 135 deletions
diff --git a/src/libcola/tests/scale_free.cpp b/src/libcola/tests/scale_free.cpp
deleted file mode 100644
index d9c2934e7..000000000
--- a/src/libcola/tests/scale_free.cpp
+++ /dev/null
@@ -1,135 +0,0 @@
-/*
- * vim: ts=4 sw=4 et tw=0 wm=0
- *
- * libcola - A library providing force-directed network layout using the
- * stress-majorization method subject to separation constraints.
- *
- * Copyright (C) 2006-2008 Monash University
- *
- * This library is free software; you can redistribute it and/or
- * modify it under the terms of the GNU Lesser General Public
- * License as published by the Free Software Foundation; either
- * version 2.1 of the License, or (at your option) any later version.
- *
- * This library is distributed in the hope that it will be useful,
- * but WITHOUT ANY WARRANTY; without even the implied warranty of
- * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
- * Lesser General Public License for more details.
- *
- * You should have received a copy of the GNU Lesser General Public
- * License along with this library in the file LICENSE; if not,
- * write to the Free Software Foundation, Inc., 59 Temple Place,
- * Suite 330, Boston, MA 02111-1307 USA
- *
-*/
-
-// generates a random graph with a power-law distribution of node degrees
-// the algorithm is pretty much the Barabasi-Albert model but with an extra step
-// so that I can guarantee the number of nodes in the graph.
-// We begin with an initial chain of m nodes.
-// Then, we add new nodes, connecting them to existing nodes v with probability
-// degree(v)/(sum all degrees)
-// We perform an extra step, looking for disconnected nodes, connecting them
-// to other nodes again with probability dependent on degree
-#include<iostream>
-#include<vector>
-#include <cmath>
-#include <time.h>
-#include <valarray>
-
-#include "graphlayouttest.h"
-
-void scale_free(const unsigned n = 50) {
- const unsigned m = 3; // begin with a chain of m nodes
- const double g=0.08; // edge density, every time we add a node we add 1+g*rand(0-1) edges
- // connecting them to existing nodes with probability based on degree
- unsigned d[n];
- double defaultEdgeLength=100;
- vector<Edge> es;
- CompoundConstraints cx,cy;
- for(unsigned i=0;i<n;i++) {
- d[i]=0;
- }
- unsigned sumdegree=0;
- for(unsigned i=1;i<m;i++) {
- es.push_back(make_pair(i-1,i));
- d[i]++;d[i-1]++;
- sumdegree+=2;
- }
- srand(3);
- for(unsigned i=2;i<n;i++) {
- for(unsigned j=0;j<i;j++) {
- double p=(double)d[j]/(double)sumdegree;
- double r=(double)rand()/(double)RAND_MAX;
- if(r*p>g) {
- es.push_back(make_pair(j,i));
- d[j]++;d[i]++;
- sumdegree+=2;
- }
- }
- }
- for(unsigned i=0;i<n;i++) {
- if(d[i]==0) {
- double maxP=0;
- unsigned end=0;
- for(unsigned j=0;j<n;j++) {
- if(j==i) continue;
- double p=(double)d[j]/(double)sumdegree;
- double r=(double)rand()/(double)RAND_MAX;
- if(r*p>maxP) {
- maxP=r*p;
- end=j;
- }
- }
- if(end>i) {
- es.push_back(make_pair(i,end));
- } else {
- es.push_back(make_pair(end,i));
- }
- d[end]++;d[i]++;
- sumdegree+=2;
- }
- }
- valarray<double> eweights(es.size());
- sort(d,d+n);
- unsigned degree=0,ctr=0;
- printf("degree distribution:\n");
- for(unsigned i=0;i<n;i++) {
- if(degree!=d[i]) {
- if(degree!=0) printf("%d,%d\n",degree,ctr);
- degree=d[i];
- ctr=0;
- }
- ctr++;
- }
- printf("%d,%d\n",degree,ctr);
- for(unsigned i=0;i<es.size();i++) {
- unsigned a=es[i].first, b=es[i].second;
- eweights[i]=1;
- cy.push_back(new SeparationConstraint(a,b,defaultEdgeLength/3));
- }
-
- cout << "|V|="<<n<<endl;
- double width=1000;
- double height=1000;
- vector<pair<double,double> > startpos(n);
- for(unsigned i=0;i<n;i++) {
- double x=getRand(width), y=getRand(height);
- startpos[i]=make_pair(x,y);
- }
- const char *f="scalefree";
- //run_test(startpos,es,defaultEdgeLength,cx,cy,CG,false,f,"cg");
- //run_test(startpos,es,defaultEdgeLength,cx,cy,IP,false,f,"ip");
- //run_test(startpos,es,defaultEdgeLength,cx,cy,UGP,false,f,"ugp");
- //run_test(startpos,es,defaultEdgeLength,cx,cy,SGP,false,f,"sgp");
- //run_test(startpos,es,defaultEdgeLength,cx,cy,IP,true,f,"cip");
- //run_test(startpos,es,defaultEdgeLength,cx,cy,UGP,true,f,"cugp");
- run_test(startpos,es,defaultEdgeLength,cx,cy,SGP,true,f,"csgp");
-}
-int main() {
- for(unsigned i=1;i<2;i++) {
- scale_free(i*50);
- }
- return 0;
-}
-// vim: filetype=cpp:expandtab:shiftwidth=4:tabstop=4:softtabstop=4:textwidth=99 :