45 #include "Galeri_Maps.h" 46 #include "Galeri_CrsMatrices.h" 47 #include "Galeri_Utils.h" 48 #include "Epetra_MultiVector.h" 50 #include "Teuchos_UnitTestHarness.hpp" 51 #include "Teuchos_RCP.hpp" 52 #include "Epetra_ConfigDefs.h" 53 #include "Epetra_Vector.h" 54 #include "Epetra_RowMatrix.h" 55 #include "Epetra_MultiVector.h" 56 #include "Epetra_CrsMatrix.h" 57 #include "Epetra_Map.h" 60 #include "Epetra_MpiComm.h" 62 #include "Epetra_SerialComm.h" 65 #include "Teuchos_Array.hpp" 79 Epetra_MpiComm comm(MPI_COMM_WORLD);
80 Epetra_Map RowMap(
N, 0, comm);
81 Epetra_CrsMatrix Matrix(Copy, RowMap, 1);
82 for(
int i = 0; i <
N; i++){
87 Matrix.InsertGlobalValues(i, 1, values, indices);
89 Matrix.FillComplete();
96 Epetra_MpiComm Comm(MPI_COMM_WORLD);
100 Epetra_CrsMatrix* Matrix;
102 Teuchos::ParameterList GaleriList;
104 int nx = 10 * Comm.NumProc();
105 int ny = 10 * Comm.NumProc();
106 GaleriList.set(
"nx", nx);
107 GaleriList.set(
"ny", ny);
113 Map = Galeri::CreateMap(
"Cartesian2D", Comm, GaleriList);
116 Matrix = Galeri::CreateCrsMatrix(
"Biharmonic2D", Map, GaleriList);
119 Ifpack_Hypre preconditioner(Matrix);
123 Epetra_MultiVector X(Matrix->RowMatrixRowMap(), 2);
124 Epetra_MultiVector KnownX(Matrix->RowMatrixRowMap(), 2);
126 Epetra_MultiVector B(Matrix->RowMatrixRowMap(), 2);
127 Matrix->Multiply(
false, KnownX, B);
133 Teuchos::ParameterList list(
"Preconditioner List");
154 preconditioner.SetParameters(list);
155 preconditioner.Compute();
160 catch (Galeri::Exception& rhs)
162 if (Comm.MyPID() == 0)
164 cerr <<
"Caught exception: ";
TEUCHOS_UNIT_TEST(Ifpack_Hypre, Construct)