48 #include "Galeri_Maps.h" 49 #include "Galeri_CrsMatrices.h" 91 int main(
int argc,
char *argv[])
95 MPI_Init(&argc, &argv);
115 ParameterList GaleriList;
116 GaleriList.set(
"nx", nx);
117 GaleriList.set(
"ny", ny);
118 GaleriList.set(
"mx", 1);
119 GaleriList.set(
"my", Comm.
NumProc());
121 Epetra_Map* Map = CreateMap(
"Cartesian2D", Comm, GaleriList);
160 std::string SolverType =
"Klu";
167 std::cerr <<
"Specified solver is not available" << std::endl;
173 return(EXIT_SUCCESS);
189 Teuchos::ParameterList List;
190 List.set(
"PrintTiming",
true);
191 List.set(
"PrintStatus",
true);
193 Solver->SetParameters(List);
202 if (Comm.
MyPID() == 0)
203 std::cout <<
"Starting symbolic factorization..." << std::endl;
204 Solver->SymbolicFactorization();
207 if (Comm.
MyPID() == 0)
208 std::cout <<
"Starting numeric factorization..." << std::endl;
209 Solver->NumericFactorization();
212 if (Comm.
MyPID() == 0)
213 std::cout <<
"Starting solution phase..." << std::endl;
264 return(EXIT_SUCCESS);
int PutScalar(double ScalarConstant)
Factory for binding a third party direct solver to an Epetra_LinearProblem.
int CreateCrsMatrix(const char *in_filename, const Epetra_Comm &Comm, Epetra_Map *&readMap, const bool transpose, const bool distribute, bool &symmetric, Epetra_CrsMatrix *&Matrix)
Amesos_BaseSolver * Create(const char *ClassType, const Epetra_LinearProblem &LinearProblem)
Amesos Create method.
int main(int argc, char *argv[])
Amesos_BaseSolver: A pure virtual class for direct solution of real-valued double-precision operators...