50 #include "Epetra_MpiComm.h" 52 #include "Epetra_SerialComm.h" 54 #include "Epetra_CrsMatrix.h" 55 #include "Epetra_Vector.h" 56 #include "Epetra_LinearProblem.h" 57 #include "Epetra_Map.h" 58 #include "Galeri_Maps.h" 59 #include "Galeri_CrsMatrices.h" 60 #include "Teuchos_ParameterList.hpp" 61 #include "Teuchos_RefCountPtr.hpp" 65 #ifdef HAVE_IFPACK_AMESOS 76 int main(
int argc,
char *argv[])
82 MPI_Init(&argc,&argv);
83 Epetra_MpiComm Comm( MPI_COMM_WORLD );
85 Epetra_SerialComm Comm;
90 GaleriList.
set(
"n",
n);
91 Teuchos::RefCountPtr<Epetra_Map> Map =
Teuchos::rcp( Galeri::CreateMap64(
"Linear", Comm, GaleriList) );
92 Teuchos::RefCountPtr<Epetra_CrsMatrix>
A =
Teuchos::rcp( Galeri::CreateCrsMatrix(
"Minij", &*Map, GaleriList) );
95 Teuchos::RefCountPtr<Ifpack_Preconditioner> Prec;
98 assert (Prec != Teuchos::null);
104 assert (Prec != Teuchos::null);
110 assert (Prec != Teuchos::null);
116 assert (Prec != Teuchos::null);
122 assert (Prec != Teuchos::null);
128 assert (Prec != Teuchos::null);
134 assert (Prec != Teuchos::null);
140 assert (Prec != Teuchos::null);
146 assert (Prec != Teuchos::null);
152 assert (Prec != Teuchos::null);
158 assert (Prec != Teuchos::null);
164 assert (Prec != Teuchos::null);
169 #ifdef HAVE_IFPACK_AMESOS 171 assert (Prec != Teuchos::null);
177 assert (Prec != Teuchos::null);
184 assert (Prec != Teuchos::null);
190 assert (Prec != Teuchos::null);
196 assert (Prec != Teuchos::null);
205 #if defined (HAVE_IFPACK_SUPPORTGRAPH) && defined (HAVE_IFPACK_AMESOS) 207 assert (Prec !=Teuchos::null);
217 #ifdef HAVE_IFPACK_SUPPORTGRAPH 219 assert (Prec != Teuchos::null);
227 if (Comm.MyPID() == 0)
228 cout <<
"Test `PrecondititonerFactory_LL.exe' passed!" << endl;
236 return(EXIT_SUCCESS);
ParameterList & set(std::string const &name, T const &value, std::string const &docString="", RCP< const ParameterEntryValidator > const &validator=null)
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
static Ifpack_Preconditioner * Create(EPrecType PrecType, Epetra_RowMatrix *Matrix, const int overlap=0, bool overrideSerialDefault=false)
Creates an instance of Ifpack_Preconditioner given the enum value of the preconditioner type (can not...
#define IFPACK_CHK_ERR(ifpack_err)
int main(int argc, char *argv[])