25int main(
int argc,
char *argv[])
28 MPI_Init(&argc, &argv);
29 Epetra_MpiComm Comm(MPI_COMM_WORLD);
31 Epetra_SerialComm Comm;
34 Teuchos::ParameterList GaleriList;
35 GaleriList.set(
"n", 5);
37 Epetra_Map* Map = CreateMap(
"Random", Comm, GaleriList);
42 Epetra_MultiVector x(*Map,NumVectors);
43 Epetra_MultiVector x_exact(*Map,NumVectors);
44 Epetra_MultiVector b(*Map,NumVectors);
52 Epetra_LinearProblem Problem(&A, &x, &b);
59 double norm = ComputeNorm(Matrix, &x_exact, &b);
60 if (Comm.MyPID() == 0)
61 std::cout <<
"norm = " << norm << std::endl;
int CreateCrsMatrix(const char *in_filename, const Epetra_Comm &Comm, Epetra_Map *&readMap, const bool transpose, const bool distribute, bool &symmetric, Epetra_CrsMatrix *&Matrix)
int main(int argc, char *argv[])