EpetraExt  Development
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
GenSQP_YUEpetraVector.cpp
Go to the documentation of this file.
1 /*
2 //@HEADER
3 // ***********************************************************************
4 //
5 // EpetraExt: Epetra Extended - Linear Algebra Services Package
6 // Copyright (2011) Sandia Corporation
7 //
8 // Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
9 // the U.S. Government retains certain rights in this software.
10 //
11 // Redistribution and use in source and binary forms, with or without
12 // modification, are permitted provided that the following conditions are
13 // met:
14 //
15 // 1. Redistributions of source code must retain the above copyright
16 // notice, this list of conditions and the following disclaimer.
17 //
18 // 2. Redistributions in binary form must reproduce the above copyright
19 // notice, this list of conditions and the following disclaimer in the
20 // documentation and/or other materials provided with the distribution.
21 //
22 // 3. Neither the name of the Corporation nor the names of the
23 // contributors may be used to endorse or promote products derived from
24 // this software without specific prior written permission.
25 //
26 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
27 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
28 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
29 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
30 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
31 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
32 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
33 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
34 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
35 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
36 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
37 //
38 // Questions? Contact Michael A. Heroux (maherou@sandia.gov)
39 //
40 // ***********************************************************************
41 //@HEADER
42 */
43 
45 
46 namespace GenSQP {
47 
48 YUEpetraVector::YUEpetraVector( const Teuchos::RCP<Epetra_MultiVector> &y_epetra_vec,
49  const Teuchos::RCP<Epetra_MultiVector> &u_epetra_vec )
50  :y_epetra_vec_(y_epetra_vec), u_epetra_vec_(u_epetra_vec)
51  {}
52 
53 // Overridden from Vector
54 
55 double YUEpetraVector::innerProd( const Vector &x ) const
56 {
57  double ydot[1];
58  double udot[1];
59  YUEpetraVector &ex = Teuchos::dyn_cast<YUEpetraVector>(const_cast <Vector&>(x));
60  y_epetra_vec_->Dot( *ex.y_epetra_vec_, ydot );
61  if (u_epetra_vec_.get() == 0)
62  udot[0] = 0.0;
63  else
64  u_epetra_vec_->Dot( *ex.u_epetra_vec_, udot );
65  return (ydot[0] + udot[0]);
66 }
67 
68 void YUEpetraVector::linComb( const double &alpha, const Vector &x, const double &beta )
69 {
70  YUEpetraVector &ex = Teuchos::dyn_cast<YUEpetraVector>(const_cast <Vector&>(x));
71  y_epetra_vec_->Update( alpha, *ex.y_epetra_vec_, beta );
72  if (u_epetra_vec_.get() != 0)
73  u_epetra_vec_->Update( alpha, *ex.u_epetra_vec_, beta );
74 }
75 
76 void YUEpetraVector::Scale( const double &alpha )
77 {
78  y_epetra_vec_->Scale( alpha );
79  if (u_epetra_vec_.get() != 0)
80  u_epetra_vec_->Scale( alpha );
81 }
82 
83 void YUEpetraVector::Set( const double &alpha )
84 {
85  y_epetra_vec_->PutScalar( alpha );
86  if (u_epetra_vec_.get() != 0)
87  u_epetra_vec_->PutScalar( alpha );
88 }
89 
90 void YUEpetraVector::Set( const double &alpha, const Vector &x )
91 {
92  YUEpetraVector &ex = Teuchos::dyn_cast<YUEpetraVector>(const_cast <Vector&>(x));
93  y_epetra_vec_->Scale( alpha, *ex.y_epetra_vec_ );
94  if (u_epetra_vec_.get() != 0)
95  u_epetra_vec_->Scale( alpha, *ex.u_epetra_vec_ );
96 }
97 
98 Teuchos::RCP<Vector> YUEpetraVector::createVector() const
99 {
100  Teuchos::RCP<Epetra_MultiVector> yptr =
101  Teuchos::rcp(new Epetra_MultiVector(y_epetra_vec_->Map(),1,false));
102  Teuchos::RCP<Epetra_MultiVector> uptr = Teuchos::null;
103  if (u_epetra_vec_.get() != 0)
104  uptr = Teuchos::rcp(new Epetra_MultiVector(u_epetra_vec_->Map(),1,false));
105 
106  return Teuchos::rcp( new YUEpetraVector( yptr, uptr ));
107 }
108 
109 Teuchos::RCP<const Epetra_MultiVector> YUEpetraVector::getYVector() const
110 {
111  return y_epetra_vec_;
112 }
113 
114 Teuchos::RCP<const Epetra_MultiVector> YUEpetraVector::getUVector() const
115 {
116  return u_epetra_vec_;
117 }
118 
119 } // namespace GenSQP
Teuchos::RCP< const Epetra_MultiVector > getYVector() const
Returns a reference counted pointer to the private y_epetra_vec data container (&quot;state variables&quot;)...
Teuchos::RCP< Vector > createVector() const
Clone to make a new (uninitialized) vector.
double innerProd(const Vector &x) const
Returns inner(*this,x).
Provides the interface to generic abstract vector libraries.
void linComb(const double &alpha, const Vector &x, const double &beta)
y = alpha*x + beta*y where y == *this.
void Scale(const double &alpha)
y = alpha*y where y == *this.
void Set(const double &alpha)
y = alpha where y == *this.
YUEpetraVector(const Teuchos::RCP< Epetra_MultiVector > &y_epetra_vec, const Teuchos::RCP< Epetra_MultiVector > &u_epetra_vec)
The GenSQP::Vector / (y,u) Epetra_MultiVector adapter class.
Teuchos::RCP< const Epetra_MultiVector > getUVector() const
Returns a reference counted pointer to the private u_epetra_vec data container (&quot;control variables&quot;)...