60int main(
int argc,
char *argv[])
64 MPI_Init(&argc,&argv);
70 Teuchos::ParameterList GaleriList;
74 GaleriList.set(
"n", nx * nx);
75 GaleriList.set(
"nx", nx);
76 GaleriList.set(
"ny", nx);
77 Teuchos::RefCountPtr<Epetra_Map> Map = Teuchos::rcp( Galeri::CreateMap(
"Linear", Comm, GaleriList) );
78 Teuchos::RefCountPtr<Epetra_RowMatrix> A = Teuchos::rcp( Galeri::CreateCrsMatrix(
"Laplace2D", &*Map, GaleriList) );
84 Teuchos::ParameterList List;
92 std::string PrecType =
"ILU";
96 Teuchos::RefCountPtr<Ifpack_Preconditioner> Prec = Teuchos::rcp( Factory.Create(PrecType, &*A, OverlapLevel) );
97 assert(Prec != Teuchos::null);
100 List.set(
"fact: drop tolerance", 1e-9);
101 List.set(
"fact: level-of-fill", 1);
105 List.set(
"schwarz: combine mode",
"Add");
143 Solver.SetAztecOption(AZ_solver,AZ_gmres);
144 Solver.SetAztecOption(AZ_output,32);
147 Solver.SetPrecOperator(&*Prec);
150 Solver.Iterate(1550,1e-8);
158 return(EXIT_SUCCESS);