44 #include "Epetra_MpiComm.h" 46 #include "Epetra_SerialComm.h" 48 #include "Epetra_Map.h" 49 #include "Epetra_CrsMatrix.h" 54 #include "Teuchos_RefCountPtr.hpp" 57 int main(
int argc,
char *argv[])
63 MPI_Init(&argc,&argv);
64 Epetra_MpiComm Comm(MPI_COMM_WORLD);
66 Epetra_SerialComm Comm;
70 if (Comm.NumProc() != 1) {
74 if (Comm.MyPID() == 0)
75 cout <<
"Please run this test with one process only" << endl;
88 #if !defined(EPETRA_NO_32BIT_GLOBAL_INDICES) || !defined(EPETRA_NO_64BIT_GLOBAL_INDICES) 89 Epetra_Map Map(-1,NumPoints,0,Comm);
94 std::vector<int> Indices(NumPoints);
95 std::vector<double> Values(NumPoints);
97 Teuchos::RefCountPtr<Epetra_CrsMatrix> A = Teuchos::rcp(
new Epetra_CrsMatrix(Copy,Map,0) );
98 for (
int i = 0 ; i < NumPoints ; ++i) {
102 NumEntries = NumPoints;
103 for (
int j = 0 ; j < NumPoints ; ++j) {
116 #if !defined(EPETRA_NO_32BIT_GLOBAL_INDICES) || !defined(EPETRA_NO_64BIT_GLOBAL_INDICES) 117 A->InsertGlobalValues(i, NumEntries, &Values[0], &Indices[0]);
127 Teuchos::RefCountPtr<Ifpack_RCMReordering> Reorder = Teuchos::rcp(
new Ifpack_RCMReordering() );
143 return(EXIT_SUCCESS);
int main(int argc, char *argv[])
Ifpack_RCMReordering: reverse Cuthill-McKee reordering.
#define IFPACK_CHK_ERR(ifpack_err)
Ifpack_ReorderFilter: a class for light-weight reorder of local rows and columns of an Epetra_RowMatr...