38#ifndef IFPACK_SUBDOMAINFILTER_H
39#define IFPACK_SUBDOMAINFILTER_H
43#ifdef HAVE_IFPACK_PARALLEL_SUBDOMAIN_SOLVERS
53#include "Teuchos_RCP.hpp"
83 Ifpack_SubdomainFilter(
const Teuchos::RCP<const Epetra_RowMatrix>& Matrix,
int subdomainId);
88 ~Ifpack_SubdomainFilter();
104 NumEntries = NumEntries_[MyRow];
109 virtual int MaxNumEntries()
const
111 return(MaxNumEntries_);
129 virtual inline int ExtractMyRowCopy(
int MyRow,
int Length,
int & NumEntries,
double *Values,
int * Indices)
const;
138 virtual int ExtractDiagonalCopy(
Epetra_Vector & Diagonal)
const;
156 if (TransA ==
true) {
165 virtual int Solve(
bool Upper,
bool Trans,
bool UnitDiagonal,
const Epetra_MultiVector& X,
206 virtual bool Filled()
const
215 virtual double NormInf()
const
224 virtual double NormOne()
const
229#ifndef EPETRA_NO_32BIT_GLOBAL_INDICES
231 virtual int NumGlobalNonzeros()
const
233 return(NumGlobalNonzeros_);
237 virtual int NumGlobalRows()
const
240 return(NumGlobalRows_);
244 virtual int NumGlobalCols()
const
246 return(NumGlobalCols_);
250 virtual int NumGlobalDiagonals()
const
252 return(NumGlobalRows_);
257 virtual long long NumGlobalNonzeros64()
const
259 return(NumGlobalNonzeros_);
263 virtual long long NumGlobalRows64()
const
265 return(NumGlobalRows_);
269 virtual long long NumGlobalCols64()
const
271 return(NumGlobalRows_);
275 virtual long long NumGlobalDiagonals64()
const
277 return(NumGlobalRows_);
281 virtual int NumMyNonzeros()
const
283 return(NumMyNonzeros_);
287 virtual int NumMyRows()
const
293 virtual int NumMyCols()
const
299 virtual int NumMyDiagonals()
const
305 virtual bool LowerTriangular()
const
307 return(Matrix_->LowerTriangular());
311 virtual bool UpperTriangular()
const
313 return(Matrix_->UpperTriangular());
317 virtual const Epetra_Map & RowMatrixRowMap()
const
323 virtual const Epetra_Map & RowMatrixColMap()
const
335 virtual const Epetra_Import* Importer()
const {
return(&*Importer_);}
337 virtual const Epetra_Export* Exporter()
const {
return(&*Exporter_);}
342 int SetOwnership(
bool ownership)
348 int SetUseTranspose(
bool UseTranspose_in)
350 UseTranspose_ = UseTranspose_in;
355 bool UseTranspose()
const
357 return(UseTranspose_);
361 bool HasNormInf()
const
387const char* Label()
const{
392 void UpdateImportVector(
int NumVectors)
const;
393 void UpdateExportVector(
int NumVectors)
const;
396 Teuchos::RCP<const Epetra_RowMatrix> Matrix_;
400 Teuchos::RCP<Epetra_MpiComm> SubComm_;
401 MPI_Comm subdomainMPIComm_;
404 Teuchos::RCP<Epetra_SerialComm> SubComm_;
407 Teuchos::RCP<Epetra_Map> Map_;
409 Teuchos::RCP<Epetra_Map> colMap_;
421 int NumGlobalNonzeros_;
427 std::vector<int> NumEntries_;
429 mutable std::vector<int> Indices_;
431 mutable std::vector<double> Values_;
436 Teuchos::RCP<Epetra_Vector> Diagonal_;
457 Teuchos::RCP<Epetra_Import> Importer_;
458 Teuchos::RCP<Epetra_Export> Exporter_;
#define IFPACK_RETURN(ifpack_err)
#define IFPACK_CHK_ERR(ifpack_err)
int NumMyRowEntries(void *A, int Row, int *numEntries)
Ifpack_OverlappingRowMatrix: matrix with ghost rows, based on Epetra_RowMatrix.