16 #include "Teuchos_ParameterList.hpp" 17 #include "Teuchos_RCP.hpp" 18 #include "Galeri_ReadHB.h" 24 int main(
int argc,
char *argv[])
27 MPI_Init(&argc, &argv);
50 std::string matrix_file=
"bcsstm05_ns.rua";
53 Galeri::ReadHB(matrix_file.c_str(), Comm, readMap,
54 readA, readx, readb, readxexact);
58 std::cout <<
"Caught exception, maybe file name is incorrect" << std::endl;
85 assert (OriginalMap.
SameAs(*readMap));
95 xexact.
Export(*readxexact, exporter,
Add);
118 std::string solverType =
"Klu";
119 Teuchos::RCP<Amesos_BaseSolver> solver = Teuchos::rcp( factory.
Create(solverType, problem) );
122 int symRet = solver->SymbolicFactorization();
124 std::cout <<
"Processor "<< map.
Comm().
MyPID() <<
" : Symbolic factorization did not complete!" << std::endl;
128 int numRet = solver->NumericFactorization();
130 std::cout <<
"Processor "<< map.
Comm().
MyPID() <<
" : Numeric factorization did not complete!" << std::endl;
135 Comm.
MinAll( &numRet, &minRet, 1 );
139 std::cout << std::endl <<
"End Result: TEST FAILED" << std::endl;
143 std::cout << std::endl <<
"End Result: TEST PASSED" << std::endl;
bool SameAs(const Epetra_BlockMap &Map) const
int MinAll(double *PartialMins, double *GlobalMins, int Count) const
const Epetra_Map & RowMatrixRowMap() const
virtual int MyPID() const=0
int FillComplete(bool OptimizeDataStorage=true)
int Export(const Epetra_SrcDistObject &A, const Epetra_Import &Importer, Epetra_CombineMode CombineMode, const Epetra_OffsetIndex *Indexor=0)
int main(int argc, char *argv[])
bool GlobalIndicesInt() const
Factory for binding a third party direct solver to an Epetra_LinearProblem.
const Epetra_Comm & Comm() const
Amesos_BaseSolver * Create(const char *ClassType, const Epetra_LinearProblem &LinearProblem)
Amesos Create method.
int NumGlobalElements() const
const int NumericallySingularMatrixError
bool GlobalIndicesLongLong() const