61int main(
int argc,
char *argv[])
66 MPI_Init(&argc,&argv);
72 Teuchos::ParameterList GaleriList;
76 GaleriList.set(
"nx", nx);
77 GaleriList.set(
"ny", nx * Comm.
NumProc());
78 GaleriList.set(
"mx", 1);
79 GaleriList.set(
"my", Comm.
NumProc());
80 Teuchos::RefCountPtr<Epetra_Map> Map = Teuchos::rcp( Galeri::CreateMap(
"Cartesian2D", Comm, GaleriList) );
81 Teuchos::RefCountPtr<Epetra_RowMatrix> A = Teuchos::rcp( Galeri::CreateCrsMatrix(
"Laplace2D", &*Map, GaleriList) );
87 Teuchos::ParameterList List;
96 Ifpack_AdditiveSchwarz<Ifpack_ICT> Prec(&*A);
100 List.set(
"fact: ict level-of-fill", 1.0);
102 List.set(
"fact: absolute threshold", 0.0);
103 List.set(
"fact: relative threshold", 1.0);
104 List.set(
"fact: relaxation value", 0.0);
109 List.set(
"schwarz: filter singletons",
true);
143 Solver.SetAztecOption(AZ_solver,AZ_cg_condnum);
144 Solver.SetAztecOption(AZ_output,32);
147 Solver.SetPrecOperator(&Prec);
152 Solver.Iterate(1550,1e-5);
161 return (EXIT_SUCCESS);