FEApp_SGMeanPrecOp.hpp

Go to the documentation of this file.

// $Id: FEApp_SGMeanPrecOp.hpp,v 1.1 2008/08/01 22:57:12 etphipp Exp $ 
// $Source: /space/CVS/Trilinos/packages/sacado/example/FEApp/FEApp_SGMeanPrecOp.hpp,v $ 
// @HEADER
// ***********************************************************************
// 
//                           Sacado Package
//                 Copyright (2006) Sandia Corporation
// 
// Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
// the U.S. Government retains certain rights in this software.
// 
// 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., 59 Temple Place, Suite 330, Boston, MA 02111-1307
// USA
// Questions? Contact David M. Gay (dmgay@sandia.gov) or Eric T. Phipps
// (etphipp@sandia.gov).
// 
// ***********************************************************************
// @HEADER

#ifndef FEAPP_SGMEANPRECOP_HPP
#define FEAPP_SGMEANPRECOP_HPP

#include "FEApp_TemplateTypes.hpp"

#if SGFAD_ACTIVE

#include "Teuchos_RCP.hpp"

#include "Epetra_Operator.h"
#include "Epetra_Map.h"
#include "Epetra_CrsMatrix.h"
#include "Epetra_MultiVector.h"
#include "EpetraExt_BlockMultiVector.h"
#include "Ifpack.h"

namespace FEApp {
    
  class SGMeanPrecOp : public Epetra_Operator {
      
  public:

    SGMeanPrecOp(const Teuchos::RCP<const Epetra_Map>& base_map,
     const Teuchos::RCP<const Epetra_Map>& sg_map,
     unsigned int num_blocks,
     const Teuchos::RCP<Epetra_CrsMatrix>& mean_jac,
     const Teuchos::RCP<Teuchos::ParameterList>& precParams);
    
    virtual ~SGMeanPrecOp();

    int reset();

    Teuchos::RCP<Epetra_CrsMatrix> getMeanJacobian();
    
    virtual int SetUseTranspose(bool UseTranspose);
    
    virtual int Apply(const Epetra_MultiVector& Input, 
          Epetra_MultiVector& Result) const;

    virtual int ApplyInverse(const Epetra_MultiVector& X, 
           Epetra_MultiVector& Y) const;
    
    virtual double NormInf() const;
    
    virtual const char* Label () const;
  
    virtual bool UseTranspose() const;
    
    virtual bool HasNormInf() const;

    virtual const Epetra_Comm & Comm() const;

    virtual const Epetra_Map& OperatorDomainMap () const;

    virtual const Epetra_Map& OperatorRangeMap () const;

  private:
    
    SGMeanPrecOp(const SGMeanPrecOp&);
    
    SGMeanPrecOp& operator=(const SGMeanPrecOp&);
    
  protected:
    
    string label;
    
    Teuchos::RCP<const Epetra_Map> base_map;

    Teuchos::RCP<const Epetra_Map> sg_map;

    Teuchos::RCP<Epetra_CrsMatrix> mean_jac;

    Teuchos::RCP<Teuchos::ParameterList> precParams;

    bool useTranspose;

    unsigned int num_blocks;

    mutable Teuchos::RCP<EpetraExt::BlockMultiVector> sg_input;

    mutable Teuchos::RCP<EpetraExt::BlockMultiVector> sg_result;

    mutable std::vector< Teuchos::RCP<Epetra_MultiVector> > input_block;

    mutable std::vector< Teuchos::RCP<Epetra_MultiVector> > result_block;

    Teuchos::RCP<Ifpack_Preconditioner> ifpackPrec;

  }; // class SGMeanPrecOp
  
} // namespace FEApp

#endif // SGFAD_ACTIVE

#endif // FEAPP_SGMEANPRECOP_HPP

---