46#include "Epetra_Map.h"
47#include "Epetra_MultiVector.h"
49#include "Epetra_Comm.h"
50#include "Epetra_SerialComm.h"
53# include "Epetra_MpiComm.h"
56#include "BelosConfigDefs.hpp"
57#include "BelosMVOPTester.hpp"
58#include "BelosEpetraAdapter.hpp"
60#ifdef HAVE_EPETRA_THYRA
62# include "Thyra_EpetraThyraWrappers.hpp"
63# include "Thyra_EpetraLinearOp.hpp"
67main (
int argc,
char *argv[])
71 using Teuchos::rcp_implicit_cast;
76 MPI_Init (&argc, &argv);
77 Epetra_MpiComm comm (MPI_COMM_WORLD);
80 Epetra_SerialComm comm;
82 const int myRank = comm.MyPID ();
85 const int globalNumRows = 100;
86 const int blockSize = 3;
87 const int indexBase = 0;
89 RCP<const Epetra_Map> range_epetra (
new Epetra_Map (globalNumRows, indexBase, comm));
90 RCP<Epetra_MultiVector> X_epetra (
new Epetra_MultiVector (*range_epetra, blockSize));
97#ifdef HAVE_EPETRA_THYRA
99 RCP<const Thyra::VectorSpaceBase<double> > range_thyra =
100 Thyra::create_VectorSpace (range_epetra);
102 RCP<Thyra::MultiVectorBase<double> > X_thyra =
103 Thyra::create_MultiVector (rcp_implicit_cast<Epetra_MultiVector> (X_epetra), range_thyra);
118 std::cout <<
"End Result: TEST PASSED" << std::endl;
124 std::cout <<
"End Result: TEST FAILED" << std::endl;
int main(int argc, char *argv[])
Stub adaptor from Thyra::MultiVectorBase to TSQR.