DALinearEqn.H

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
 
    DAFoam  : Discrete Adjoint with OpenFOAM
    Version : v4
 
    Description:
        Solve the linear euqation using Petsc
 
\*---------------------------------------------------------------------------*/
 
#ifndef DALinearEqn_H
#define DALinearEqn_H
 
#include "fvOptions.H"
#include "surfaceFields.H"
#include "DAOption.H"
#include "DAUtility.H"
#include "DAStateInfo.H"
#include "DAModel.H"
#include "DAIndex.H"
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
namespace Foam
{
 
/*---------------------------------------------------------------------------*\
                       Class DALinearEqn Declaration
\*---------------------------------------------------------------------------*/
 
class DALinearEqn
{
 
private:
    DALinearEqn(const DALinearEqn&);
 
    void operator=(const DALinearEqn&);
 
protected:
    const fvMesh& mesh_;
 
    const DAOption& daOption_;
 
public:
    DALinearEqn(
        const fvMesh& mesh,
        const DAOption& daOption);
 
    virtual ~DALinearEqn()
    {
    }
 
    // Members
 
    void createMLRKSP(
        const Mat jacMat,
        const Mat jacPCMat,
        KSP ksp);
 
    label solveLinearEqn(
        const KSP ksp,
        const Vec rhsVec,
        Vec solVec);
 
    static PetscErrorCode myKSPMonitor(
        KSP,
        PetscInt,
        PetscReal,
        void*);
 
    double getRunTime();
 
    label getPrintInterval();
};
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
} // End namespace Foam
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#endif
 
// ************************************************************************* //