42 #include "Galeri_Maps.h" 43 #include "Galeri_CrsMatrices.h" 44 #include "Galeri_Utils.h" 55 #include "Teuchos_ParameterList.hpp" 63 int 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);
int PutScalar(double ScalarConstant)
virtual int Multiply(bool TransA, const Epetra_MultiVector &X, Epetra_MultiVector &Y) const=0
int main(int argc, char *argv[])