/*************************************************************************
* This function takes a graph and creates a sequence of coarser graphs
**************************************************************************/
GraphType *MCCoarsen2Way(CtrlType *ctrl, GraphType *graph)
{
int i, clevel;
GraphType *cgraph;
IFSET(ctrl->dbglvl, DBG_TIME, starttimer(ctrl->CoarsenTmr));
cgraph = graph;
clevel = 0;
do {
if (ctrl->dbglvl&DBG_COARSEN) {
printf("%6d %7d %10d [%d] [%6.4f", cgraph->nvtxs, cgraph->nedges,
idxsum(cgraph->nvtxs, cgraph->adjwgtsum), ctrl->CoarsenTo, ctrl->nmaxvwgt);
for (i=0; i<graph->ncon; i++)
printf(" %5.3f", ssum_strd(cgraph->nvtxs, cgraph->nvwgt+i, cgraph->ncon));
printf("]\n");
}
switch (ctrl->CType) {
case MATCH_RM:
MCMatch_RM(ctrl, cgraph);
break;
case MATCH_HEM:
if (clevel < 1)
MCMatch_RM(ctrl, cgraph);
else
MCMatch_HEM(ctrl, cgraph);
break;
case MATCH_SHEM:
if (clevel < 1)
MCMatch_RM(ctrl, cgraph);
else
MCMatch_SHEM(ctrl, cgraph);
break;
case MATCH_SHEMKWAY:
MCMatch_SHEM(ctrl, cgraph);
break;
case MATCH_SHEBM_ONENORM:
MCMatch_SHEBM(ctrl, cgraph, 1);
break;
case MATCH_SHEBM_INFNORM:
MCMatch_SHEBM(ctrl, cgraph, -1);
break;
case MATCH_SBHEM_ONENORM:
MCMatch_SBHEM(ctrl, cgraph, 1);
break;
case MATCH_SBHEM_INFNORM:
MCMatch_SBHEM(ctrl, cgraph, -1);
break;
default:
errexit("Unknown CType: %d\n", ctrl->CType);
}
cgraph = cgraph->coarser;
clevel++;
} while (cgraph->nvtxs > ctrl->CoarsenTo && cgraph->nvtxs < COARSEN_FRACTION2*cgraph->finer->nvtxs && cgraph->nedges > cgraph->nvtxs/2);
if (ctrl->dbglvl&DBG_COARSEN) {
printf("%6d %7d %10d [%d] [%6.4f", cgraph->nvtxs, cgraph->nedges,
idxsum(cgraph->nvtxs, cgraph->adjwgtsum), ctrl->CoarsenTo, ctrl->nmaxvwgt);
for (i=0; i<graph->ncon; i++)
printf(" %5.3f", ssum_strd(cgraph->nvtxs, cgraph->nvwgt+i, cgraph->ncon));
printf("]\n");
}
IFSET(ctrl->dbglvl, DBG_TIME, stoptimer(ctrl->CoarsenTmr));
return cgraph;
}
/*************************************************************************
* This function takes a graph and creates a sequence of coarser graphs
**************************************************************************/
GraphType *MCCoarsen2Way(CtrlType *ctrl, GraphType *graph)
{
idxtype i, clevel;
GraphType *cgraph;
IFSET(ctrl->dbglvl, DBG_TIME, gk_startcputimer(ctrl->CoarsenTmr));
cgraph = graph;
clevel = 0;
do {
if (ctrl->dbglvl&DBG_COARSEN) {
mprintf("%6D %7D %10D [%D] [%6.4f", cgraph->nvtxs, cgraph->nedges,
idxsum(cgraph->nvtxs, cgraph->adjwgtsum, 1), ctrl->CoarsenTo, ctrl->nmaxvwgt);
for (i=0; i<graph->ncon; i++)
mprintf(" %5.3f", gk_fsum(cgraph->nvtxs, cgraph->nvwgt+i, cgraph->ncon));
mprintf("]\n");
}
if (cgraph->nedges == 0) {
MCMatch_RM(ctrl, cgraph);
}
else {
switch (ctrl->CType) {
case MTYPE_RM:
MCMatch_RM(ctrl, cgraph);
break;
case MTYPE_HEM:
if (clevel < 1)
MCMatch_RM(ctrl, cgraph);
else
MCMatch_HEM(ctrl, cgraph);
break;
case MTYPE_SHEM:
if (clevel < 1)
MCMatch_RM(ctrl, cgraph);
else
MCMatch_SHEM(ctrl, cgraph);
break;
case MTYPE_SHEMKWAY:
MCMatch_SHEM(ctrl, cgraph);
break;
case MTYPE_SHEBM_ONENORM:
MCMatch_SHEBM(ctrl, cgraph, 1);
break;
case MTYPE_SHEBM_INFNORM:
MCMatch_SHEBM(ctrl, cgraph, -1);
break;
case MTYPE_SBHEM_ONENORM:
MCMatch_SBHEM(ctrl, cgraph, 1);
break;
case MTYPE_SBHEM_INFNORM:
MCMatch_SBHEM(ctrl, cgraph, -1);
break;
default:
errexit("Unknown CType: %d\n", ctrl->CType);
}
}
cgraph = cgraph->coarser;
clevel++;
} while (cgraph->nvtxs > ctrl->CoarsenTo && cgraph->nvtxs < COARSEN_FRACTION2*cgraph->finer->nvtxs && cgraph->nedges > cgraph->nvtxs/2);
if (ctrl->dbglvl&DBG_COARSEN) {
mprintf("%6D %7D %10D [%D] [%6.4f", cgraph->nvtxs, cgraph->nedges,
idxsum(cgraph->nvtxs, cgraph->adjwgtsum, 1), ctrl->CoarsenTo, ctrl->nmaxvwgt);
for (i=0; i<graph->ncon; i++)
mprintf(" %5.3f", gk_fsum(cgraph->nvtxs, cgraph->nvwgt+i, cgraph->ncon));
mprintf("]\n");
}
IFSET(ctrl->dbglvl, DBG_TIME, gk_stopcputimer(ctrl->CoarsenTmr));
return cgraph;
}
