16 typedef ScalarType NormType;
18 std::vector<std::pair<std::string, std::string> > solvers;
19 solvers.push_back (std::make_pair (
"A",
"1"));
20 solvers.push_back (std::make_pair (
"A",
"2"));
21 solvers.push_back (std::make_pair (
"B",
"3"));
22 solvers.push_back (std::make_pair (
"B",
"4"));
23 solvers.push_back (std::make_pair (
"C",
"5"));
24 solvers.push_back (std::make_pair (
"C",
"6"));
26 for (
size_t k = 0; k < solvers.size (); ++k) {
27 const std::string packageName = solvers[k].first;
28 const std::string solverName = solvers[k].second;
29 cout <<
"Package \"" << packageName <<
"\", solver \"" << solverName
32 Trilinos::Details::getLinearSolver<MV, OP, NormType> (packageName, solverName);
33 if (solver.get () == NULL) {
34 std::ostringstream err;
35 err <<
"Solver \"" << solvers[k].second <<
"\" from package \""
36 << solvers[k].first <<
"\" does not exist!";
37 throw std::logic_error (err.str ());
43 solver->setMatrix (
A);
51 cout <<
"Finished solver->solve(X, B)" << endl << endl;