SpoolesOO.h

Go to the documentation of this file.

// @HEADER
// ***********************************************************************
// 
//                Amesos: Direct Sparse Solver Package
//                 Copyright (2004) Sandia Corporation
// 
// Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
// license for use of this work by or on behalf of the U.S. Government.
// 
// This library is free software; you can redistribute it and/or modify
// it under the terms of the GNU Lesser General Public License as
// published by the Free Software Foundation; either version 2.1 of the
// License, or (at your option) any later version.
//  
// This library is distributed in the hope that it will be useful, but
// WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
// Lesser General Public License for more details.
//  
// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301
// USA
// Questions? Contact Michael A. Heroux (maherou@sandia.gov) 
// 
// ***********************************************************************
// @HEADER
 
#ifndef _SPOOLESOO_H_
#define _SPOOLESOO_H_
 
class Epetra_Comm;
class Epetra_BlockMap;
class Epetra_MultiVector;
class Epetra_RowMatrix;
//  #include "Epetra_LinearProblem.h"
#include "Epetra_Operator.h"
#include "Epetra_Object.h"
 


 
class SpoolesOO {
    
  public:
  SpoolesOO(Epetra_RowMatrix * A, Epetra_MultiVector * X, Epetra_MultiVector * B);
 
  SpoolesOO();
 
  virtual ~SpoolesOO(void);
  
  int SetUserMatrix(Epetra_RowMatrix * UserMatrix);
 
  int SetLHS(Epetra_MultiVector * X);
 
  int SetRHS(Epetra_MultiVector * B);
 
  Epetra_RowMatrix * GetUserMatrix() const {return(UserMatrix_);};
 
  Epetra_MultiVector * GetLHS() const {return(X_);};
 
  Epetra_MultiVector * GetRHS() const {return(B_);};
 
  bool GetTrans( ) const { return Transpose_ ;} ;
 
  void SetTrans( bool trans ) { Transpose_ = trans ;} ; 
 
  int SetSpoolesDefaults();
 
  int Solve() ;
 
 protected:
 
  Epetra_Operator * UserOperator_;
  Epetra_RowMatrix * UserMatrix_;
  Epetra_Operator * PrecOperator_;
  Epetra_RowMatrix * PrecMatrix_;
  Epetra_MultiVector * X_;
  Epetra_MultiVector * B_;
 
  
  bool Transpose_ ;
 
  int x_LDA_;
  double *x_;
  int b_LDA_;
  double *b_;
  bool inConstructor_;
};
 
 
#endif /* _SPOOLESOO_H_ */