FEApp_SGMatrixFreeOp.hpp

Go to the documentation of this file.

// $Id: FEApp_SGMatrixFreeOp.hpp,v 1.1 2008/08/01 22:57:12 etphipp Exp $ 
// $Source: /space/CVS/Trilinos/packages/sacado/example/FEApp/FEApp_SGMatrixFreeOp.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_SGMATRIXFREEOP_HPP
#define FEAPP_SGMATRIXFREEOP_HPP

#include "FEApp_TemplateTypes.hpp"

#if SGFAD_ACTIVE

#include "Teuchos_RCP.hpp"

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

namespace FEApp {
    
  class SGMatrixFreeOp : public Epetra_Operator {
      
  public:

    typedef SGType::expansion_type::tp_type tp_type;

    SGMatrixFreeOp(const Teuchos::RCP<const Epetra_Map>& base_map,
       const Teuchos::RCP<const Epetra_Map>& sg_map,
       const Teuchos::RCP<const tp_type >& Cijk,
       const Teuchos::RCP<std::vector< Teuchos::RCP<Epetra_CrsMatrix> > >& jacs);
    
    virtual ~SGMatrixFreeOp();

    virtual std::vector< Teuchos::RCP<Epetra_CrsMatrix> >&
    getJacobianBlocks();
    
    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:
    
    SGMatrixFreeOp(const SGMatrixFreeOp&);
    
    SGMatrixFreeOp& operator=(const SGMatrixFreeOp&);
    
  protected:
    
    string label;
    
    Teuchos::RCP<const Epetra_Map> base_map;

    Teuchos::RCP<const Epetra_Map> sg_map;

    Teuchos::RCP<const tp_type > Cijk;

    Teuchos::RCP<std::vector< Teuchos::RCP<Epetra_CrsMatrix> > > jacs;

    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;

    mutable Teuchos::RCP<Epetra_MultiVector> tmp;

  }; // class SGMatrixFreeOp
  
} // namespace FEApp

#endif // SGFAD_ACTIVE

#endif // FEAPP_SGMATRIXFREEOP_HPP

---