EpetraExt Development
Loading...
Searching...
No Matches
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
46namespace GenSQP {
47
48YUEpetraVector::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
55double 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
68void 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
76void 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
83void 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
90void 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
98Teuchos::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
109Teuchos::RCP<const Epetra_MultiVector> YUEpetraVector::getYVector() const
110{
111 return y_epetra_vec_;
112}
113
114Teuchos::RCP<const Epetra_MultiVector> YUEpetraVector::getUVector() const
115{
116 return u_epetra_vec_;
117}
118
119} // namespace GenSQP
Provides the interface to generic abstract vector libraries.
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 ("control variables").
double innerProd(const Vector &x) const
Returns inner(*this,x).
Teuchos::RCP< Vector > createVector() const
Clone to make a new (uninitialized) vector.
void Scale(const double &alpha)
y = alpha*y where y == *this.
void linComb(const double &alpha, const Vector &x, const double &beta)
y = alpha*x + beta*y where y == *this.
Teuchos::RCP< const Epetra_MultiVector > getYVector() const
Returns a reference counted pointer to the private y_epetra_vec data container ("state variables").
YUEpetraVector(const Teuchos::RCP< Epetra_MultiVector > &y_epetra_vec, const Teuchos::RCP< Epetra_MultiVector > &u_epetra_vec)
void Set(const double &alpha)
y = alpha where y == *this.