SR1 Class Reference

The quasi-Newton Symmetric Rank One formula. More...

Inheritance diagram for SR1:
[legend]
Collaboration diagram for SR1:
[legend]

Public Member Functions

 TypeName ("SR1")
 Runtime type information. More...
 
 SR1 (const fvMesh &mesh, const dictionary &dict)
 Construct from components. More...
 
virtual ~SR1 ()=default
 Destructor. More...
 
void computeCorrection ()
 Compute design variables correction. More...
 
virtual void updateOldCorrection (const scalarField &oldCorrection)
 
virtual void write ()
 Write old info to dict. More...
 
- Public Member Functions inherited from updateMethod
 TypeName ("updateMethod")
 Runtime type information. More...
 
 declareRunTimeSelectionTable (autoPtr, updateMethod, dictionary,(const fvMesh &mesh, const dictionary &dict),(mesh, dict))
 
 updateMethod (const fvMesh &mesh, const dictionary &dict)
 Construct from components. More...
 
virtual ~updateMethod ()=default
 Destructor. More...
 
void setObjectiveDeriv (const scalarField &derivs)
 Set objective derivative. More...
 
void setConstraintDeriv (const PtrList< scalarField > &derivs)
 Set constraints derivative. More...
 
void setObjectiveValue (const scalar value)
 Set constraints derivative. More...
 
void setConstraintValues (const scalarField &values)
 Set constraints derivative. More...
 
void setStep (const scalar eta)
 Set step for optimisation methods. More...
 
void setGlobalSum (const bool useGlobalSum)
 Set globalSum variable. More...
 
virtual void computeCorrection ()=0
 Return the correction of the design variables. More...
 
scalarFieldreturnCorrection ()
 Return the correction of the design variables. More...
 
void writeCorrection ()
 
virtual scalar computeMeritFunction ()
 
virtual scalar meritFunctionDirectionalDerivative ()
 
boolinitialEtaSet ()
 Return whether initial eta was set. More...
 
virtual void updateOldCorrection (const scalarField &oldCorrection)
 
virtual void write ()
 Write useful quantities to files. More...
 

Protected Member Functions

void allocateMatrices ()
 Allocate matrices in the first optimisation cycle. More...
 
void updateHessian ()
 Update approximation of the inverse Hessian. More...
 
void update ()
 Update design variables. More...
 
void readFromDict ()
 Read old info from dict. More...
 
- Protected Member Functions inherited from updateMethod
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)
 
scalar globalSum (const scalarField &field)
 Compute either global or local sum, based on globalSum flag. More...
 
scalar globalSum (tmp< scalarField > &tfield)
 Compute either global or local sum, based on globalSum flag. More...
 
dictionary coeffsDict ()
 Return optional dictionary with parameters specific to each method. More...
 

Protected Attributes

scalar etaHessian_
 Step for the Newton method. More...
 
label nSteepestDescent_
 Number of first steepest descent steps. More...
 
scalar ratioThreshold_
 For stability check. More...
 
labelList activeDesignVars_
 Map to active design variables. More...
 
SquareMatrix< scalar > HessianInv_
 
SquareMatrix< scalar > HessianInvOld_
 The previous Hessian inverse. More...
 
scalarField derivativesOld_
 The previous derivatives. More...
 
scalarField correctionOld_
 The previous correction. More...
 
label counter_
 Optimisation counter. More...
 
- Protected Attributes inherited from updateMethod
const fvMeshmesh_
 
const dictionary dict_
 
IOdictionary optMethodIODict_
 Used to output values useful for continuation runs. More...
 
scalarField objectiveDerivatives_
 Derivatives of the objective functions. More...
 
PtrList< scalarFieldconstraintDerivatives_
 Derivatives of the constraints. More...
 
scalar objectiveValue_
 Objective value. More...
 
scalarField cValues_
 Constraint values. More...
 
scalarField correction_
 Design variables correction. More...
 
scalarField cumulativeCorrection_
 
scalar eta_
 Step multiplying the correction. More...
 
bool initialEtaSet_
 Is initially set? More...
 
word correctionFolder_
 Folder storing the corrections to file. More...
 
bool globalSum_
 Whether to use gSum or sum in the inner products. More...
 

