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#ifndef GALERI_TEST_USE_LONGLONG_GO
92 Map = CreateMap(
"Cartesian2D", Comm, GaleriList);
94 Map = CreateMap64(
"Cartesian2D", Comm, GaleriList);
96 Matrix = CreateCrsMatrix(
"Laplace2D", Map, GaleriList);
101 Matrix->
Multiply(
false, ExactSolution, RHS);
113 double ResidualNorm = ComputeNorm(Matrix, &LHS, &RHS);
115 if (Comm.MyPID() == 0)
116 cout << ResidualNorm << endl;
121 catch (Galeri::Exception& rhs)
123 if (Comm.MyPID() == 0)
125 cerr <<
"Caught exception: ";
134 return(EXIT_SUCCESS);