5#include "Epetra_MpiComm.h"
7#include "Epetra_SerialComm.h"
10#include "Epetra_Vector.h"
11#include "Epetra_Util.h"
14#include "Teuchos_ParameterList.hpp"
15#include "Galeri_Maps.h"
16#include "Galeri_CrsMatrices.h"
17#include "Galeri_Utils.h"
19using namespace Galeri;
25int main(
int argc,
char *argv[])
28 MPI_Init(&argc, &argv);
29 Epetra_MpiComm Comm(MPI_COMM_WORLD);
31 Epetra_SerialComm Comm;
34 Teuchos::ParameterList GaleriList;
35 GaleriList.set(
"n", 5);
37 Epetra_Map* Map = CreateMap(
"Random", Comm, GaleriList);
42 Epetra_MultiVector x(*Map,NumVectors);
43 Epetra_MultiVector x_exact(*Map,NumVectors);
44 Epetra_MultiVector b(*Map,NumVectors);
52 Epetra_LinearProblem Problem(&A, &x, &b);
59 double norm = ComputeNorm(Matrix, &x_exact, &b);
60 if (Comm.MyPID() == 0)
61 std::cout <<
"norm = " << norm << std::endl;
#define AMESOS_CHK_ERR(a)
int CreateCrsMatrix(const char *in_filename, const Epetra_Comm &Comm, Epetra_Map *&readMap, const bool transpose, const bool distribute, bool &symmetric, Epetra_CrsMatrix *&Matrix)
int main(int argc, char *argv[])
Amesos_Klu: A serial, unblocked code ideal for getting started and for very sparse matrices,...
int SymbolicFactorization()
Performs SymbolicFactorization on the matrix A.
int NumericFactorization()
Performs NumericFactorization on the matrix A.
int Solve()
Solves A X = B (or AT x = B)
Amesos_TestRowMatrix: a class to test Epetra_RowMatrix based codes.
virtual int Multiply(bool TransA, const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
Returns the result of a Epetra_RowMatrix multiplied by a Epetra_MultiVector X in Y.