src/eigsolv/LinSolvers.h

Go to the documentation of this file.

//**************************************************************************
//
//                                 NOTICE
//
// This software is a result of the research described in the report
//
// " A comparison of algorithms for modal analysis in the absence 
//   of a sparse direct method", P. Arbenz, R. Lehoucq, and U. Hetmaniuk,
//  Sandia National Laboratories, Technical report SAND2003-1028J.
//
// It is based on the Epetra, AztecOO, and ML packages defined in the Trilinos
// framework ( http://software.sandia.gov/trilinos/ ).
//
// The distribution of this software follows also the rules defined in Trilinos.
// This notice shall be marked on any reproduction of this software, in whole or
// in part.
//
// 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 program 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.
//
// Code Authors: U. Hetmaniuk (ulhetma@sandia.gov), R. Lehoucq (rblehou@sandia.gov)
//
//**************************************************************************
/**************************************************************************
 *                                                                         *
 *               Swiss Federal Institute of Technology (ETH)               *
 *                       CH-8092 Zuerich, Switzerland                      *
 *                                                                         *
 *                       (C) 1999 All Rights Reserved                      *
 *                                                                         *
 *                                NOTICE                                   *
 *                                                                         *
 *  Permission to use, copy, modify, and distribute this software and      *
 *  its documentation for any purpose and without fee is hereby granted    *
 *  provided that the above copyright notice appear in all copies and      *
 *  that both the copyright notice and this permission notice appear in    *
 *  supporting documentation.                                              *
 *                                                                         *
 *  Neither the Swiss Federal Institute of Technology nor the author make  *
 *  any representations about the suitability of this software for any     *
 *  purpose.  This software is provided ``as is'' without express or       *
 *  implied warranty.                                                      *
 *                                                                         *
 ***************************************************************************/

#ifndef LIN_SOLVERS_H
#define LIN_SOLVERS_H

#include "Epetra_ConfigDefs.h"

#include "Epetra_Operator.h"
#include "Epetra_Vector.h"

class LinSolvers {

private:

    const Epetra_Operator *A;
    const Epetra_Operator *Prec;

    double *work_s;

    int solverType;

    static const int LINSOLV_CGS;       // not implemented
    static const int LINSOLV_MINRES;    // not implemented
    static const int LINSOLV_SYMMLQ;    // not implemented
    static const int LINSOLV_QMR;

    // Don't define these functions
    LinSolvers(const LinSolvers &L);
    LinSolvers & operator=(const LinSolvers &L);

protected:

    int QMRS(const Epetra_Vector &b, Epetra_Vector &x, double _tol, int _itMax);

public:

    LinSolvers(const Epetra_Operator *AA, const Epetra_Operator *PP, double *wTmp,
               int _type = LINSOLV_QMR);

    ~LinSolvers() { }

    int solve(const Epetra_Vector &b, Epetra_Vector &x, double _tol, int _itMax);

};

#endif

---