Vector Matrix::rowsums() const {
Vector vector(ncol());
SafeCall(igraph_matrix_rowsum(ptr(), vector.ptr()));
return vector;
}
/*__________________________________________________________________________ MAIN CYCLE (FLTK CYCLE) ____*/
void mainidle_cb(void*){ //this routine updates the program.
//thus, it computes the EVOLUTION
double shooted;
double dens, err;
double totdens, toterr, totimerr;
char s[100];
// ---- running controls AND PRINTING
if(
(amstepping==0 && runningcontrol==1 && graphisloaded==1 && ticks<=maxtime ) ||
(amstepping==1 && runningcontrol==1 && graphisloaded==1 && ticks<=maxtime && tickstep<=step-1)
)
{
Evolution(deltat);
//PRINTS
if((int)printdatabutton->value()==1){
//if have steady state
if(usesteady==1){
igraph_matrix_t activation;
igraph_matrix_init(&activation,nodesnumber,totrun);
igraph_matrix_null(&activation);
igraph_vector_t correlation;
igraph_vector_init(&correlation,(nodesnumber*nodesnumber)); igraph_vector_null(&correlation);
fprintf(output1,"%i ", ticks);
fprintf(output2,"%i ", ticks);
fprintf(output5,"%i ", ticks);
fprintf(output6,"%i ", ticks);
totdens=0; toterr=0;
for(int i=0;i<nodesnumber;++i){
shooted=0;
dens=0;
err=0;
for(int j=0; j<totrun; ++j){
dens=dens+MATRIX(density,i,j);
err=err+((VECTOR(statstate)[i]-MATRIX(density,i,j))*(VECTOR(statstate)[i]-MATRIX(density,i,j)));
shooted=shooted+MATRIX(loss,i,j);
if(MATRIX(loss,i,j)!=0){++MATRIX(activation,i,j);}
}
dens=dens/totrun;
err=sqrt(err)/totrun;
shooted=shooted/totrun;
totdens=totdens+dens;
toterr=toterr+err;
fprintf(output1,"%f ",dens);
fprintf(output2,"%f ",err);
fprintf(output5,"%f ",shooted);
}
totdens=totdens/nodesnumber;
toterr=toterr/nodesnumber;
fprintf(output1,"%f ",totdens);
fprintf(output2,"%f ",toterr);
// printf("\n\n ACTIVATION MATRIX \n \n"); print_matrix_ur(&activation,stdout); printf("\n\n");
// ---- CORRELATION ---
igraph_vector_t meanactivation;
igraph_vector_init(&meanactivation,nodesnumber);
igraph_vector_null(&meanactivation);
//calculate mean activation
for (int j=0; j<totrun; ++j) {
for (int i=0; i<nodesnumber; ++i) {
VECTOR(meanactivation)[i]=VECTOR(meanactivation)[i]+MATRIX(activation,i,j);
}
}
igraph_vector_scale(&meanactivation,1./totrun);
//calculate actual correlation
for (int x=0; x<nodesnumber ; ++x) {
for(int y=0; y<nodesnumber; ++y){
double prod=0;
for (int j=0; j<totrun; ++j) {
prod=prod+ ( MATRIX(activation,x,j)*MATRIX(activation,y,j) );
}
prod=prod/totrun;
VECTOR(correlation)[(x*nodesnumber+y)] = prod - (VECTOR(meanactivation)[x]*VECTOR(meanactivation)[y]);
}
}
igraph_vector_destroy(&meanactivation);
for (int i=0; i<(nodesnumber*nodesnumber); ++i) {
fprintf(output6,"%f ",VECTOR(correlation)[i]);
}
//calculate error on run
igraph_matrix_t distl1; igraph_matrix_t distimel1;
igraph_matrix_init(&distl1,nodesnumber,totrun); igraph_matrix_init(&distimel1,nodesnumber,totrun);
igraph_matrix_null(&distl1); igraph_matrix_null(&distimel1);
igraph_vector_t rundistl1; igraph_vector_t rundistimel1;
igraph_vector_init(&rundistl1,totrun); igraph_vector_init(&rundistimel1,totrun);
igraph_vector_null(&rundistl1); igraph_vector_null(&rundistimel1);
toterr=0; totimerr=0;
//for every run
for(int j=0;j<totrun;++j){
//i evaluate the distance between the state and the stationary state (toterr) and the distance between old and new density (totimerr)
for(int i=0; i<nodesnumber; ++i)
{
//L1 DISTANCE WRT STATSTATE & DENSITY
MATRIX(distl1,i,j)=fabs(VECTOR(statstate)[i]-MATRIX(density,i,j));
//L1 DISTANCE WRT OLD DENSITY & DENSITY
MATRIX(distimel1,i,j)=fabs(MATRIX(densityold,i,j)-MATRIX(density,i,j));
}
}
igraph_matrix_rowsum(&distl1,&rundistl1); igraph_matrix_rowsum(&distimel1,&rundistimel1);
igraph_vector_scale(&rundistl1,(1./nodesnumber)); igraph_vector_scale(&rundistimel1,(1./nodesnumber));
toterr= (double)( igraph_vector_sum(&rundistl1) ) / (double)totrun ;
totimerr= (double)( igraph_vector_sum(&rundistimel1)) / (double)totrun;
igraph_vector_destroy(&rundistl1); igraph_vector_destroy(&rundistimel1);
igraph_matrix_destroy(&distl1); igraph_matrix_destroy(&distimel1);
fprintf(output2,"%f %f",toterr, totimerr);
fprintf(output1,"\n");
fprintf(output2,"\n");
fprintf(output5,"\n");
fprintf(output6,"\n");
//if i have BRIDGES ("clustered" graph), I print the traffic on the BRIDGES, using "output3" file
if(isclustered==1){
//for each bridge
fprintf(output3,"%i ",ticks);
for(int nbri=0; nbri<igraph_matrix_ncol(&bridgeslinks); ++nbri){
for(int i=0; i<igraph_matrix_nrow(&bridgeslinks);++i){
double tfl=0;
for(int j=0; j<totrun; ++j){
int beid;
beid=(int)MATRIX(bridgeslinks,i,nbri);
tfl=tfl+MATRIX(flux,beid,j);
}
fprintf(output3,"%f ",tfl);
}
}
fprintf(output3,"\n");
}
igraph_matrix_destroy(&activation);
igraph_vector_destroy(&correlation);
}
//if i HAVENT STEADY STATE
else {
fprintf(output1,"%i ", ticks);
fprintf(output5,"%i ", ticks);
for(int i=0;i<nodesnumber;++i){
shooted=0;
dens=0;
for(int j=0; j<totrun; ++j){
dens=dens+MATRIX(density,i,j);
shooted=shooted+MATRIX(loss,i,j);
}
dens=dens/totrun;
shooted=shooted/totrun;
fprintf(output1,"%f " ,dens);
fprintf(output5,"%f " ,shooted);
}
fprintf(output1,"\n");
fprintf(output5,"\n");
}
}
}
if((ticks==maxtime || tickstep==step) && runningcontrol==1){
run();
}
// ---- no graph loaded
if(graphisloaded==0 && rewrite==1) {
runbutton->deactivate();
sprintf(s,"No\nnetwork\nloaded");
databuff->text(s);
}
// ---- graph loaded
else {
runbutton->activate();
if(islattice==1 && rewrite==1){
if(istoro==1){
sprintf(s,"Nodes=%i\nToroidal\nLattice\n%iD Side=%i",nodesnumber, latticedim, latticeside);
}
else{
sprintf(s,"Nodes=%i\nLattice\n%iD Side=%i",nodesnumber, latticedim, latticeside);
}
databuff->text(s);
}
else if(rewrite==1){
sprintf(s,"Nodes=%i",nodesnumber);
databuff->text(s);
}
}
//have path
if(havepath==1 && rewrite==1){pathbuff->text(path); }
else if(havepath==0 && rewrite==1){pathbuff->text("No Path");}
if(error==1 && rewrite==1){
sprintf(s,errorstring);
databuff->text(s);
}
if (ticks<=maxtime){
scene->redraw();
datascene->redraw();
}
rewrite=0;
//Fl::repeat_timeout(1.0, mainidle_cb);
}
int main() {
igraph_matrix_t m, m2;
igraph_vector_t v;
long int i, j, i2, j2;
igraph_real_t r1, r2;
igraph_matrix_init(&m, 4, 3);
byrow(&m);
/* igraph_matrix_e */
printf("igraph_matrix_e\n");
apply(m, printf("%i ", (int)igraph_matrix_e(&m, i, j)), printf("\n"));
/* igraph_matrix_e_ptr */
printf("igraph_matrix_e_ptr\n");
apply(m, printf("%i ", (int)igraph_matrix_e_ptr(&m, i, j)[0]), printf("\n"));
/* igraph_matrix_set */
printf("igraph_matrix_set\n");
apply(m, igraph_matrix_set(&m, i, j, i), (void) 0 );
print_matrix(&m);
apply(m, igraph_matrix_set(&m, i, j, j), (void) 0 );
print_matrix(&m);
/* igraph_matrix_fill */
printf("igraph_matrix_fill\n");
igraph_matrix_fill(&m, 42);
print_matrix(&m);
igraph_matrix_fill(&m, -42.1);
print_matrix(&m);
/* igraph_matrix_update */
printf("igraph_matrix_update\n");
igraph_matrix_init(&m2, 0, 0);
byrow(&m);
igraph_matrix_update(&m2, &m);
print_matrix(&m2);
/* igraph_matrix_rbind */
printf("igraph_matrix_rbind\n");
igraph_matrix_rbind(&m2, &m);
print_matrix(&m2);
printf("\n");
igraph_matrix_resize(&m, 0, igraph_matrix_ncol(&m2));
igraph_matrix_rbind(&m2, &m);
print_matrix(&m2);
printf("\n");
igraph_matrix_rbind(&m, &m2);
print_matrix(&m);
/* igraph_matrix_cbind */
printf("igraph_matrix_cbind\n");
igraph_matrix_resize(&m, 4, 3);
igraph_matrix_resize(&m2, 4, 2);
byrow(&m);
byrow(&m2);
igraph_matrix_cbind(&m, &m2);
print_matrix(&m);
/* igraph_matrix_swap */
printf("igraph_matrix_swap\n");
igraph_matrix_update(&m, &m2);
igraph_matrix_null(&m);
igraph_matrix_swap(&m, &m2);
print_matrix(&m);
print_matrix(&m2);
/* igraph_matrix_get_row */
/* igraph_matrix_set_row */
printf("igraph_matrix_get_row\n");
printf("igraph_matrix_set_row\n");
igraph_vector_init(&v, 0);
for (i=0; i<igraph_matrix_nrow(&m); i++) {
igraph_matrix_get_row(&m, &v, i);
igraph_matrix_set_row(&m2, &v, i);
}
print_matrix(&m2);
/* igraph_matrix_set_col */
printf("igraph_matrix_set_col\n");
igraph_matrix_null(&m2);
for (i=0; i<igraph_matrix_ncol(&m); i++) {
igraph_matrix_get_col(&m, &v, i);
igraph_matrix_set_col(&m2, &v, i);
}
print_matrix(&m2);
/* igraph_matrix_swap_rows */
printf("igraph_matrix_swap_rows\n");
igraph_matrix_swap_rows(&m2, 0, 0);
igraph_matrix_swap_rows(&m2, 0, 2);
print_matrix(&m2);
/* igraph_matrix_swap_cols */
printf("igraph_matrix_swap_cols\n");
igraph_matrix_swap_cols(&m2, 0, 0);
igraph_matrix_swap_cols(&m2, 0, 1);
print_matrix(&m2);
/* igraph_matrix_add_constant */
printf("igraph_matrix_add_constant\n");
igraph_matrix_add_constant(&m2, 0);
print_matrix(&m2);
igraph_matrix_add_constant(&m2, -1);
print_matrix(&m2);
/* igraph_matrix_add */
printf("igraph_matrix_add\n");
byrow(&m2);
byrow(&m);
igraph_matrix_add(&m2, &m);
print_matrix(&m2);
/* igraph_matrix_sub */
printf("igraph_matrix_sub\n");
igraph_matrix_sub(&m2, &m);
print_matrix(&m2);
/* igraph_matrix_mul_elements */
printf("igraph_matrix_mul_elements\n");
igraph_matrix_mul_elements(&m2, &m);
print_matrix(&m2);
/* igraph_matrix_div_elements */
printf("igraph_matrix_div_elements\n");
igraph_matrix_fill(&m, 2);
igraph_matrix_div_elements(&m2, &m);
print_matrix(&m2);
/* igraph_matrix_min */
printf("igraph_matrix_min\n");
if (igraph_matrix_min(&m2) != 0) {
return 1;
}
if (igraph_matrix_min(&m) != 2) {
return 1;
}
/* igraph_matrix_which_min */
printf("igraph_matrix_which_min\n");
igraph_matrix_which_min(&m2, &i, &j);
if (i != 0 || j != 0) { return 2; }
MATRIX(m2,0,1) = -1;
igraph_matrix_which_min(&m2, &i, &j);
if (i != 0 || j != 1) { return 2; }
MATRIX(m2,3,1) = -2;
igraph_matrix_which_min(&m2, &i, &j);
if (i != 3 || j != 1) { return 2; }
/* igraph_matrix_which_max */
printf("igraph_matrix_which_max\n");
MATRIX(m2,3,0) = 100;
igraph_matrix_which_max(&m2, &i, &j);
if (i != 3 || j != 0) { return 3; }
/* igraph_matrix_minmax */
printf("igraph_matrix_minmax\n");
igraph_matrix_minmax(&m2, &r1, &r2);
printf("%g %g\n", r1, r2);
/* igraph_matrix_which_minmax */
printf("igraph_matrix_which_minmax\n");
igraph_matrix_which_minmax(&m2, &i, &j, &i2, &j2);
if (i != 3 || j != 1 || i2 != 3 || j2 != 0) { return 4; }
/* igraph_matrix_isnull */
printf("igraph_matrix_isnull\n");
if (igraph_matrix_isnull(&m2)) { return 5; }
igraph_matrix_null(&m);
if (!igraph_matrix_isnull(&m)) { return 5; }
igraph_matrix_resize(&m2, 5, 0);
if (!igraph_matrix_isnull(&m2)) { return 5; }
/* igraph_matrix_empty */
printf("igraph_matrix_empty\n");
if (!igraph_matrix_empty(&m2)) { return 6; }
igraph_matrix_resize(&m2, 5, 5);
if (igraph_matrix_empty(&m2)) { return 6; }
/* igraph_matrix_is_symmetric */
printf("igraph_matrix_is_symmetric\n");
byrow(&m2);
if (igraph_matrix_is_symmetric(&m2)) { return 7; }
igraph_matrix_update(&m, &m2);
igraph_matrix_transpose(&m);
igraph_matrix_add(&m, &m2);
if (!igraph_matrix_is_symmetric(&m)) { return 7; }
/* igraph_matrix_prod */
printf("igraph_matrix_prod\n");
igraph_matrix_resize(&m, 3,2);
byrow(&m);
igraph_matrix_add_constant(&m, 1);
print_matrix(&m);
printf("product: %g\n", igraph_matrix_prod(&m));
/* igraph_matrix_rowsum */
printf("igraph_matrix_rowsum\n");
igraph_matrix_rowsum(&m, &v);
print_vector(&v);
/* igraph_matrix_colsum */
printf("igraph_matrix_colsum\n");
igraph_matrix_colsum(&m, &v);
print_vector(&v);
/* igraph_matrix_contains */
printf("igraph_matrix_contains\n");
if (igraph_matrix_contains(&m, 0)) { return 8; }
if (igraph_matrix_contains(&m, 6.0001)) { return 8; }
if (igraph_matrix_contains(&m, 7)) { return 8; }
if (!igraph_matrix_contains(&m, 1)) { return 8; }
if (!igraph_matrix_contains(&m, 6)) { return 8; }
/* igraph_matrix_search */
printf("igraph_matrix_search\n");
if (!igraph_matrix_search(&m, 0, 6.0, &i2, &i, &j)) { return 9; }
if (i2 != 5 || i != 2 || j != 1) { return 9; }
/* igraph_matrix_remove_row */
printf("igraph_matrix_remove_row\n");
igraph_matrix_remove_row(&m, 1);
print_matrix(&m);
igraph_matrix_resize(&m,5,4);
byrow(&m);
igraph_matrix_remove_row(&m, 4);
print_matrix(&m);
igraph_matrix_remove_row(&m, 0);
print_matrix(&m);
/* igraph_matrix_select_cols */
printf("igraph_matrix_select_cols\n");
igraph_matrix_resize(&m, 6, 5);
apply(m, igraph_matrix_set(&m, i, j, j), (void) 0 );
igraph_vector_resize(&v, 3);
VECTOR(v)[0]=0; VECTOR(v)[1]=4; VECTOR(v)[2]=2;
igraph_matrix_select_cols(&m, &m2, &v);
print_matrix(&m2);
igraph_vector_resize(&v, 1);
igraph_matrix_select_cols(&m, &m2, &v);
print_matrix(&m2);
igraph_vector_clear(&v);
igraph_matrix_select_cols(&m, &m2, &v);
if (!igraph_matrix_empty(&m2)) { return 9; }
igraph_vector_destroy(&v);
igraph_matrix_destroy(&m2);
igraph_matrix_destroy(&m);
if (IGRAPH_FINALLY_STACK_SIZE() != 0) return 10;
return 0;
}
