///Generate graph using the AGM model. CProbV = vector of Pc
PUNGraph TAGM::GenAGM(TVec<TIntV>& CmtyVV, const TFltV& CProbV, TRnd& Rnd, const double PNoCom) {
PUNGraph G = TUNGraph::New(100 * CmtyVV.Len(), -1);
printf("AGM begins\n");
for (int i = 0; i < CmtyVV.Len(); i++) {
TIntV& CmtyV = CmtyVV[i];
for (int u = 0; u < CmtyV.Len(); u++) {
if ( G->IsNode(CmtyV[u])) {
continue;
}
G->AddNode(CmtyV[u]);
}
double Prob = CProbV[i];
RndConnectInsideCommunity(G, CmtyV, Prob, Rnd);
}
if (PNoCom > 0.0) { //if we want to connect nodes that do not share any community
TIntSet NIDS;
for (int c = 0; c < CmtyVV.Len(); c++) {
for (int u = 0; u < CmtyVV[c].Len(); u++) {
NIDS.AddKey(CmtyVV[c][u]);
}
}
TIntV NIDV;
NIDS.GetKeyV(NIDV);
RndConnectInsideCommunity(G,NIDV,PNoCom,Rnd);
}
printf("AGM completed (%d nodes %d edges)\n",G->GetNodes(),G->GetEdges());
G->Defrag();
return G;
}
PUNGraph TAGM::GenAGM(TVec<TIntV>& CmtyVV, const double& DensityCoef, const double& ScaleCoef, TRnd& Rnd){
TFltV CProbV;
double Prob;
for(int i=0;i<CmtyVV.Len();i++) {
Prob = ScaleCoef*pow(double(CmtyVV[i].Len()),-DensityCoef);
if(Prob>1){Prob = 1;}
CProbV.Add(Prob);
}
PUNGraph G = TUNGraph::New();
printf("AGM begins\n");
for(int i=0;i<CmtyVV.Len();i++) {
TIntV& CmtyV = CmtyVV[i];
for(int u=0;u<CmtyV.Len();u++) {
G->AddNode(CmtyV[u]);
}
Prob = CProbV[i];
printf("\r%d(%d)/%d",i,CmtyVV[i].Len(),CmtyVV.Len());
RndConnectInsideCommunity(G,CmtyV,Prob,Rnd);
}
printf("AGM completed (%d nodes %d edges)\n",G->GetNodes(),G->GetEdges());
return G;
}
