summaryrefslogtreecommitdiff
path: root/benchmarks/CUDA/DG/3rdParty/ParMetis-3.1/ParMETISLib/kmetis.c
diff options
context:
space:
mode:
Diffstat (limited to 'benchmarks/CUDA/DG/3rdParty/ParMetis-3.1/ParMETISLib/kmetis.c')
-rw-r--r--benchmarks/CUDA/DG/3rdParty/ParMetis-3.1/ParMETISLib/kmetis.c274
1 files changed, 0 insertions, 274 deletions
diff --git a/benchmarks/CUDA/DG/3rdParty/ParMetis-3.1/ParMETISLib/kmetis.c b/benchmarks/CUDA/DG/3rdParty/ParMetis-3.1/ParMETISLib/kmetis.c
deleted file mode 100644
index 2ad20e3..0000000
--- a/benchmarks/CUDA/DG/3rdParty/ParMetis-3.1/ParMETISLib/kmetis.c
+++ /dev/null
@@ -1,274 +0,0 @@
-/*
- * Copyright 1997, Regents of the University of Minnesota
- *
- * kmetis.c
- *
- * This is the entry point of Moc_PARMETIS_PartGraphKway
- *
- * Started 10/19/96
- * George
- *
- * $Id: kmetis.c,v 1.9 2003/07/31 16:27:27 karypis Exp $
- *
- */
-
-#include <parmetislib.h>
-
-
-/***********************************************************************************
-* This function is the entry point of the parallel k-way multilevel partitionioner.
-* This function assumes nothing about the graph distribution.
-* It is the general case.
-************************************************************************************/
-void ParMETIS_V3_PartKway(idxtype *vtxdist, idxtype *xadj, idxtype *adjncy, idxtype *vwgt,
- idxtype *adjwgt, int *wgtflag, int *numflag, int *ncon, int *nparts,
- float *tpwgts, float *ubvec, int *options, int *edgecut, idxtype *part,
- MPI_Comm *comm)
-{
- int h, i;
- int nvtxs = -1, npes, mype;
- CtrlType ctrl;
- WorkSpaceType wspace;
- GraphType *graph;
- float avg, maximb, *mytpwgts;
- int moptions[10];
- int seed, dbglvl = 0;
- int iwgtflag, inumflag, incon, inparts, ioptions[10];
- float *itpwgts, iubvec[MAXNCON];
-
- MPI_Comm_size(*comm, &npes);
- MPI_Comm_rank(*comm, &mype);
-
-
- /********************************/
- /* Try and take care bad inputs */
- /********************************/
- if (options != NULL && options[0] == 1)
- dbglvl = options[PMV3_OPTION_DBGLVL];
-
- CheckInputs(STATIC_PARTITION, npes, dbglvl, wgtflag, &iwgtflag, numflag, &inumflag, ncon,
- &incon, nparts, &inparts, tpwgts, &itpwgts, ubvec, iubvec, NULL, NULL,
- options, ioptions, part, comm);
-
-
- /*********************************/
- /* Take care the nparts = 1 case */
- /*********************************/
- if (inparts <= 1) {
- idxset(vtxdist[mype+1]-vtxdist[mype], 0, part);
- *edgecut = 0;
- return;
- }
-
- /******************************/
- /* Take care of npes = 1 case */
- /******************************/
- if (npes == 1 && inparts > 1) {
- moptions[0] = 0;
- nvtxs = vtxdist[1];
-
- if (incon == 1) {
- METIS_WPartGraphKway(&nvtxs, xadj, adjncy, vwgt, adjwgt, &iwgtflag, &inumflag,
- &inparts, itpwgts, moptions, edgecut, part);
- }
- else {
- /* ADD: this is because METIS does not support tpwgts for all constraints */
- mytpwgts = fmalloc(inparts, "mytpwgts");
- for (i=0; i<inparts; i++)
- mytpwgts[i] = itpwgts[i*incon];
-
- moptions[7] = -1;
- METIS_mCPartGraphRecursive2(&nvtxs, &incon, xadj, adjncy, vwgt, adjwgt, &iwgtflag,
- &inumflag, &inparts, mytpwgts, moptions, edgecut, part);
-
- free(mytpwgts);
- }
-
- return;
- }
-
-
- if (inumflag == 1)
- ChangeNumbering(vtxdist, xadj, adjncy, part, npes, mype, 1);
-
- /*****************************/
- /* Set up control structures */
- /*****************************/
- if (ioptions[0] == 1) {
- dbglvl = ioptions[PMV3_OPTION_DBGLVL];
- seed = ioptions[PMV3_OPTION_SEED];
- }
- else {
- dbglvl = GLOBAL_DBGLVL;
- seed = GLOBAL_SEED;
- }
- SetUpCtrl(&ctrl, inparts, dbglvl, *comm);
- ctrl.CoarsenTo = amin(vtxdist[npes]+1, 25*incon*amax(npes, inparts));
- ctrl.seed = (seed == 0) ? mype : seed*mype;
- ctrl.sync = GlobalSEMax(&ctrl, seed);
- ctrl.partType = STATIC_PARTITION;
- ctrl.ps_relation = -1;
- ctrl.tpwgts = itpwgts;
- scopy(incon, iubvec, ctrl.ubvec);
-
- graph = Moc_SetUpGraph(&ctrl, incon, vtxdist, xadj, vwgt, adjncy, adjwgt, &iwgtflag);
-
- PreAllocateMemory(&ctrl, graph, &wspace);
-
- IFSET(ctrl.dbglvl, DBG_TIME, InitTimers(&ctrl));
- IFSET(ctrl.dbglvl, DBG_TIME, MPI_Barrier(ctrl.gcomm));
- IFSET(ctrl.dbglvl, DBG_TIME, starttimer(ctrl.TotalTmr));
-
- /*******************************************/
- /* Check for funny cases */
- /* -graph with no edges */
- /* -graph with self edges */
- /* -graph with poor vertex distribution */
- /* -graph with less than 2*npe nodes */
- /*******************************************/
- if (vtxdist[npes] < SMALLGRAPH || vtxdist[npes] < npes*20 || GlobalSESum(&ctrl, graph->nedges) == 0) {
- IFSET(ctrl.dbglvl, DBG_INFO, rprintf(&ctrl, "Partitioning a graph of size %d serially\n", vtxdist[npes]));
- PartitionSmallGraph(&ctrl, graph, &wspace);
- }
- else {
- /***********************/
- /* Partition the graph */
- /***********************/
- Moc_Global_Partition(&ctrl, graph, &wspace);
- ParallelReMapGraph(&ctrl, graph, &wspace);
- }
-
- IFSET(ctrl.dbglvl, DBG_TIME, MPI_Barrier(ctrl.gcomm));
- IFSET(ctrl.dbglvl, DBG_TIME, stoptimer(ctrl.TotalTmr));
-
- idxcopy(graph->nvtxs, graph->where, part);
- *edgecut = graph->mincut;
-
- /*******************/
- /* Print out stats */
- /*******************/
- IFSET(ctrl.dbglvl, DBG_TIME, PrintTimingInfo(&ctrl));
- IFSET(ctrl.dbglvl, DBG_TIME, MPI_Barrier(ctrl.gcomm));
-
- if (ctrl.dbglvl&DBG_INFO) {
- rprintf(&ctrl, "Final %d-way CUT: %6d \tBalance: ", inparts, graph->mincut);
- avg = 0.0;
- for (h=0; h<incon; h++) {
- maximb = 0.0;
- for (i=0; i<inparts; i++)
- maximb = amax(maximb, graph->gnpwgts[i*incon+h]/itpwgts[i*incon+h]);
- avg += maximb;
- rprintf(&ctrl, "%.3f ", maximb);
- }
- rprintf(&ctrl, " avg: %.3f\n", avg/(float)incon);
- }
-
- GKfree((void **)&itpwgts, (void **)&graph->lnpwgts, (void **)&graph->gnpwgts, (void **)&graph->nvwgt, LTERM);
- FreeInitialGraphAndRemap(graph, iwgtflag);
- FreeWSpace(&wspace);
- FreeCtrl(&ctrl);
-
- if (inumflag == 1)
- ChangeNumbering(vtxdist, xadj, adjncy, part, npes, mype, 0);
-
-}
-
-
-
-/*************************************************************************
-* This function is the driver to the multi-constraint partitioning algorithm.
-**************************************************************************/
-void Moc_Global_Partition(CtrlType *ctrl, GraphType *graph, WorkSpaceType *wspace)
-{
- int i, ncon, nparts;
- float ftmp, ubavg, lbavg, lbvec[MAXNCON];
-
- ncon = graph->ncon;
- nparts = ctrl->nparts;
- ubavg = savg(graph->ncon, ctrl->ubvec);
-
- SetUp(ctrl, graph, wspace);
-
- if (ctrl->dbglvl&DBG_PROGRESS) {
- rprintf(ctrl, "[%6d %8d %5d %5d] [%d] [", graph->gnvtxs, GlobalSESum(ctrl, graph->nedges),
- GlobalSEMin(ctrl, graph->nvtxs), GlobalSEMax(ctrl, graph->nvtxs), ctrl->CoarsenTo);
- for (i=0; i<ncon; i++)
- rprintf(ctrl, " %.3f", GlobalSEMinFloat(ctrl,graph->nvwgt[samin_strd(graph->nvtxs, graph->nvwgt+i, ncon)*ncon+i]));
- rprintf(ctrl, "] [");
- for (i=0; i<ncon; i++)
- rprintf(ctrl, " %.3f", GlobalSEMaxFloat(ctrl, graph->nvwgt[samax_strd(graph->nvtxs, graph->nvwgt+i, ncon)*ncon+i]));
- rprintf(ctrl, "]\n");
- }
-
- if (graph->gnvtxs < 1.3*ctrl->CoarsenTo ||
- (graph->finer != NULL &&
- graph->gnvtxs > graph->finer->gnvtxs*COARSEN_FRACTION)) {
-
- /* Done with coarsening. Find a partition */
- graph->where = idxmalloc(graph->nvtxs+graph->nrecv, "graph->where");
- Moc_InitPartition_RB(ctrl, graph, wspace);
-
- if (ctrl->dbglvl&DBG_PROGRESS) {
- Moc_ComputeParallelBalance(ctrl, graph, graph->where, lbvec);
- rprintf(ctrl, "nvtxs: %10d, balance: ", graph->gnvtxs);
- for (i=0; i<graph->ncon; i++)
- rprintf(ctrl, "%.3f ", lbvec[i]);
- rprintf(ctrl, "\n");
- }
-
- /* In case no coarsening took place */
- if (graph->finer == NULL) {
- Moc_ComputePartitionParams(ctrl, graph, wspace);
- Moc_KWayFM(ctrl, graph, wspace, NGR_PASSES);
- }
- }
- else {
- Moc_GlobalMatch_Balance(ctrl, graph, wspace);
-
- Moc_Global_Partition(ctrl, graph->coarser, wspace);
-
- Moc_ProjectPartition(ctrl, graph, wspace);
- Moc_ComputePartitionParams(ctrl, graph, wspace);
-
- if (graph->ncon > 1 && graph->level < 3) {
- for (i=0; i<ncon; i++) {
- ftmp = ssum_strd(nparts, graph->gnpwgts+i, ncon);
- if (ftmp != 0.0)
- lbvec[i] = (float)(nparts) *
- graph->gnpwgts[samax_strd(nparts, graph->gnpwgts+i, ncon)*ncon+i]/ftmp;
- else
- lbvec[i] = 1.0;
- }
- lbavg = savg(graph->ncon, lbvec);
-
- if (lbavg > ubavg + 0.035) {
- if (ctrl->dbglvl&DBG_PROGRESS) {
- Moc_ComputeParallelBalance(ctrl, graph, graph->where, lbvec);
- rprintf(ctrl, "nvtxs: %10d, cut: %8d, balance: ", graph->gnvtxs, graph->mincut);
- for (i=0; i<graph->ncon; i++)
- rprintf(ctrl, "%.3f ", lbvec[i]);
- rprintf(ctrl, "\n");
- }
-
- Moc_KWayBalance(ctrl, graph, wspace, graph->ncon);
- }
- }
-
- Moc_KWayFM(ctrl, graph, wspace, NGR_PASSES);
-
- if (ctrl->dbglvl&DBG_PROGRESS) {
- Moc_ComputeParallelBalance(ctrl, graph, graph->where, lbvec);
- rprintf(ctrl, "nvtxs: %10d, cut: %8d, balance: ", graph->gnvtxs, graph->mincut);
- for (i=0; i<graph->ncon; i++)
- rprintf(ctrl, "%.3f ", lbvec[i]);
- rprintf(ctrl, "\n");
- }
-
- if (graph->level != 0)
- GKfree((void **)&graph->lnpwgts, (void **)&graph->gnpwgts, LTERM);
- }
-
- return;
-}
-
-
generated by cgit v1.3 (git 2.53.0) at 2026-07-10 23:54:21 +0000