src/eigsolv/QuadraticEigsolvOperators.h

Go to the documentation of this file.

#ifndef QUADRATICEIGSOLVOPERATORS_H
#define QUADRATICEIGSOLVOPERATORS_H

#include <iostream>
#include <fstream>
#include <sstream>
#include <cstdlib>
#include <cassert>
#include <unistd.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_Export.h"
#include "Epetra_LinearProblem.h"
#include "Epetra_Map.h"
#include <Teuchos_ParameterList.hpp>

#include "AbstractEigsolvOperators.h"
#include "HierarchicalBasisPrec.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 QuadraticEigsolvOperators : public AbstractEigsolvOperators {

public:

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

    ~QuadraticEigsolvOperators();

    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:

    QuadraticEigsolvOperators(const QuadraticEigsolvOperators& leo);
    QuadraticEigsolvOperators& operator=(const QuadraticEigsolvOperators &ref);   

    int buildPreconditioner(double sigma);

    int buildHierarchicalPrec(double sigma);

    int buildHierarchicalHPrec();

    Epetra_CrsMatrix* getH11();

    void matrixAssembly(FemaxMesh& mesh_mixed);

    const Epetra_Comm& Comm;
    Epetra_FECrsMatrix* A;
    Epetra_FECrsMatrix* M;
    Epetra_CrsMatrix* K;
    Epetra_CrsMatrix* K11;
    Epetra_CrsMatrix* K12;
    Epetra_CrsMatrix* K22;
    Epetra_CrsMatrix* H;
    Epetra_CrsMatrix* H11;
    Epetra_CrsMatrix* H12;
    Epetra_CrsMatrix* H22;
    Epetra_FECrsMatrix* Y;
    Epetra_CrsMatrix* Y11;
    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;
    
    int non;
    int noe;
};

#endif

---