63int main(
int argc,
char* argv[])
66 MPI_Init(&argc, &argv);
81 Teuchos::ParameterList GaleriList;
83 GaleriList.set(
"nx", 10 * Comm.NumProc());
84 GaleriList.set(
"ny", 10);
86 GaleriList.set(
"mx", Comm.NumProc());
87 GaleriList.set(
"my", 1);
91 Map = CreateMap(
"Cartesian2D", Comm, GaleriList);
92 Matrix = CreateCrsMatrix(
"Laplace2D", Map, GaleriList);
97 Matrix->
Multiply(
false, ExactSolution, RHS);
109 double ResidualNorm = ComputeNorm(Matrix, &LHS, &RHS);
111 if (Comm.MyPID() == 0)
112 cout << ResidualNorm << endl;
117 catch (Galeri::Exception& rhs)
119 if (Comm.MyPID() == 0)
128 return(EXIT_SUCCESS);
int main(int argc, char *argv[])