Additional Inherited Members

- Static Public Member Functions inherited from updateMethod
static autoPtr< updateMethodNew (const fvMesh &mesh, const dictionary &dict)
 Return a reference to the selected turbulence model. More...
 

Detailed Description

The quasi-Newton Symmetric Rank One formula.

Source files

Definition at line 54 of file SR1.H.

Constructor & Destructor Documentation

◆ SR1()

SR1 ( const fvMesh mesh,
const dictionary dict 
)

Construct from components.

Definition at line 159 of file SR1.C.

References SR1::activeDesignVars_, updateMethod::coeffsDict(), Foam::endl(), Foam::Info, SR1::readFromDict(), and readIfPresent().

Here is the call graph for this function:

◆ ~SR1()

virtual ~SR1 ( )
virtualdefault

Destructor.

Member Function Documentation

◆ allocateMatrices()

void allocateMatrices ( )
protected

Allocate matrices in the first optimisation cycle.

Definition at line 49 of file SR1.C.

References SR1::activeDesignVars_, UList< T >::empty(), forAll, SR1::HessianInv_, SR1::HessianInvOld_, Foam::identity(), updateMethod::objectiveDerivatives_, UList< T >::size(), and Foam::Zero.

Here is the call graph for this function:

◆ updateHessian()

void updateHessian ( )
protected

Update approximation of the inverse Hessian.

Definition at line 71 of file SR1.C.

References Foam::endl(), s(), Foam::sqr(), Foam::sqrt(), WarningInFunction, y, and Foam::Zero.

Here is the call graph for this function:

◆ update()

void update ( )
protected

Update design variables.

Definition at line 102 of file SR1.C.

References Foam::endl(), forAll, Foam::Info, Field< Type >::map(), and Foam::Zero.

Here is the call graph for this function:

◆ readFromDict()

void readFromDict ( )
protected

Read old info from dict.

Definition at line 135 of file SR1.C.

References Foam::identity(), n, and Foam::Zero.

Referenced by SR1::SR1().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ TypeName()

TypeName ( "SR1"  )

Runtime type information.

◆ computeCorrection()

void computeCorrection ( )
virtual

Compute design variables correction.

Implements updateMethod.

Definition at line 201 of file SR1.C.

References update().

Here is the call graph for this function:

◆ updateOldCorrection()

void updateOldCorrection ( const scalarField oldCorrection)
virtual

Update old correction. Useful for quasi-Newton methods coupled with line search

Reimplemented from updateMethod.

Definition at line 217 of file SR1.C.

References updateMethod::updateOldCorrection().

Here is the call graph for this function:

◆ write()

void write ( )
virtual

Write old info to dict.

Reimplemented from updateMethod.

Definition at line 224 of file SR1.C.

References updateMethod::write().

Here is the call graph for this function:

Member Data Documentation

◆ etaHessian_

scalar etaHessian_
protected

Step for the Newton method.

Definition at line 63 of file SR1.H.

◆ nSteepestDescent_

label nSteepestDescent_
protected

Number of first steepest descent steps.

Definition at line 66 of file SR1.H.

◆ ratioThreshold_

scalar ratioThreshold_
protected

For stability check.

Definition at line 69 of file SR1.H.

◆ activeDesignVars_

labelList activeDesignVars_
protected

Map to active design variables.

Definition at line 72 of file SR1.H.

Referenced by SR1::allocateMatrices(), and SR1::SR1().

◆ HessianInv_

SquareMatrix<scalar> HessianInv_
protected

The Hessian inverse. Should have the size of the active design variables

Definition at line 76 of file SR1.H.

Referenced by SR1::allocateMatrices().

◆ HessianInvOld_

SquareMatrix<scalar> HessianInvOld_
protected

The previous Hessian inverse.

Definition at line 79 of file SR1.H.

Referenced by SR1::allocateMatrices().

◆ derivativesOld_

scalarField derivativesOld_
protected

The previous derivatives.

Definition at line 82 of file SR1.H.

◆ correctionOld_

scalarField correctionOld_
protected

The previous correction.

Definition at line 85 of file SR1.H.

◆ counter_

label counter_
protected

Optimisation counter.

Definition at line 88 of file SR1.H.


The documentation for this class was generated from the following files: