//Compilation: g++ example.cpp specgraph.cpp rand.cpp or simply 'make'

#include<iostream>
#include "specgraph.h"
#include "rand.h"

using namespace std;

int main()
{
 randomize();
 //Parameters of the algorithm
 specgraph_probaerror=0.01;
 specgraph_maxiter=100000;
 specgraph_stopthreshold=0.001;
	
 //Graph construction; n = number of vertices
 const int n=500;
 Graph G(n);
 int i,j,k;

 //The cyclic graph
 //for(i=0;i<n;i++)G.AddSymEdge(i,(i+1)%n);

 //The complete graph
 //for(i=0;i<n;i++)for(j=i+1;j<n;j++)G.AddSymEdge(i,j);

 //A random graph of average degree 2d
 const int d=10;
		for(k=0;k<d*n;k++){
		i=alea(n);j=alea(n);G.AddSymEdge(i,j);
		}

 //Other graphs:
 // for(i=1;i<n;i++)G.AddSymEdge(0,i);
 //for(i=1;i<n;i+=2)G.AddSymEdge(i,i+1);

 //cout<<"Bounds for lambda1: "<<endl;
 //G.Lambda1();
 cout<<"Bounds for specrad: "<<endl;
 G.SpecRad();

 return 0;
}