00001 // @HEADER 00002 // *********************************************************************** 00003 // 00004 // Anasazi: Block Eigensolvers Package 00005 // Copyright (2004) Sandia Corporation 00006 // 00007 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive 00008 // license for use of this work by or on behalf of the U.S. Government. 00009 // 00010 // This library is free software; you can redistribute it and/or modify 00011 // it under the terms of the GNU Lesser General Public License as 00012 // published by the Free Software Foundation; either version 2.1 of the 00013 // License, or (at your option) any later version. 00014 // 00015 // This library is distributed in the hope that it will be useful, but 00016 // WITHOUT ANY WARRANTY; without even the implied warranty of 00017 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00018 // Lesser General Public License for more details. 00019 // 00020 // You should have received a copy of the GNU Lesser General Public 00021 // License along with this library; if not, write to the Free Software 00022 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 00023 // USA 00024 // Questions? Contact Michael A. Heroux (maherou@sandia.gov) 00025 // 00026 // *********************************************************************** 00027 // @HEADER 00028 00029 #ifndef ANASAZI_EIGENPROBLEM_H 00030 #define ANASAZI_EIGENPROBLEM_H 00031 00037 #include "AnasaziConfigDefs.hpp" 00038 #include "AnasaziTypes.hpp" 00039 #include "Teuchos_SerialDenseMatrix.hpp" 00040 #include "Teuchos_RCP.hpp" 00041 00042 00048 namespace Anasazi { 00049 00050 template<class ScalarType, class MV, class OP> 00051 class Eigenproblem { 00052 00053 public: 00054 00056 00057 00059 Eigenproblem() {}; 00060 00062 virtual ~Eigenproblem() {}; 00064 00066 00067 00075 virtual void setOperator( const Teuchos::RCP<const OP> &Op ) = 0; 00076 00078 virtual void setA( const Teuchos::RCP<const OP> &A ) = 0; 00079 00081 virtual void setM( const Teuchos::RCP<const OP> &M ) = 0; 00082 00084 virtual void setPrec( const Teuchos::RCP<const OP> &Prec ) = 0; 00085 00090 virtual void setInitVec( const Teuchos::RCP<MV> &InitVec ) = 0; 00091 00098 virtual void setAuxVecs( const Teuchos::RCP<const MV> &AuxVecs ) = 0; 00099 00101 virtual void setNEV( int nev ) = 0; 00102 00108 virtual void setHermitian( bool isSym ) = 0; 00109 00123 virtual bool setProblem() = 0; 00124 00131 virtual void setSolution(const Eigensolution<ScalarType,MV> &sol) = 0; 00132 00134 00136 00137 00139 virtual Teuchos::RCP<const OP> getOperator() const = 0; 00140 00142 virtual Teuchos::RCP<const OP> getA() const = 0; 00143 00145 virtual Teuchos::RCP<const OP> getM() const = 0; 00146 00148 virtual Teuchos::RCP<const OP> getPrec() const = 0; 00149 00151 virtual Teuchos::RCP<const MV> getInitVec() const = 0; 00152 00154 virtual Teuchos::RCP<const MV> getAuxVecs() const = 0; 00155 00157 virtual int getNEV() const = 0; 00158 00160 virtual bool isHermitian() const = 0; 00161 00163 virtual bool isProblemSet() const = 0; 00164 00170 virtual const Eigensolution<ScalarType,MV> & getSolution() const = 0; 00171 00173 }; 00174 00175 } // end Anasazi namespace 00176 #endif 00177 00178 // end AnasaziEigenproblem.hpp