46 const int numvecs,
const int stride)
60 const std::vector<int>& index )
61 :
Epetra_MultiVector(CV, P_vec, &(const_cast<std::vector<int> &>(index))[0], index.size())
124 int numvecs = index.size();
126 std::vector<int> index2( numvecs );
127 for(
int i=0; i<numvecs; i++)
130 TEUCHOS_TEST_FOR_EXCEPTION( tmp_vec==NULL, std::invalid_argument,
"Anasazi::EpetraMultiVec::SetBlocks() cast of MultiVec<double> to EpetraMultiVec failed.");
132 temp_vec.
MvAddMv( 1.0, A_vec, 0.0, A_vec );
135 temp_vec.
MvAddMv( 1.0, A, 0.0, A );
152 TEUCHOS_TEST_FOR_EXCEPTION( A_vec==NULL, std::invalid_argument,
"Anasazi::EpetraMultiVec::SetBlocks() cast of MultiVec<double> to EpetraMultiVec failed.");
155 Multiply(
'N',
'N', alpha, *A_vec, B_Pvec, beta ) != 0,
156 EpetraMultiVecFailure,
"Anasazi::EpetraMultiVec::MvTimesMatAddMv() call to Epetra_MultiVec::Multiply() returned a nonzero value.");
169 TEUCHOS_TEST_FOR_EXCEPTION( A_vec==NULL, std::invalid_argument,
"Anasazi::EpetraMultiVec::MvAddMv() cast of MultiVec<double> to EpetraMultiVec failed.");
171 TEUCHOS_TEST_FOR_EXCEPTION( B_vec==NULL, std::invalid_argument,
"Anasazi::EpetraMultiVec::MvAddMv() cast of MultiVec<double> to EpetraMultiVec failed.");
174 Update( alpha, *A_vec, beta, *B_vec, 0.0 ) != 0,
175 EpetraMultiVecFailure,
"Anasazi::EpetraMultiVec::MvAddMv() call to Epetra_MultiVec::Update() returned a nonzero value.");
186 #ifdef HAVE_ANASAZI_EXPERIMENTAL
195 Epetra_MultiVector B_Pvec(Epetra_DataAccess::View, LocalMap, B.values(), B.stride(), B.numCols());
198 B_Pvec.Multiply(
'T',
'N', alpha, *A_vec, *
this, 0.0 ) != 0,
199 EpetraMultiVecFailure,
"Anasazi::EpetraMultiVec::MvTransMv() call to Epetra_MultiVec::Multiply() returned a nonzero value.");
210 #ifdef HAVE_ANASAZI_EXPERIMENTAL
216 TEUCHOS_TEST_FOR_EXCEPTION( A_vec==NULL, std::invalid_argument,
"Anasazi::EpetraMultiVec::MvDot() cast of MultiVec<double> to EpetraMultiVec failed.");
218 if (( (
int)b.size() >= A_vec->
NumVectors() ) ) {
220 this->
Dot( *A_vec, &b[0] ) != 0,
221 EpetraMultiVecFailure,
"Anasazi::EpetraMultiVec::MvDot() call to Epetra_MultiVec::Dot() returned a nonzero value.");
235 "Anasazi::EpetraMultiVec::MvScale() alpha argument size was inconsistent with number of vectors in mv.");
237 std::vector<int> tmp_index( 1, 0 );
238 for (
int i=0; i<numvecs; i++) {
241 temp_vec.
Scale( alpha[i] ) != 0,
242 EpetraMultiVecFailure,
"Anasazi::EpetraMultiVec::MvScale() call to Epetra_MultiVec::Scale() returned a nonzero value.");
277 TEUCHOS_TEST_FOR_EXCEPTION( vec_X==NULL, std::invalid_argument,
"Anasazi::EpetraOp::Apply() cast of MultiVec<double> to Epetra_MultiVector failed.");
278 TEUCHOS_TEST_FOR_EXCEPTION( vec_Y==NULL, std::invalid_argument,
"Anasazi::EpetraOp::Apply() cast of MultiVec<double> to Epetra_MultiVector failed.");
280 int info = Epetra_Op->
Apply( *vec_X, *vec_Y );
282 "Anasazi::EpetraOp::Apply(): Error returned from Epetra_Operator::Apply()" );
298 : isAInverse( isAInverse_ ), Epetra_AOp(AOp), Epetra_MOp(MOp)
319 TEUCHOS_TEST_FOR_EXCEPTION( vec_X==NULL, std::invalid_argument,
"Anasazi::EpetraGenOp::Apply() cast of MultiVec<double> to Epetra_MultiVector failed.");
320 TEUCHOS_TEST_FOR_EXCEPTION( vec_Y==NULL, std::invalid_argument,
"Anasazi::EpetraGenOp::Apply() cast of MultiVec<double> to Epetra_MultiVector failed.");
326 info = Epetra_MOp->
Apply( *vec_X, temp_Y );
328 "Anasazi::EpetraGenOp::Apply(): Error returned from Epetra_Operator::Apply()" );
334 info = Epetra_AOp->
Apply( temp_Y, *vec_Y );
337 "Anasazi::EpetraGenOp::Apply(): Error returned from Epetra_Operator::Apply()" );
349 info = Epetra_MOp->
Apply( X, temp_Y );
350 if (info!=0)
return info;
356 info = Epetra_AOp->
Apply( temp_Y, Y );
371 info = Epetra_AOp->
Apply( X, temp_Y );
374 if (info!=0)
return info;
393 : Epetra_Op(Op), isTrans_(isTrans)
413 vec_X->NumVectors() );
415 TEUCHOS_TEST_FOR_EXCEPTION( vec_X==NULL , std::invalid_argument,
"Anasazi::EpetraSymOp::Apply() cast of MultiVec<double> to Epetra_MultiVector failed.");
416 TEUCHOS_TEST_FOR_EXCEPTION( vec_Y==NULL , std::invalid_argument,
"Anasazi::EpetraSymOp::Apply() cast of MultiVec<double> to Epetra_MultiVector failed.");
417 TEUCHOS_TEST_FOR_EXCEPTION( temp_vec==NULL, std::invalid_argument,
"Anasazi::EpetraSymOp::Apply() allocation Epetra_MultiVector failed.");
428 "Anasazi::EpetraSymOp::Apply(): Error returned from Epetra_Operator::Apply()" );
434 info=Epetra_Op->
Apply( *vec_X, *temp_vec );
438 "Anasazi::EpetraSymOp::Apply(): Error returned from Epetra_Operator::Apply()" );
446 "Anasazi::EpetraSymOp::Apply(): Error returned from Epetra_Operator::Apply()" );
450 info=Epetra_Op->
Apply( *temp_vec, *vec_Y );
454 "Anasazi::EpetraSymOp::Apply(): Error returned from Epetra_Operator::Apply()" );
461 "Anasazi::EpetraSymOp::Apply(): Error returned from Epetra_Operator::Apply()" );
474 if (info!=0) {
return info; }
479 info=Epetra_Op->
Apply( X, temp_vec );
480 if (info!=0) {
return info; }
484 if (info!=0) {
return info; }
487 info=Epetra_Op->
Apply( temp_vec, Y );
488 if (info!=0) {
return info; }
505 if (info!=0) {
return info; }
511 if (info!=0) {
return info; }
515 if (info!=0) {
return info; }
518 info=Epetra_Op->
Apply( temp_vec, Y );
519 if (info!=0) {
return info; }
537 : Epetra_MV(MV), isTrans_(isTrans)
560 info = temp_vec.
Multiply(
'T',
'N', 1.0, *Epetra_MV, *vec_X, 0.0 );
562 "Anasazi::EpetraSymMVOp::Apply(): Error returned from Epetra_MultiVector::Multiply()" );
565 info = vec_Y->
Multiply(
'N',
'N', 1.0, *Epetra_MV, temp_vec, 0.0 );
567 "Anasazi::EpetraSymMVOp::Apply(): Error returned from Epetra_MultiVector::Multiply()" );
574 info = temp_vec.
Multiply(
'N',
'N', 1.0, *Epetra_MV, *vec_X, 0.0 );
576 "Anasazi::EpetraSymMVOp::Apply(): Error returned from Epetra_MultiVector::Multiply()" );
579 info = vec_Y->
Multiply(
'T',
'N', 1.0, *Epetra_MV, temp_vec, 0.0 );
581 "Anasazi::EpetraSymMVOp::Apply(): Error returned from Epetra_MultiVector::Multiply()" );
596 : Epetra_MV(MV), Epetra_OP(OP)
600 Epetra_OP->
Apply( *Epetra_MV, *Epetra_WMV );
616 info = temp_vec.
Multiply(
'N',
'N', 1.0, *Epetra_WMV, *vec_X, 0.0 );
618 "Anasazi::EpetraWSymMVOp::Apply(): Error returned from Epetra_MultiVector::Multiply()" );
621 info = vec_Y->
Multiply(
'T',
'N', 1.0, *Epetra_MV, temp_vec, 0.0 );
623 "Anasazi::EpetraWSymMVOp::Apply(): Error returned from Epetra_MultiVector::Multiply()" );
637 : Epetra_MV(MV), Epetra_OP(OP)
641 Epetra_OP->
Apply( *Epetra_MV, *Epetra_WMV );
657 info = temp_vec.
Multiply(
'N',
'N', 1.0, *Epetra_WMV, *vec_X, 0.0 );
659 "Anasazi::EpetraW2SymMVOp::Apply(): Error returned from Epetra_MultiVector::Multiply()" );
662 info = vec_Y->
Multiply(
'T',
'N', 1.0, *Epetra_WMV, temp_vec, 0.0 );
664 "Anasazi::EpetraW2SymMVOp::Apply(): Error returned from Epetra_MultiVector::Multiply()" );
const Epetra_BlockMap & Map() const
int Dot(const Epetra_MultiVector &A, double *Result) const
EpetraMultiVecAccessor is an interfaceto allow any Anasazi::MultiVec implementation that is based on ...
void Apply(const MultiVec< double > &X, MultiVec< double > &Y) const
Apply method [inherited from Anasazi::Operator class].
void Apply(const MultiVec< double > &X, MultiVec< double > &Y) const
Apply method.
void MvTransMv(double alpha, const MultiVec< double > &A, Teuchos::SerialDenseMatrix< int, double > &B) const
Compute a dense matrix B through the matrix-matrix multiply .
ScalarType * values() const
MultiVec< double > * CloneViewNonConst(const std::vector< int > &index)
Creates a new EpetraMultiVec that shares the selected contents of *this.
void SetBlock(const MultiVec< double > &A, const std::vector< int > &index)
Copy the vectors in A to a set of vectors in *this.
virtual int SetUseTranspose(bool UseTranspose)=0
int ApplyInverse(const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
Apply inverse method [inherited from Epetra_Operator class].
virtual int ApplyInverse(const Epetra_MultiVector &X, Epetra_MultiVector &Y) const=0
#define TEUCHOS_TEST_FOR_EXCEPTION(throw_exception_test, Exception, msg)
void Apply(const MultiVec< double > &X, MultiVec< double > &Y) const
Apply method [inherited from Anasazi::Operator class].
~EpetraGenOp()
Destructor.
virtual const Epetra_Map & OperatorDomainMap() const=0
Namespace Anasazi contains the classes, structs, enums and utilities used by the Anasazi package...
virtual int Apply(const Epetra_MultiVector &X, Epetra_MultiVector &Y) const=0
EpetraSymMVOp(const Teuchos::RCP< const Epetra_MultiVector > &MV, bool isTrans=false)
Basic constructor for applying operator [default] or .
virtual const Epetra_Map & OperatorRangeMap() const=0
int Scale(double ScalarValue)
void MvDot(const MultiVec< double > &A, std::vector< double > &b) const
Compute a vector b where the components are the individual dot-products, i.e. where A[i] is the i-th...
void MvTimesMatAddMv(double alpha, const MultiVec< double > &A, const Teuchos::SerialDenseMatrix< int, double > &B, double beta)
Update *this with .
ConjType
Enumerated types used to specify conjugation arguments.
const MultiVec< double > * CloneView(const std::vector< int > &index) const
Creates a new EpetraMultiVec that shares the selected contents of *this.
virtual int GetNumberVecs() const =0
The number of vectors (i.e., columns) in the multivector.
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
OrdinalType numRows() const
MultiVec< double > * CloneCopy() const
Creates a new EpetraMultiVec and copies contents of *this into the new vector (deep copy)...
EpetraSymOp(const Teuchos::RCP< Epetra_Operator > &Op, bool isTrans=false)
Basic constructor for applying operator [default] or .
int ApplyInverse(const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
Apply inverse method [inherited from Epetra_Operator class].
MultiVec< double > * Clone(const int numvecs) const
Creates a new empty EpetraMultiVec containing numvecs columns.
const Epetra_Map & OperatorDomainMap() const
Returns the Epetra_Map object associated with the domain of this operator.
const Epetra_Comm & Comm() const
void MvAddMv(double alpha, const MultiVec< double > &A, double beta, const MultiVec< double > &B)
Replace *this with .
~EpetraSymOp()
Destructor.
void MvScale(double alpha)
Scale each element of the vectors in *this with alpha.
void Apply(const MultiVec< double > &X, MultiVec< double > &Y) const
Apply method.
int Update(double ScalarA, const Epetra_MultiVector &A, double ScalarThis)
EpetraW2SymMVOp(const Teuchos::RCP< const Epetra_MultiVector > &MV, const Teuchos::RCP< Epetra_Operator > &OP)
Basic constructor for applying operator .
int Multiply(char TransA, char TransB, double ScalarAB, const Epetra_MultiVector &A, const Epetra_MultiVector &B, double ScalarThis)
EpetraWSymMVOp(const Teuchos::RCP< const Epetra_MultiVector > &MV, const Teuchos::RCP< Epetra_Operator > &OP)
Basic constructor for applying operator .
EpetraMultiVecFailure is thrown when a return value from an Epetra call on an Epetra_MultiVector is n...
EpetraGenOp(const Teuchos::RCP< Epetra_Operator > &AOp, const Teuchos::RCP< Epetra_Operator > &MOp, bool isAInverse=true)
Basic constructor for applying operator [default] or .
void Apply(const MultiVec< double > &X, MultiVec< double > &Y) const
This method takes the Anasazi::MultiVec X and applies the operator to it resulting in the Anasazi::Mu...
OrdinalType stride() const
void Apply(const MultiVec< double > &X, MultiVec< double > &Y) const
Apply method.
OrdinalType numCols() const
EpetraMultiVec(const Epetra_BlockMap &Map_in, const int numvecs)
Basic EpetraMultiVec constructor.
EpetraOp(const Teuchos::RCP< Epetra_Operator > &Op)
Basic constructor. Accepts reference-counted pointer to an Epetra_Operator.
Declarations of Anasazi multi-vector and operator classes using Epetra_MultiVector and Epetra_Operato...
Exceptions thrown to signal error in operator application.
const Epetra_Map & OperatorDomainMap() const
Returns the Epetra_Map object associated with the domain of this operator.
Basic adapter class for Anasazi::MultiVec that uses Epetra_MultiVector.