EpetraExt Package Browser (Single Doxygen Collection) Development
Loading...
Searching...
No Matches
EpetraExt_BTF_LinearProblem.h
Go to the documentation of this file.
1//@HEADER
2// ***********************************************************************
3//
4// EpetraExt: Epetra Extended - Linear Algebra Services Package
5// Copyright (2011) Sandia Corporation
6//
7// Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
8// the U.S. Government retains certain rights in this software.
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 Michael A. Heroux (maherou@sandia.gov)
38//
39// ***********************************************************************
40//@HEADER
41
42#ifndef EpetraExt_LINEARPROBLEM_BTF_H
43#define EpetraExt_LINEARPROBLEM_BTF_H
44
45#include <EpetraExt_Transform.h>
46
47#include <vector>
48#include <map>
49#include <set>
50
54class Epetra_Vector;
56class Epetra_CrsGraph;
57class Epetra_Map;
58class Epetra_BlockMap;
60
61namespace EpetraExt {
62
63class LinearProblem_BTF : public SameTypeTransform<Epetra_LinearProblem> {
64
65 public:
66
68
69 LinearProblem_BTF( double thres = 0.0,
70 int verbose = 0 )
71 : NewMap_(0),
72 NewMatrix_(0),
73 NewGraph_(0),
74 NewLHS_(0),
75 NewRHS_(0),
76 NewProblem_(0),
77 OrigGraph_(0),
78 OrigMatrix_(0),
79 OrigRHS_(0),
80 OrigLHS_(0),
81 OrigProblem_(0),
82 threshold_(thres),
83 verbose_(verbose),
84 changedLP_(false)
85 {}
86
88
89 bool fwd();
90 bool rvs();
91
92 bool changedLP() { return changedLP_; }
93
94 private:
95
96 void deleteNewObjs_();
97
99
101
104
107
115
116 std::vector<int> OldGlobalElements_;
117
118 std::vector< std::set<int> > ZeroElements_;
119
120 std::vector< std::vector<Epetra_SerialDenseMatrix*> > Blocks_;
121 std::vector<int> BlockDim_;
122 std::vector<int> BlockCnt_;
123 std::map<int,int> BlockRowMap_;
124 std::map<int,int> SubBlockRowMap_;
125 std::map<int,int> BlockColMap_;
126 std::map<int,int> SubBlockColMap_;
127
128 std::vector< std::vector<int> > NewBlockRows_;
129
130 const double threshold_;
131 const int verbose_;
132
134};
135
136} //namespace EpetraExt
137
138#endif //EpetraExt_LINEARPROBLEM_BTF_H
std::vector< std::vector< Epetra_SerialDenseMatrix * > > Blocks_
NewTypeRef operator()(OriginalTypeRef orig)
Analysis of transform operation on original object and construction of new object.
std::vector< std::vector< int > > NewBlockRows_
bool rvs()
Reverse transfer of data from new object created in the operator() method call to the orig object inp...
LinearProblem_BTF(double thres=0.0, int verbose=0)
bool fwd()
Forward transfer of data from orig object input in the operator() method call to the new object creat...
std::vector< std::set< int > > ZeroElements_
EpetraExt::BlockCrsMatrix: A class for constructing a distributed block matrix.