1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
|
/* checkComponents.c */
#include "../Graph.h"
#include "../../timings.h"
/*--------------------------------------------------------------------*/
int
main ( int argc, char *argv[] )
/*
--------------------------------------------------
read in a Graph and check for connected components
created -- 96fmar02, cca
--------------------------------------------------
*/
{
char *inGraphFileName ;
double t1, t2 ;
int icomp, msglvl, ncomp, rc, v ;
int *counts, *weights ;
IV *mapIV ;
Graph *g ;
FILE *msgFile ;
if ( argc != 4 ) {
fprintf(stdout,
"\n\n usage : %s msglvl msgFile inGraphFile"
"\n msglvl -- message level"
"\n msgFile -- message file"
"\n inGraphFile -- input file, must be *.graphf or *.graphb"
"\n", argv[0]) ;
return(0) ;
}
msglvl = atoi(argv[1]) ;
if ( strcmp(argv[2], "stdout") == 0 ) {
msgFile = stdout ;
} else if ( (msgFile = fopen(argv[2], "a")) == NULL ) {
fprintf(stderr, "\n fatal error in %s"
"\n unable to open file %s\n",
argv[0], argv[2]) ;
return(-1) ;
}
inGraphFileName = argv[3] ;
fprintf(msgFile,
"\n %s "
"\n msglvl -- %d"
"\n msgFile -- %s"
"\n inFile -- %s"
"\n",
argv[0], msglvl, argv[2], inGraphFileName) ;
fflush(msgFile) ;
/*
------------------------
read in the Graph object
------------------------
*/
if ( strcmp(inGraphFileName, "none") == 0 ) {
fprintf(msgFile, "\n no file to read from") ;
exit(0) ;
}
g = Graph_new() ;
MARKTIME(t1) ;
rc = Graph_readFromFile(g, inGraphFileName) ;
MARKTIME(t2) ;
fprintf(msgFile, "\n CPU %9.5f : read in graph from file %s",
t2 - t1, inGraphFileName) ;
if ( rc != 1 ) {
fprintf(msgFile, "\n return value %d from Graph_readFromFile(%p,%s)",
rc, g, inGraphFileName) ;
exit(-1) ;
}
fprintf(msgFile, "\n\n after reading Graph object from file %s",
inGraphFileName) ;
if ( msglvl > 2 ) {
Graph_writeForHumanEye(g, msgFile) ;
} else {
Graph_writeStats(g, msgFile) ;
}
fflush(msgFile) ;
/*
---------------------------------------
get the map from vertices to components
---------------------------------------
*/
MARKTIME(t1) ;
mapIV = Graph_componentMap(g) ;
MARKTIME(t2) ;
ncomp = 1 + IV_max(mapIV) ;
fprintf(msgFile, "\n CPU %9.5f : find %d components",
t2 - t1, ncomp) ;
if ( msglvl > 1 ) {
fprintf(msgFile, "\n\n component map") ;
IV_writeForHumanEye(mapIV, msgFile) ;
}
/*
------------------------------------
get the statistics on the components
------------------------------------
*/
counts = IVinit(ncomp, 0) ;
weights = IVinit(ncomp, 0) ;
MARKTIME(t1) ;
Graph_componentStats(g, IV_entries(mapIV), counts, weights) ;
MARKTIME(t2) ;
fprintf(msgFile, "\n CPU %9.5f : compute component statistics",
t2 - t1) ;
fprintf(msgFile, "\n\n component # vertices weight") ;
for ( icomp = 0 ; icomp < ncomp ; icomp++ ) {
fprintf(msgFile, "\n %7d %7d %7d",
icomp, counts[icomp], weights[icomp]) ;
}
/*
----------------
free the storage
----------------
*/
Graph_free(g) ;
IV_free(mapIV) ;
IVfree(counts) ;
IVfree(weights) ;
fprintf(msgFile, "\n") ;
fclose(msgFile) ;
return(1) ; }
/*--------------------------------------------------------------------*/
|
