/*
Function: symbolFaxDgraph
Unused
+*/
int symbolFaxDgraph(SCOTCH_Dgraph * dgraph,
MPI_Comm pastix_comm)
{
PASTIX_INT baseval;
PASTIX_INT vertglobalnbr;
PASTIX_INT vertlocalnbr;
PASTIX_INT vertmaxnbr;
PASTIX_INT vertgstnbr;
PASTIX_INT *vertlocaltab;
PASTIX_INT edgelocalnbr;
PASTIX_INT edgeglobalnbr;
PASTIX_INT edgemaxnbr;
PASTIX_INT *edgelocaltab;
PASTIX_INT *edgeghosttab;
SCOTCH_dgraphData (dgraph,
(SCOTCH_Num *)&baseval,
(SCOTCH_Num *)&vertglobalnbr,
(SCOTCH_Num *)&vertlocalnbr,
(SCOTCH_Num *)&vertmaxnbr,
(SCOTCH_Num *)&vertgstnbr,
(SCOTCH_Num **)&vertlocaltab,
NULL,
NULL,
NULL,
(SCOTCH_Num *)&edgelocalnbr,
(SCOTCH_Num *)&edgeglobalnbr,
(SCOTCH_Num *)&edgemaxnbr,
(SCOTCH_Num **)&edgelocaltab,
(SCOTCH_Num **)&edgeghosttab,
NULL,
&pastix_comm);
return EXIT_SUCCESS;
}
int
main (
int argc,
char * argv[])
{
MPI_Comm proccomm;
int procglbnbr; /* Number of processes sharing graph data */
int proclocnum; /* Number of this process */
SCOTCH_Num vertglbnbr;
SCOTCH_Num vertlocnbr;
SCOTCH_Num baseval;
SCOTCH_Dgraph grafdat;
SCOTCH_Dgraph coargrafdat;
SCOTCH_Num coarvertglbnbr;
SCOTCH_Num coarvertlocnbr;
double coarrat;
FILE * file;
#ifdef SCOTCH_PTHREAD
int thrdlvlreqval;
int thrdlvlproval;
#endif /* SCOTCH_PTHREAD */
int i;
errorProg (argv[0]);
#ifdef SCOTCH_PTHREAD
thrdlvlreqval = MPI_THREAD_MULTIPLE;
if (MPI_Init_thread (&argc, &argv, thrdlvlreqval, &thrdlvlproval) != MPI_SUCCESS)
errorPrint ("main: Cannot initialize (1)");
if (thrdlvlreqval > thrdlvlproval)
errorPrint ("main: MPI implementation is not thread-safe: recompile without SCOTCH_PTHREAD");
#else /* SCOTCH_PTHREAD */
if (MPI_Init (&argc, &argv) != MPI_SUCCESS)
errorPrint ("main: Cannot initialize (2)");
#endif /* SCOTCH_PTHREAD */
if (argc != 2) {
errorPrint ("main: invalid number of parameters");
exit (1);
}
proccomm = MPI_COMM_WORLD;
MPI_Comm_size (proccomm, &procglbnbr); /* Get communicator data */
MPI_Comm_rank (proccomm, &proclocnum);
fprintf (stderr, "Proc %2d of %2d, pid %d\n", proclocnum, procglbnbr, getpid ());
#ifndef SCOTCH_CHECK_AUTO
if (proclocnum == 0) { /* Synchronize on keybord input */
char c;
printf ("Waiting for key press...\n");
scanf ("%c", &c);
}
#endif /* SCOTCH_CHECK_AUTO */
if (MPI_Barrier (proccomm) != MPI_SUCCESS) { /* Synchronize for debug */
errorPrint ("main: cannot communicate");
return (1);
}
if (SCOTCH_dgraphInit (&grafdat, proccomm) != 0) { /* Initialize source graph */
errorPrint ("main: cannot initialize graph (1)");
return (1);
}
file = NULL;
if ((proclocnum == 0) &&
((file = fopen (argv[1], "r")) == NULL)) {
errorPrint ("main: cannot open graph file");
return (1);
}
if (SCOTCH_dgraphLoad (&grafdat, file, 0, 0) != 0) {
errorPrint ("main: cannot load graph");
return (1);
}
if (file != NULL)
fclose (file);
if (MPI_Barrier (proccomm) != MPI_SUCCESS) { /* Synchronize for debug */
errorPrint ("main: cannot communicate");
return (1);
}
SCOTCH_dgraphData (&grafdat, NULL, &vertglbnbr, &vertlocnbr, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL);
coarrat = 0.8; /* Lazy coarsening ratio */
for (i = 0; i < 3; i ++) { /* For all test cases */
SCOTCH_Num * multloctab;
SCOTCH_Num multlocsiz;
SCOTCH_Num foldval;
char * foldstr;
char * coarstr;
int procnum;
int o;
switch (i) {
case 0 :
multlocsiz = vertlocnbr;
foldval = SCOTCH_COARSENNONE;
foldstr = "Plain coarsening";
break;
case 1 :
multlocsiz = (SCOTCH_Num) (((double) vertglbnbr * coarrat) / (double) (procglbnbr / 2)) + 1;
foldval = SCOTCH_COARSENFOLD;
foldstr = "Folding";
break;
case 2 :
multlocsiz = (SCOTCH_Num) (((double) vertglbnbr * coarrat) / (double) (procglbnbr / 2)) + 1;
foldval = SCOTCH_COARSENFOLDDUP;
foldstr = "Folding with duplication";
break;
}
if (proclocnum == 0)
printf ("%s\n", foldstr);
if ((multloctab = malloc (multlocsiz * 2 * sizeof (SCOTCH_Num))) == NULL) {
errorPrint ("main: cannot allocate multinode array");
return (1);
}
if (SCOTCH_dgraphInit (&coargrafdat, proccomm) != 0) { /* Initialize band graph */
errorPrint ("main: cannot initialize graph (2)");
return (1);
}
o = SCOTCH_dgraphCoarsen (&grafdat, 0, coarrat, foldval, &coargrafdat, multloctab);
SCOTCH_dgraphData (&coargrafdat, NULL, &coarvertglbnbr, &coarvertlocnbr, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL);
for (procnum = 0; procnum < procglbnbr; procnum ++) {
switch (o) {
case 0 :
coarstr = "coarse graph created";
break;
case 1 :
coarstr = "graph could not be coarsened";
break;
case 2 :
coarstr = "folded graph not created here";
break;
case 3 :
coarstr = "cannot create coarse graph";
break;
}
if (procnum == proclocnum)
printf ("%d: %s (%ld / %ld / %ld)\n", procnum, coarstr, (long) multlocsiz, (long) coarvertlocnbr, (long) vertlocnbr);
MPI_Barrier (proccomm);
}
if (coarvertlocnbr > multlocsiz) {
errorPrint ("main: invalid local multinode array size");
return (1);
}
SCOTCH_dgraphExit (&coargrafdat);
free (multloctab);
}
SCOTCH_dgraphExit (&grafdat);
MPI_Finalize ();
exit (0);
}
int
main (
int argc,
char * argv[])
{
MPI_Comm proccomm;
int procglbnbr; /* Number of processes sharing graph data */
int proclocnum; /* Number of this process */
SCOTCH_Num vertglbnbr;
SCOTCH_Num vertlocnbr;
SCOTCH_Num * fronloctab;
SCOTCH_Num baseval;
SCOTCH_Dgraph grafdat;
SCOTCH_Dgraph bandgrafdat;
SCOTCH_Num bandvertglbnbr;
SCOTCH_Num bandvertlocnbr;
SCOTCH_Num * bandvlblloctab;
FILE * file;
int procnum;
#ifdef SCOTCH_PTHREAD
int thrdlvlreqval;
int thrdlvlproval;
#endif /* SCOTCH_PTHREAD */
SCOTCH_errorProg (argv[0]);
#ifdef SCOTCH_PTHREAD
thrdlvlreqval = MPI_THREAD_MULTIPLE;
if (MPI_Init_thread (&argc, &argv, thrdlvlreqval, &thrdlvlproval) != MPI_SUCCESS)
SCOTCH_errorPrint ("main: Cannot initialize (1)");
if (thrdlvlreqval > thrdlvlproval)
SCOTCH_errorPrint ("main: MPI implementation is not thread-safe: recompile without SCOTCH_PTHREAD");
#else /* SCOTCH_PTHREAD */
if (MPI_Init (&argc, &argv) != MPI_SUCCESS)
SCOTCH_errorPrint ("main: Cannot initialize (2)");
#endif /* SCOTCH_PTHREAD */
if (argc != 2) {
SCOTCH_errorPrint ("usage: %s graph_file", argv[0]);
exit (EXIT_FAILURE);
}
proccomm = MPI_COMM_WORLD;
MPI_Comm_size (proccomm, &procglbnbr); /* Get communicator data */
MPI_Comm_rank (proccomm, &proclocnum);
fprintf (stderr, "Proc %2d of %2d, pid %d\n", proclocnum, procglbnbr, getpid ());
#ifdef SCOTCH_CHECK_NOAUTO
if (proclocnum == 0) { /* Synchronize on keybord input */
char c;
printf ("Waiting for key press...\n");
scanf ("%c", &c);
}
#endif /* SCOTCH_CHECK_NOAUTO */
if (MPI_Barrier (proccomm) != MPI_SUCCESS) { /* Synchronize for debug */
SCOTCH_errorPrint ("main: cannot communicate (1)");
exit (EXIT_FAILURE);
}
if (SCOTCH_dgraphInit (&grafdat, proccomm) != 0) { /* Initialize source graph */
SCOTCH_errorPrint ("main: cannot initialize graph (1)");
exit (EXIT_FAILURE);
}
file = NULL;
if ((proclocnum == 0) &&
((file = fopen (argv[1], "r")) == NULL)) {
SCOTCH_errorPrint ("main: cannot open graph file");
exit (EXIT_FAILURE);
}
if (SCOTCH_dgraphLoad (&grafdat, file, 0, 0) != 0) {
SCOTCH_errorPrint ("main: cannot load graph");
exit (EXIT_FAILURE);
}
if (file != NULL)
fclose (file);
if (MPI_Barrier (proccomm) != MPI_SUCCESS) { /* Synchronize for debug */
SCOTCH_errorPrint ("main: cannot communicate (2)");
exit (EXIT_FAILURE);
}
SCOTCH_dgraphData (&grafdat, NULL, &vertglbnbr, &vertlocnbr, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL);
if ((fronloctab = malloc (vertlocnbr * sizeof (SCOTCH_Num))) == NULL) {
SCOTCH_errorPrint ("main: cannot allocate frontier array");
exit (EXIT_FAILURE);
}
if (SCOTCH_dgraphInit (&bandgrafdat, proccomm) != 0) { /* Initialize band graph */
SCOTCH_errorPrint ("main: cannot initialize graph (2)");
exit (EXIT_FAILURE);
}
fronloctab[0] = 0;
if (SCOTCH_dgraphBand (&grafdat, (proclocnum == 1) ? 1 : 0, fronloctab, 4, &bandgrafdat) != 0) {
SCOTCH_errorPrint ("main: cannot compute band graph");
exit (EXIT_FAILURE);
}
free (fronloctab);
SCOTCH_dgraphData (&bandgrafdat, &baseval, &bandvertglbnbr, &bandvertlocnbr, NULL, NULL, NULL, NULL, NULL, &bandvlblloctab, NULL, NULL, NULL, NULL, NULL, NULL, NULL);
for (procnum = 0; procnum < procglbnbr; procnum ++) {
SCOTCH_Num bandvertlocnum;
MPI_Barrier (proccomm);
if (procnum == proclocnum) {
if ((file = fopen ("/tmp/test_scotch_dgraph_band.map", (procnum == 0) ? "w" : "a+")) == NULL) {
SCOTCH_errorPrint ("main: cannot open mapping file");
exit (EXIT_FAILURE);
}
if (procnum == 0)
fprintf (file, "%ld\n", (long) bandvertglbnbr);
for (bandvertlocnum = 0; bandvertlocnum < bandvertlocnbr; bandvertlocnum ++)
fprintf (file, "%ld\t1\n", (long) bandvlblloctab[bandvertlocnum]);
fclose (file);
}
}
SCOTCH_dgraphExit (&bandgrafdat);
SCOTCH_dgraphExit (&grafdat);
MPI_Finalize ();
exit (EXIT_SUCCESS);
}