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);
108 typedef Thyra::TsqrAdaptor<double> tsqr_adapter_type;
118 std::cout <<
"End Result: TEST PASSED" << std::endl;
124 std::cout <<
"End Result: TEST FAILED" << std::endl;