FEI Version of the Day
fei_Solver_Belos.hpp
00001 /*
00002 // @HEADER
00003 // ************************************************************************
00004 //             FEI: Finite Element Interface to Linear Solvers
00005 //                  Copyright (2005) Sandia Corporation.
00006 //
00007 // Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation, the
00008 // U.S. Government retains certain rights in this software.
00009 //
00010 // Redistribution and use in source and binary forms, with or without
00011 // modification, are permitted provided that the following conditions are
00012 // met:
00013 //
00014 // 1. Redistributions of source code must retain the above copyright
00015 // notice, this list of conditions and the following disclaimer.
00016 //
00017 // 2. Redistributions in binary form must reproduce the above copyright
00018 // notice, this list of conditions and the following disclaimer in the
00019 // documentation and/or other materials provided with the distribution.
00020 //
00021 // 3. Neither the name of the Corporation nor the names of the
00022 // contributors may be used to endorse or promote products derived from
00023 // this software without specific prior written permission.
00024 //
00025 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
00026 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
00027 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
00028 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
00029 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
00030 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
00031 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
00032 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
00033 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
00034 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
00035 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
00036 //
00037 // Questions? Contact Alan Williams (william@sandia.gov) 
00038 //
00039 // ************************************************************************
00040 // @HEADER
00041 */
00042 
00043 #ifndef _fei_Solver_Belos_hpp_
00044 #define _fei_Solver_Belos_hpp_
00045 
00046 
00047 #include <fei_trilinos_macros.hpp>
00048 
00049 #ifdef HAVE_FEI_BELOS
00050 
00051 #include <fei_macros.hpp>
00052 #include <fei_Solver.hpp>
00053 #include <fei_Logger.hpp>
00054 
00055 #ifdef HAVE_FEI_TEUCHOS
00056 namespace Teuchos {
00057   class ParameterList;
00058 }
00059 #endif
00060 
00061 class Epetra_CrsMatrix;
00062 namespace Belos {
00063 template<class Scalar,class MV,class OP> class SolverManager;
00064 }
00065 
00066 #ifdef HAVE_FEI_ML
00067 #include <ml_include.h>
00068 #include <ml_epetra_preconditioner.h>
00069 #endif
00070 
00073 class Solver_Belos : public fei::Solver, private fei::Logger {
00074  public:
00078   Solver_Belos();
00079 
00082   virtual ~Solver_Belos();
00083 
00097   int solve(fei::LinearSystem* linearSystem,
00098       fei::Matrix* preconditioningMatrix,
00099       int numParams,
00100       const char* const* solverParams,
00101       int& iterationsTaken,
00102       int& status);
00103 
00120   int solve(fei::LinearSystem* linearSystem,
00121       fei::Matrix* preconditioningMatrix,
00122       const fei::ParameterSet& parameterSet,
00123       int& iterationsTaken,
00124       int& status);
00125 
00126   Teuchos::ParameterList& get_ParameterList();
00127 
00128   void setMaxIters(int maxits) {maxIters_ = maxits;}
00129   void setTolerance(double tol) {tolerance_ = tol;}
00130   void setUseTranspose(bool useTrans) {useTranspose_ = useTrans;}
00131 
00132   void setUseML(bool useml);
00133 
00134  private:
00135   int setup_ml_operator(AztecOO& azoo, Epetra_CrsMatrix* A);
00136 
00137  private:
00138   Teuchos::RCP<Belos::SolverManager<double,Epetra_MultiVector,Epetra_Operator> > belos_solver_manager_;
00139   double tolerance_;
00140   int maxIters_;
00141   bool useTranspose_;
00142   Teuchos::RCP<Teuchos::ParameterList> paramlist_;
00143 
00144   bool useML_;
00145 #ifdef HAVE_FEI_ML
00146   ML_Epetra::MultiLevelPreconditioner* ml_prec_;
00147   bool ml_defaults_set_;
00148   int *ml_aztec_options_;
00149   double *ml_aztec_params_;
00150 #endif
00151 
00152   std::string name_;
00153   std::string dbgprefix_;
00154 }; //class Solver_Belos
00155 
00156 #endif // HAVE_FEI_BELOS
00157 
00158 #endif
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Friends