5 #include "Epetra_MpiComm.h" 7 #include "Epetra_SerialComm.h" 9 #include "Epetra_Map.h" 10 #include "Epetra_Vector.h" 11 #include "Epetra_Util.h" 15 #include "Galeri_Maps.h" 16 #include "Galeri_CrsMatrices.h" 17 #include "Galeri_Utils.h" 25 int main(
int argc,
char *argv[])
28 MPI_Init(&argc, &argv);
29 Epetra_MpiComm
Comm(MPI_COMM_WORLD);
31 Epetra_SerialComm
Comm;
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);
46 A.Multiply(
false,x_exact,b);
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;
Amesos_Klu: A serial, unblocked code ideal for getting started and for very sparse matrices...
ParameterList & set(std::string const &name, T const &value, std::string const &docString="", RCP< const ParameterEntryValidator > const &validator=null)
Amesos_TestRowMatrix: a class to test Epetra_RowMatrix based codes.
#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[])