00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042 #ifndef LOCA_EPETRA_LOWRANKUPDATEROWMATRIX_H
00043 #define LOCA_EPETRA_LOWRANKUPDATEROWMATRIX_H
00044
00045 #include "Teuchos_RCP.hpp"
00046
00047 #include "LOCA_Epetra_LowRankUpdateOp.H"
00048 #include "Epetra_RowMatrix.h"
00049
00050 namespace LOCA {
00051
00052 namespace Epetra {
00053
00074 class LowRankUpdateRowMatrix :
00075 public LOCA::Epetra::LowRankUpdateOp,
00076 public virtual Epetra_RowMatrix {
00077
00078 public:
00079
00081
00091 LowRankUpdateRowMatrix(
00092 const Teuchos::RCP<LOCA::GlobalData>& global_data,
00093 const Teuchos::RCP<Epetra_RowMatrix>& jacRowMatrix,
00094 const Teuchos::RCP<Epetra_MultiVector>& U_multiVec,
00095 const Teuchos::RCP<Epetra_MultiVector>& V_multiVec,
00096 bool setup_for_solve,
00097 bool include_UV_terms);
00098
00100 virtual ~LowRankUpdateRowMatrix();
00101
00103 virtual const Epetra_BlockMap & Map() const;
00104
00106
00108 virtual int NumMyRowEntries(int MyRow, int & NumEntries) const;
00109
00110
00112 virtual int MaxNumEntries() const;
00113
00115 virtual int ExtractMyRowCopy(int MyRow, int Length, int & NumEntries,
00116 double *Values, int * Indices) const;
00117
00119 virtual int ExtractDiagonalCopy(Epetra_Vector & Diagonal) const;
00120
00122
00124
00129 virtual int Multiply(bool TransA, const Epetra_MultiVector& X,
00130 Epetra_MultiVector& Y) const;
00131
00136 virtual int Solve(bool Upper, bool Trans, bool UnitDiagonal,
00137 const Epetra_MultiVector& X,
00138 Epetra_MultiVector& Y) const;
00139
00144 virtual int InvRowSums(Epetra_Vector& x) const;
00145
00147 virtual int LeftScale(const Epetra_Vector& x);
00148
00153 virtual int InvColSums(Epetra_Vector& x) const;
00154
00156 virtual int RightScale(const Epetra_Vector& x);
00157
00159
00161
00166 virtual bool Filled() const;
00167
00169 virtual double NormInf() const;
00170
00172 virtual double NormOne() const;
00173
00175 virtual int NumGlobalNonzeros() const;
00176
00178 virtual int NumGlobalRows() const;
00179
00181 virtual int NumGlobalCols() const;
00182
00187 virtual int NumGlobalDiagonals() const;
00188
00193 virtual int NumMyNonzeros() const;
00194
00196 virtual int NumMyRows() const;
00197
00199 virtual int NumMyCols() const;
00200
00205 virtual int NumMyDiagonals() const;
00206
00211 virtual bool LowerTriangular() const;
00212
00217 virtual bool UpperTriangular() const;
00218
00223 virtual const Epetra_Map & RowMatrixRowMap() const;
00224
00229 virtual const Epetra_Map & RowMatrixColMap() const;
00230
00235 virtual const Epetra_Import * RowMatrixImporter() const;
00236
00238
00239 protected:
00240
00242 double computeUV(int MyRow, int MyCol) const;
00243
00244 private:
00245
00247 LowRankUpdateRowMatrix(const LowRankUpdateRowMatrix&);
00248
00250 LowRankUpdateRowMatrix& operator=(const LowRankUpdateRowMatrix&);
00251
00252 protected:
00253
00255 Teuchos::RCP<Epetra_RowMatrix> J_rowMatrix;
00256
00258 Teuchos::RCP<Epetra_MultiVector> nonconst_U;
00259
00261 Teuchos::RCP<Epetra_MultiVector> nonconst_V;
00262
00264 bool includeUV;
00265
00267 int m;
00268
00270 const Epetra_BlockMap& U_map;
00271
00273 const Epetra_BlockMap& V_map;
00274
00276 const Epetra_BlockMap& row_map;
00277
00278 };
00279 }
00280 }
00281
00282 #endif