matrix_repart.c

Go to the documentation of this file.

/*****************************************************************************
 * Copyright (c) 2019 FrontISTR Commons
 * This software is released under the MIT License, see LICENSE.txt
 *****************************************************************************/
 
#include <stdio.h>
#include <stdlib.h>
#include "hecmw_util.h"
 
#ifdef HECMW_WITH_METIS
#include "metis.h"
#endif
 
#include "matrix2graph.h"
#include "separator.h"
Separator_result *separator;
 
void bi_part_directive(int *neqns, int *nttbr, int *irow, int *jcol,
                       int *num_graph1, int *num_graph2, int *num_separator) {
#ifdef HECMW_WITH_METIS
  graph_type *graph;
  int num_of_row, num_of_col, num_of_nzero;
 
  int *perm, *iperm;
 
  num_of_row   = *neqns;
  num_of_col   = *neqns;
  num_of_nzero = *nttbr;
 
  graph = (graph_type *)malloc(sizeof(graph_type));
  if (graph == NULL) separator_memory_exit("graph");
 
  fprintf(stderr, "Start transforming matrix to graph\n");
  matrix2graph(num_of_row, num_of_col, num_of_nzero, irow, jcol, graph);
  fprintf(stderr,
          "Graph Information "
          "---------------------------------------------------\n");
  fprintf(stderr, "#Vertices: %d, #Edges: %d\n\n", graph->nvtxs,
          graph->nedges / 2);
  perm  = (int *)calloc(num_of_col, sizeof(int));
  iperm = (int *)calloc(num_of_col, sizeof(int));
  if ((perm == NULL) || (iperm == NULL))
    separator_memory_exit("matrix_repart: perm, iperm");
 
#if defined(METIS_VER_MAJOR) && (METIS_VER_MAJOR == 5)
  {
    int options[METIS_NOPTIONS];
    METIS_SetDefaultOptions(options);
    options[METIS_OPTION_COMPRESS] = 0;
    METIS_NodeND(&graph->nvtxs, graph->xadj, graph->adjncy, NULL, NULL, perm,
                 iperm);
  }
#else
  {
    int options[8];
    int num_flag;
    /* the following are options. see METIS manual for METIS_NODEND() */
    options[0] = 1; /* specify parameters */
    options[1] = 3; /* default */
    options[2] = 1; /* default */
    options[3] = 2; /* default */
    options[4] = 0; /* default */
    options[5] = 0; /* do not try to compress the matrix */
    options[6] = 0; /* default */
    options[7] = 1; /* default */
 
    num_flag = 0;
    METIS_NodeND(&graph->nvtxs, graph->xadj, graph->adjncy, &num_flag, options,
                 perm, iperm);
  }
#endif
  /* copy to separator */
 
  *num_graph1    = separator->num_of_lgraph;
  *num_graph2    = separator->num_of_rgraph;
  *num_separator = separator->num_of_separator;
 
#else
  fprintf(
      stderr,
      "Error: Direct Parallel Solver not available. Please install Metis.\n");
  HECMW_abort(HECMW_comm_get_comm());
#endif
 
  return;
}
 
void bi_part_directive_(int *neqns, int *nttbr, int *irow, int *jcol,
                        int *num_graph1, int *num_graph2, int *num_separator) {
  bi_part_directive(neqns, nttbr, irow, jcol, num_graph1, num_graph2,
                    num_separator);
}
 
void bi_part_directive__(int *neqns, int *nttbr, int *irow, int *jcol,
                         int *num_graph1, int *num_graph2, int *num_separator) {
  bi_part_directive(neqns, nttbr, irow, jcol, num_graph1, num_graph2,
                    num_separator);
}
 
void BI_PART_DIRECTIVE(int *neqns, int *nttbr, int *irow, int *jcol,
                       int *num_graph1, int *num_graph2, int *num_separator) {
  bi_part_directive(neqns, nttbr, irow, jcol, num_graph1, num_graph2,
                    num_separator);
}