EpetraExt Package Browser (Single Doxygen Collection) Development
test/SymmRCM/cxx_main.cpp
Go to the documentation of this file.
00001 //@HEADER
00002 // ***********************************************************************
00003 // 
00004 //     EpetraExt: Epetra Extended - Linear Algebra Services Package
00005 //                 Copyright (2001) Sandia Corporation
00006 // 
00007 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
00008 // license for use of this work by or on behalf of the U.S. Government.
00009 // 
00010 // This library is free software; you can redistribute it and/or modify
00011 // it under the terms of the GNU Lesser General Public License as
00012 // published by the Free Software Foundation; either version 2.1 of the
00013 // License, or (at your option) any later version.
00014 //  
00015 // This library is distributed in the hope that it will be useful, but
00016 // WITHOUT ANY WARRANTY; without even the implied warranty of
00017 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
00018 // Lesser General Public License for more details.
00019 //  
00020 // You should have received a copy of the GNU Lesser General Public
00021 // License along with this library; if not, write to the Free Software
00022 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
00023 // USA
00024 // Questions? Contact Michael A. Heroux (maherou@sandia.gov) 
00025 // 
00026 // ***********************************************************************
00027 //@HEADER
00028 
00029 // SymmRCM Test routine
00030 #include <Epetra_ConfigDefs.h>
00031 #include "EpetraExt_Version.h"
00032 
00033 #ifdef EPETRA_MPI
00034 #include "Epetra_MpiComm.h"
00035 #include <mpi.h>
00036 #endif
00037 #include "Epetra_SerialComm.h"
00038 #include "Epetra_Time.h"
00039 #include "Epetra_BlockMap.h"
00040 #include "Epetra_CrsGraph.h"
00041 #include "Epetra_CrsMatrix.h"
00042 #include "Epetra_Vector.h"
00043 
00044 #include "EpetraExt_SymmRCM_CrsGraph.h"
00045 #include "EpetraExt_LPTrans_From_GraphTrans.h"
00046 
00047 #include "../epetra_test_err.h"
00048 
00049 int main(int argc, char *argv[]) {
00050 
00051   int i, ierr=0, returnierr=0;
00052 
00053 #ifdef EPETRA_MPI
00054 
00055   // Initialize MPI
00056 
00057   MPI_Init(&argc,&argv);
00058   int size, rank; // Number of MPI processes, My process ID
00059 
00060   MPI_Comm_size(MPI_COMM_WORLD, &size);
00061   MPI_Comm_rank(MPI_COMM_WORLD, &rank);
00062 
00063 #else
00064 
00065   int size = 1; // Serial case (not using MPI)
00066   int rank = 0;
00067 
00068 #endif
00069 
00070   bool verbose = false;
00071 
00072   // Check if we should print results to standard out
00073   if (argc>1) if (argv[1][0]=='-' && argv[1][1]=='v') verbose = true;
00074 
00075 
00076 #ifdef EPETRA_MPI
00077   Epetra_MpiComm Comm(MPI_COMM_WORLD);
00078 #else
00079   Epetra_SerialComm Comm;
00080 #endif
00081   if (!verbose) Comm.SetTracebackMode(0); // This should shut down any error traceback reporting
00082 
00083   int MyPID = Comm.MyPID();
00084   int NumProc = Comm.NumProc();
00085 
00086   bool verbose1 = verbose;
00087 
00088   if (verbose) verbose = (MyPID==0);
00089 
00090   if (verbose)
00091     cout << EpetraExt::EpetraExt_Version() << endl << endl;
00092 
00093   Comm.Barrier();
00094 
00095   if (verbose1) cout << Comm << endl << flush;
00096                                                                                 
00097   Comm.Barrier();
00098 
00099   int NumMyElements = 10000;
00100   int NumGlobalElements = NumMyElements*NumProc+EPETRA_MIN(NumProc,3);
00101   if (MyPID < 3) NumMyElements++;
00102   int IndexBase = 0;
00103   int ElementSize = 7;
00104   bool DistributedGlobal = (NumGlobalElements>NumMyElements);
00105   
00106   Epetra_BlockMap * Map;
00107 
00108   cout << "NOT IMPLEMENTED YET!!!\n";
00109 
00110 #ifdef EPETRA_MPI
00111   MPI_Finalize();
00112 #endif
00113 
00114   return returnierr;
00115 }
00116 
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines