73int main(
int argc,
char *argv[]) {
75 int i, ierr=0, returnierr=0;
81 MPI_Init(&argc,&argv);
84 MPI_Comm_size(MPI_COMM_WORLD, &size);
85 MPI_Comm_rank(MPI_COMM_WORLD, &rank);
97 if (argc>1)
if (argv[1][0]==
'-' && argv[1][1]==
'v') verbose =
true;
107 int MyPID = Comm.
MyPID();
110 bool verbose1 = verbose;
111 if( verbose ) verbose = (MyPID==0);
118 if( verbose1 ) cout << Comm << endl << flush;
120 int NumMyElements = 3;
121 int NumGlobalElements = NumProc*NumMyElements;
124 Epetra_Map Map( NumGlobalElements, NumMyElements, 0, Comm );
125 if( verbose1 ) cout << Map << endl << flush;
129 int PIDFac = 10*MyPID;
130 int index = PIDFac+2;
138 if( verbose1 ) cout << Graph << endl << flush;
140 EpetraExt::Transform_Composite<Epetra_LinearProblem> CompTrans;
147 EpetraExt::SameTypeTransform<Epetra_LinearProblem> * RI_LPTrans = RI_Trans;
148 CompTrans.addTransform( RI_LPTrans );
151 EpetraExt::SameTypeTransform<Epetra_LinearProblem> *
153 CompTrans.addTransform( RCM_LPTrans );
157 EpetraExt::SameTypeTransform<Epetra_LinearProblem> *
159 CompTrans.addTransform( Overlap_LPTrans );
173 vector<double> valA(3);
174 valA[0]=0; valA[1]=1; valA[2]=2;
int main(int argc, char *argv[])