updateMethod.H

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
    \\  /    A nd           | www.openfoam.com
     \\/     M anipulation  |
-------------------------------------------------------------------------------
    Copyright (C) 2007-2023 PCOpt/NTUA
    Copyright (C) 2013-2023 FOSS GP
    Copyright (C) 2019 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
    This file is part of OpenFOAM.

    OpenFOAM is free software: you can redistribute it and/or modify it
    under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    OpenFOAM 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.  See the GNU General Public License
    for more details.

    You should have received a copy of the GNU General Public License
    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.

Class
    Foam::updateMethod

Description
    Abstract base class for optimisation methods

SourceFiles
    updateMethod.C

\*---------------------------------------------------------------------------*/

#ifndef updateMethod_H
#define updateMethod_H

#include "localIOdictionary.H"
#include "designVariables.H"
#include "runTimeSelectionTables.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

namespace Foam
{

/*---------------------------------------------------------------------------*\
                      Class updateMethod Declaration
\*---------------------------------------------------------------------------*/

class updateMethod
:
    public localIOdictionary
{
protected:

    // Protected data

        const fvMesh& mesh_;

        const dictionary dict_;

        //- Design variables
        autoPtr<designVariables>& designVars_;

        //- Number of constraints
        label nConstraints_;

        //- Map to active design variables
        const labelList& activeDesignVars_;

        //- Derivatives of the objective functions
        scalarField objectiveDerivatives_;

        //- Derivatives of the constraints
        PtrList<scalarField> constraintDerivatives_;

        //- Objective value
        scalar objectiveValue_;

        //- Old objective value
        //  Used for convergence checking
        scalar objectiveValueOld_;

        //- Constraint values
        scalarField cValues_;

        //- Design variables correction
        scalarField correction_;

        //- Cumulative design variables correction throughout the optimisation
        //- loop
        scalarField cumulativeCorrection_;

        //- Step multiplying the correction
        scalar eta_;

        //- Optimisation cycle count
        label counter_;

        //- Is initially set?
        bool initialEtaSet_;

        //- Folder storing the corrections to file
        //  For some optimisation methods with a very high number of
        //  design variables (e.g. topology), it doesn't make much sense
        //  to write all updates in the updateMethodDict. Hence, a
        //  separate file is used to write the corrections, in case they are
        //  needed for post-processing
        word correctionFolder_;

        //- Whether to use gSum or sum in the inner products
        bool globalSum_;

        // Scalar -- matrix multiplications
        const scalarField leftMult
        (
            const scalarField&,
            const SquareMatrix<scalar>&
        );

        const scalarField rightMult
        (
            const SquareMatrix<scalar>&,
            const scalarField&
        );

        SquareMatrix<scalar> outerProd
        (
            const scalarField&,
            const scalarField&
        );

        SquareMatrix<scalar> inv(SquareMatrix<scalar> A);

        //- Compute either global or local sum, based on globalSum flag
        scalar globalSum(const scalarField& field);

        //- Compute either global or local sum, based on globalSum flag
        scalar globalSum(tmp<scalarField>& tfield);

        //- Compute either global or local sum, based on globalSum flag
        label globalSum(const label);

        //- Helper function to either read a scalarField of certain size
        //- from a dictionary, or construct a zero field
        tmp<scalarField> readOrZeroField(const word& name, const label size);


private:

    // Private Member Functions

        //- No copy construct
        updateMethod(const updateMethod&) = delete;

        //- No copy assignment
        void operator=(const updateMethod&) = delete;


public:

    //- Runtime type information
    TypeName("updateMethod");


    // Declare run-time constructor selection table

        declareRunTimeSelectionTable
        (
            autoPtr,
            updateMethod,
            dictionary,
            (
                const fvMesh& mesh,
                const dictionary& dict,
                autoPtr<designVariables>& designVars,
                const label nConstraints,
                const word& type
            ),
            (mesh, dict, designVars, nConstraints, type)
        );


    // Constructors

        //- Construct from components
        updateMethod
        (
            const fvMesh& mesh,
            const dictionary& dict,
            autoPtr<designVariables>& designVars,
            const label nConstraints,
            const word& type
        );


    // Selectors

        //- Return a reference to the selected turbulence model
        static autoPtr<updateMethod> New
        (
            const fvMesh& mesh,
            const dictionary& dict,
            autoPtr<designVariables>& designVars,
            const label nConstraints
        );


    //- Destructor
    virtual ~updateMethod() = default;


    // Member Functions

        //- Return optional dictionary with parameters specific to each method
        dictionary coeffsDict(const word& type) const;

        //- Set objective derivative
        void setObjectiveDeriv(const scalarField& derivs);

        //- Set constraints derivative
        void setConstraintDeriv(const PtrList<scalarField>& derivs);

        //- Set objective value
        void setObjectiveValue(const scalar value);

        //- Set old objective value
        void setObjectiveValueOld(const scalar value);

        //- Set values of constraints
        void setConstraintValues(const scalarField& values);

        //- Get objective value
        scalar getObjectiveValue() const;

        //- Get old objective value
        scalar getObjectiveValueOld() const;

        //- Get values of constraints
        const scalarField& getConstraintValues() const;

        //- Get optimisation cycle
        label getCycle() const;

        //- Set step for optimisation methods
        void setStep(const scalar eta);

        //- Multiply step
        void modifyStep(const scalar multiplier);

        //- Set globalSum variable.
        //  Should be set by the optimisationManager owning the updateMethod
        void setGlobalSum(const bool useGlobalSum);

        //- Set the number of constraints
        void setConstaintsNumber(const label nConstraints);

        //- Get the number of constraints
        label nConstraints() const;

        //- Return the correction of the design variables
        virtual void computeCorrection() = 0;

        //- Return the correction of the design variables
        //const scalarField& returnCorrection() const;

        //- Return the correction of the design variables
        scalarField& returnCorrection();

        void writeCorrection();

        //- Compute merit function. Could be different than the objective
        //- in the presence of constraints
        virtual scalar computeMeritFunction();

        //- Directional derivative of the merit function, in the direction of
        //- the correction. Could be different than the objective directional
        //- derivative in the presence of constraints
        virtual scalar meritFunctionDirectionalDerivative();

        //- Return whether initial eta was set
        bool& initialEtaSet();

        //- Update old correction. Useful for quasi-newton methods coupled with
        //- line search
        virtual void updateOldCorrection(const scalarField& oldCorrection);

        //- Write continuation data under uniform
        virtual bool writeData(Ostream& os) const;

        //- Write non-continuation data, usually under the optimisation folder
        virtual bool writeAuxiliaryData();
};


// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

} // End namespace Foam

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

#endif

// ************************************************************************* //