/* 
  * Copyright 1997, Regents of the University of Minnesota 
  * 
  * oemetis.c 
  * 
  * This file contains the driving routine for multilevel method 
  * 
  * Started 8/28/94 
  * George 
  * 
  * $Id: oemetis.c,v 1.1 1998/11/27 17:59:37 karypis Exp $ 
  * 
  */ 

 #include <metis.h> 



 /************************************************************************* 
 * Let the game begin 
 **************************************************************************/ 
 main(long argc, char *argv[]) 
 { 
   long i, options[10]; 
   idxtype *perm, *iperm; 
   GraphType graph; 
   char filename[256]; 
   long numflag = 0, wgtflag; 
   timer TOTALTmr, METISTmr, IOTmr, SMBTmr; 


   if (argc != 2) { 
     printf("Usage: %s <GraphFile>\n",argv[0]); 
     exit(4); 
   } 
      
   strcpy(filename, argv[1]); 

   cleartimer(TOTALTmr); 
   cleartimer(METISTmr); 
   cleartimer(IOTmr); 
   cleartimer(SMBTmr); 

   starttimer(TOTALTmr); 
   starttimer(IOTmr); 
   ReadGraph(&graph, filename, &wgtflag); 
   if (graph.nvtxs <= 0) { 
     printf("Empty graph. Nothing to do.\n"); 
     exit(4); 
   } 
   if (graph.ncon != 1) { 
     printf("Ordering can only be applied to graphs with one constraint.\n"); 
     exit(4); 
   } 
   stoptimer(IOTmr); 

   /* Ordering does not use weights! */ 
   GKfree(&graph.vwgt, &graph.adjwgt, LTERM); 

   printf("**********************************************************************\n"); 
   printf("%s", METISTITLE); 
   printf("Graph Information ---------------------------------------------------\n"); 
   printf("  Name: %s, #Vertices: %ld, #Edges: %ld\n\n", filename, graph.nvtxs, graph.nedges/2); 
   printf("Edge-Based Ordering... ----------------------------------------------\n"); 

   perm = idxmalloc(graph.nvtxs, "main: perm"); 
   iperm = idxmalloc(graph.nvtxs, "main: iperm"); 
   options[0] = 0; 

   starttimer(METISTmr); 
   METIS_EdgeND(&graph.nvtxs, graph.xadj, graph.adjncy, &numflag, options, perm, iperm); 
   stoptimer(METISTmr); 

   starttimer(IOTmr); 
   WritePermutation(filename, iperm, graph.nvtxs);  
   stoptimer(IOTmr); 

   starttimer(SMBTmr); 
   ComputeFillIn(&graph, iperm); 
   stoptimer(SMBTmr); 

   stoptimer(TOTALTmr); 

   printf("\nTiming Information --------------------------------------------------\n"); 
   printf("  I/O:                     \t %7.3f\n", gettimer(IOTmr)); 
   printf("  Ordering:                \t %7.3f   (OEMETIS time)\n", gettimer(METISTmr)); 
   printf("  Symbolic Factorization:  \t %7.3f\n", gettimer(SMBTmr)); 
   printf("  Total:                   \t %7.3f\n", gettimer(TOTALTmr)); 
   printf("**********************************************************************\n"); 


   GKfree(&graph.xadj, &graph.adjncy, &perm, &iperm, LTERM); 
 }