src/eigsolv/LinearEigsolvOperators.h

Go to the documentation of this file.

#ifndef LINEAREIGSOLVOPERATORS_H
#define LINEAREIGSOLVOPERATORS_H

#include "AbstractEigsolvOperators.h"

#ifdef HAVE_MPI
#include "mpi.h"
#include "Epetra_MpiComm.h"
#else
#include "Epetra_SerialComm.h"
#endif

#include "AztecOO.h"
#include "AztecOO_Operator.h"

#include "Epetra_ConfigDefs.h"
#include "Epetra_Operator.h"
#include "Epetra_CrsMatrix.h"
#include "Epetra_FECrsMatrix.h"
#include "Epetra_Export.h"
#include "Epetra_LinearProblem.h"
#include "Epetra_Map.h"
#include "Teuchos_ParameterList.hpp"

#include <iostream>
#include <fstream>
#include <sstream>
#include <cstdlib>
#include <cassert>
#include <unistd.h>

#include "debugmeshbuilder.h"
#include "tetmeshbuilder.h"
#include "broadcastmeshbuilder.h"
#include "netgenreader.h"
#include "femaxmesh.h"
#include "nedelecmesh.h"
#include "lagrangemesh.h"
#include "export_Epetra_CrsMatrix.h"

#define USE_PM_MATRIXMATRIX 1

using namespace std;
using namespace mesh;


class LinearEigsolvOperators : public AbstractEigsolvOperators {

public:

  LinearEigsolvOperators(FemaxMesh& mesh_mixed, Teuchos::ParameterList& params, const Epetra_Comm& MyComm);
    ~LinearEigsolvOperators();

    const Epetra_Operator* getA() const; 
    const Epetra_Operator* getM() const;
    const Epetra_Operator* getAsigmaPrec(double sigma);
    const Epetra_Operator* getHSolver();
    const Epetra_Operator* getHPrec();
    Epetra_Operator* getH();
    Epetra_Operator* getY() const;
    Epetra_Operator* getC() const;

private:
    LinearEigsolvOperators(const LinearEigsolvOperators& leo);
    LinearEigsolvOperators& operator=(const LinearEigsolvOperators &ref);   
    int buildPreconditioner(double sigma);
    void matrixAssembly(FemaxMesh& mesh_mixed);

    FemaxMesh& mesh_mixed_;
    const Epetra_Comm& Comm;
    Epetra_FECrsMatrix* A;
    Epetra_FECrsMatrix* M;
    Epetra_CrsMatrix* K;
    Epetra_CrsMatrix* H;
    Epetra_FECrsMatrix* Y;
    Epetra_CrsMatrix* C;
    Epetra_Operator* Prec;

    Epetra_Operator* h_solver_op_;
    Epetra_Operator* HPrec;
    double _sigma_K;
    Teuchos::ParameterList& params_;
    
    // The next two objects are defined for using 
    // the Neumann preconditioner from AztecOO.
    // They might disappear when ML is operational.
    Epetra_LinearProblem* Aprec;
    AztecOO* AprecSolver;

};

#endif

---