Thyra Package Browser (Single Doxygen Collection)
Version of the Day
Loading...
Searching...
No Matches
adapters
epetra
src
Thyra_EpetraOperatorViewExtractorStd.cpp
Go to the documentation of this file.
1
// @HEADER
2
// ***********************************************************************
3
//
4
// Thyra: Interfaces and Support for Abstract Numerical Algorithms
5
// Copyright (2004) Sandia Corporation
6
//
7
// Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
8
// license for use of this work by or on behalf of the U.S. Government.
9
//
10
// Redistribution and use in source and binary forms, with or without
11
// modification, are permitted provided that the following conditions are
12
// met:
13
//
14
// 1. Redistributions of source code must retain the above copyright
15
// notice, this list of conditions and the following disclaimer.
16
//
17
// 2. Redistributions in binary form must reproduce the above copyright
18
// notice, this list of conditions and the following disclaimer in the
19
// documentation and/or other materials provided with the distribution.
20
//
21
// 3. Neither the name of the Corporation nor the names of the
22
// contributors may be used to endorse or promote products derived from
23
// this software without specific prior written permission.
24
//
25
// THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
26
// EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
27
// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
28
// PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
29
// CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
30
// EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
31
// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
32
// PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
33
// LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
34
// NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
35
// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
36
//
37
// Questions? Contact Roscoe A. Bartlett (bartlettra@ornl.gov)
38
//
39
// ***********************************************************************
40
// @HEADER
41
42
#include "
Thyra_EpetraOperatorViewExtractorStd.hpp
"
43
#include "
Thyra_EpetraLinearOpBase.hpp
"
44
#include "Thyra_ScaledAdjointLinearOpBase.hpp"
45
46
class
Epetra_Operator;
47
48
49
namespace
Thyra
{
50
51
52
// Overridden from EpetraOperatorViewExtractorBase
53
54
55
bool
EpetraOperatorViewExtractorStd::isCompatible
(
const
LinearOpBase<double>
&
fwdOp
)
const
56
{
57
double
wrappedScalar
= 0.0;
58
EOpTransp
wrappedTransp
= NOTRANS;
59
const
LinearOpBase<double>
*
wrappedFwdOp
=
NULL
;
60
::Thyra::unwrap
(
fwdOp
, &
wrappedScalar
, &
wrappedTransp
, &
wrappedFwdOp
);
61
const
EpetraLinearOpBase
*
eFwdOp
=
62
dynamic_cast<
const
EpetraLinearOpBase
*
>
(
wrappedFwdOp
);
63
if
(!
eFwdOp
) {
64
return
false
;
65
}
66
return
true
;
67
}
68
69
70
void
EpetraOperatorViewExtractorStd::getNonconstEpetraOpView
(
71
const
RCP
<
LinearOpBase<double>
> &
/* fwdOp */
,
72
const
Ptr<
RCP<Epetra_Operator>
> &
/* epetraOp */
,
73
const
Ptr<EOpTransp>
&
/* epetraOpTransp */
,
74
const
Ptr<EApplyEpetraOpAs>
&
/* epetraOpApplyAs */
,
75
const
Ptr<EAdjointEpetraOp>
&
/* epetraOpAdjointSupport */
,
76
const
Ptr<double>
&
/* epetraOpScalar */
77
)
const
78
{
79
TEUCHOS_TEST_FOR_EXCEPT
(
true
);
80
// ToDo: Implement once this is needed by just copying what is below and
81
// removing the 'const' in the right places!
82
}
83
84
85
void
EpetraOperatorViewExtractorStd::getEpetraOpView
(
86
const
RCP
<
const
LinearOpBase<double>
> &
fwdOp
,
87
const
Ptr<
RCP<const Epetra_Operator>
> &
epetraOp
,
88
const
Ptr<EOpTransp>
&
epetraOpTransp
,
89
const
Ptr<EApplyEpetraOpAs>
&
epetraOpApplyAs
,
90
const
Ptr<EAdjointEpetraOp>
&
epetraOpAdjointSupport
,
91
const
Ptr<double>
&
epetraOpScalar
92
)
const
93
{
94
using
Teuchos::outArg;
95
double
wrappedFwdOpScalar
= 0.0;
96
EOpTransp
wrappedFwdOpTransp
= NOTRANS;
97
Teuchos::RCP<const LinearOpBase<double>
>
wrappedFwdOp
;
98
unwrap
(
fwdOp
,&
wrappedFwdOpScalar
, &
wrappedFwdOpTransp
, &
wrappedFwdOp
);
99
Teuchos::RCP<const EpetraLinearOpBase>
epetraFwdOp
=
100
Teuchos::rcp_dynamic_cast<const EpetraLinearOpBase>
(
wrappedFwdOp
,
true
);
101
EOpTransp
epetra_epetraOpTransp
;
102
epetraFwdOp
->getEpetraOpView(
epetraOp
, outArg(
epetra_epetraOpTransp
),
103
epetraOpApplyAs
,
epetraOpAdjointSupport
);
104
*
epetraOpTransp
=
trans_trans
(
real_trans
(
epetra_epetraOpTransp
),
wrappedFwdOpTransp
);
105
*
epetraOpScalar
=
wrappedFwdOpScalar
;
106
}
107
108
109
// ToDo: Refactor unwrap(...) to not use raw pointers!
110
111
112
}
// namespace Thyra
Thyra_EpetraLinearOpBase.hpp
Thyra_EpetraOperatorViewExtractorStd.hpp
Teuchos::RCP
Thyra::EpetraLinearOpBase
Abstract base class for all LinearOpBase objects that can return an Epetra_Operator view of themselve...
Definition
Thyra_EpetraLinearOpBase.hpp:69
Thyra::EpetraOperatorViewExtractorStd::getEpetraOpView
void getEpetraOpView(const RCP< const LinearOpBase< double > > &fwdOp, const Ptr< RCP< const Epetra_Operator > > &epetraOp, const Ptr< EOpTransp > &epetraOpTransp, const Ptr< EApplyEpetraOpAs > &epetraOpApplyAs, const Ptr< EAdjointEpetraOp > &epetraOpAdjointSupport, const Ptr< double > &epetraOpScalar) const
Definition
Thyra_EpetraOperatorViewExtractorStd.cpp:85
Thyra::EpetraOperatorViewExtractorStd::isCompatible
bool isCompatible(const LinearOpBase< double > &fwdOp) const
Definition
Thyra_EpetraOperatorViewExtractorStd.cpp:55
Thyra::EpetraOperatorViewExtractorStd::getNonconstEpetraOpView
void getNonconstEpetraOpView(const RCP< LinearOpBase< double > > &fwdOp, const Ptr< RCP< Epetra_Operator > > &epetraOp, const Ptr< EOpTransp > &epetraOpTransp, const Ptr< EApplyEpetraOpAs > &epetraOpApplyAs, const Ptr< EAdjointEpetraOp > &epetraOpAdjointSupport, const Ptr< double > &epetraOpScalar) const
Definition
Thyra_EpetraOperatorViewExtractorStd.cpp:70
TEUCHOS_TEST_FOR_EXCEPT
#define TEUCHOS_TEST_FOR_EXCEPT(throw_exception_test)
Thyra
Definition
Thyra_DiagonalEpetraLinearOpWithSolveFactory.cpp:53
Generated by
1.10.0