DAJacCondRdW.H

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
 
    DAFoam  : Discrete Adjoint with OpenFOAM
    Version : v3
 
    Description:
        Child class for dRdW
 
\*---------------------------------------------------------------------------*/
 
#ifndef DAJacCondRdW_H
#define DAJacCondRdW_H
 
#include "DAJacCon.H"
#include "addToRunTimeSelectionTable.H"
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
namespace Foam
{
 
/*---------------------------------------------------------------------------*\
      Class DAJacCondRdW Declaration
\*---------------------------------------------------------------------------*/
 
class DAJacCondRdW
    : public DAJacCon
{
 
protected:
 
    Vec dRdWTPreallocOn_;
    Vec dRdWTPreallocOff_;
    Vec dRdWPreallocOn_;
    Vec dRdWPreallocOff_;
 
    void initializePetscVecs();
 
    void allocateJacobianConnections(
        Vec preallocOnProc,
        Vec preallocOffProc,
        Vec preallocOnProcT,
        Vec preallocOffProcT,
        Mat connections,
        const label row);
 
    void preallocateJacobianMatrix(
        Mat dRMat,
        const Vec preallocOnProc,
        const Vec preallocOffProc) const;
 
    void setupdRdWCon(
        const HashTable<List<List<word>>>& stateResConInfo,
        const label isPrealloc);
 
public:
    TypeName("dRdW");
    // Constructors
 
    //- Construct from components
    DAJacCondRdW(
        const word modelType,
        const fvMesh& mesh,
        const DAOption& daOption,
        const DAModel& daModel,
        const DAIndex& daIndex);
 
    //- Destructor
    virtual ~DAJacCondRdW()
    {
    }
 
    // Member functions
 
    virtual void setupJacConPreallocation(const dictionary& options);
 
    virtual void initializeJacCon(const dictionary& options);
 
    virtual void setupJacCon(const dictionary& options);
 
    virtual void preallocatedRdW(
        Mat dRMat,
        const label transposed) const;
    
    virtual void clear();
};
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
} // End namespace Foam
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#endif
 
// ************************************************************************* //