OOPS
Public Member Functions | Private Types | Private Member Functions | Private Attributes | List of all members
oops::SaddlePointLMPMatrix< MODEL, OBS > Class Template Reference

The preconditioner for the saddle-point minimizer. More...

#include <SaddlePointLMPMatrix.h>

Inheritance diagram for oops::SaddlePointLMPMatrix< MODEL, OBS >:
Inheritance graph
[legend]
Collaboration diagram for oops::SaddlePointLMPMatrix< MODEL, OBS >:
Collaboration graph
[legend]

Public Member Functions

 SaddlePointLMPMatrix (const CostFct_ &j)
 
void setup (const std::vector< SPVector_ > &, const std::vector< SPVector_ > &)
 
const int & getk () const
 
void multiply (const SPVector_ &, SPVector_ &) const
 

Private Types

typedef ControlIncrement< MODEL, OBS > CtrlInc_
 
typedef CostFctWeak< MODEL, OBS > CostFctWeak_
 
typedef CostFunction< MODEL, OBS > CostFct_
 
typedef SaddlePointVector< MODEL, OBS > SPVector_
 
typedef DualVector< MODEL, OBS > LagVector_
 

Private Member Functions

void Pinitmultiply (const SPVector_ &, SPVector_ &) const
 
void Gmultiply (const SPVector_ &, Eigen::VectorXd &) const
 
void Rmultiply (Eigen::VectorXd &, SPVector_ &) const
 
void Fmultiply (Eigen::VectorXd &, Eigen::VectorXd &) const
 

Private Attributes

const CostFctWeak_ & j_
 
const bool idmodel_
 
std::vector< SPVector_ > xyVEC_
 
std::vector< SPVector_ > pqVEC_
 
std::vector< LagVector_ > RpVEC_
 
std::vector< CtrlInc_ > RqVEC_
 
std::unique_ptr< SPVector_ > spvecinit_
 
int nvec_
 
Eigen::MatrixXd ZMat_
 
Eigen::MatrixXd FMat_
 

Detailed Description

template<typename MODEL, typename OBS>
class oops::SaddlePointLMPMatrix< MODEL, OBS >

The preconditioner for the saddle-point minimizer.

The preconditioner is obtained by using low-rank updates. Let us define the matrices:

\( R = [ 0 Rp ], G = [ Zinv*Rp' 0 ] \) \( [ Rq 0 ] [ 0 Zinv'*Rq' ] \)

\( F = I + G*P0inv*R \)

where \( P0 \) is the initial preconditioner and \( I \) is the identity matrix of order 2.

Then the preconditioner is updated from

\( Pk = P0inv - P0inv*R*Finv*G*P0inv \)

The solvers represent matrices as objects that implement a "multiply" method. This class defines objects that apply the saddle-point matrix.

Definition at line 55 of file SaddlePointLMPMatrix.h.

Member Typedef Documentation

◆ CostFct_

template<typename MODEL , typename OBS >
typedef CostFunction<MODEL, OBS> oops::SaddlePointLMPMatrix< MODEL, OBS >::CostFct_
private

Definition at line 58 of file SaddlePointLMPMatrix.h.

◆ CostFctWeak_

template<typename MODEL , typename OBS >
typedef CostFctWeak<MODEL, OBS> oops::SaddlePointLMPMatrix< MODEL, OBS >::CostFctWeak_
private

Definition at line 57 of file SaddlePointLMPMatrix.h.

◆ CtrlInc_

template<typename MODEL , typename OBS >
typedef ControlIncrement<MODEL, OBS> oops::SaddlePointLMPMatrix< MODEL, OBS >::CtrlInc_
private

Definition at line 56 of file SaddlePointLMPMatrix.h.

◆ LagVector_

template<typename MODEL , typename OBS >
typedef DualVector<MODEL, OBS> oops::SaddlePointLMPMatrix< MODEL, OBS >::LagVector_
private

Definition at line 60 of file SaddlePointLMPMatrix.h.

◆ SPVector_

template<typename MODEL , typename OBS >
typedef SaddlePointVector<MODEL, OBS> oops::SaddlePointLMPMatrix< MODEL, OBS >::SPVector_
private

Definition at line 59 of file SaddlePointLMPMatrix.h.

Constructor & Destructor Documentation

◆ SaddlePointLMPMatrix()

template<typename MODEL , typename OBS >
oops::SaddlePointLMPMatrix< MODEL, OBS >::SaddlePointLMPMatrix ( const CostFct_ &  j)
explicit

Definition at line 95 of file SaddlePointLMPMatrix.h.

Member Function Documentation

◆ Fmultiply()

template<typename MODEL , typename OBS >
void oops::SaddlePointLMPMatrix< MODEL, OBS >::Fmultiply ( Eigen::VectorXd &  x,
Eigen::VectorXd &  fx 
) const
private

Definition at line 272 of file SaddlePointLMPMatrix.h.

Here is the call graph for this function:

◆ getk()

template<typename MODEL , typename OBS >
const int& oops::SaddlePointLMPMatrix< MODEL, OBS >::getk ( ) const
inline

Definition at line 67 of file SaddlePointLMPMatrix.h.

◆ Gmultiply()

template<typename MODEL , typename OBS >
void oops::SaddlePointLMPMatrix< MODEL, OBS >::Gmultiply ( const SPVector_ &  x,
Eigen::VectorXd &  z 
) const
private

Definition at line 234 of file SaddlePointLMPMatrix.h.

Here is the call graph for this function:

◆ multiply()

template<typename MODEL , typename OBS >
void oops::SaddlePointLMPMatrix< MODEL, OBS >::multiply ( const SPVector_ &  x,
SPVector_ &  z 
) const

Definition at line 212 of file SaddlePointLMPMatrix.h.

Here is the call graph for this function:

◆ Pinitmultiply()

template<typename MODEL , typename OBS >
void oops::SaddlePointLMPMatrix< MODEL, OBS >::Pinitmultiply ( const SPVector_ &  x,
SPVector_ &  z 
) const
private

Definition at line 160 of file SaddlePointLMPMatrix.h.

Here is the call graph for this function:

◆ Rmultiply()

template<typename MODEL , typename OBS >
void oops::SaddlePointLMPMatrix< MODEL, OBS >::Rmultiply ( Eigen::VectorXd &  x,
SPVector_ &  rx 
) const
private

Definition at line 256 of file SaddlePointLMPMatrix.h.

Here is the call graph for this function:

◆ setup()

template<typename MODEL , typename OBS >
void oops::SaddlePointLMPMatrix< MODEL, OBS >::setup ( const std::vector< SPVector_ > &  xyVEC,
const std::vector< SPVector_ > &  pqVEC 
)

Definition at line 103 of file SaddlePointLMPMatrix.h.

Here is the call graph for this function:

Member Data Documentation

◆ FMat_

template<typename MODEL , typename OBS >
Eigen::MatrixXd oops::SaddlePointLMPMatrix< MODEL, OBS >::FMat_
private

Definition at line 89 of file SaddlePointLMPMatrix.h.

◆ idmodel_

template<typename MODEL , typename OBS >
const bool oops::SaddlePointLMPMatrix< MODEL, OBS >::idmodel_
private

Definition at line 81 of file SaddlePointLMPMatrix.h.

◆ j_

template<typename MODEL , typename OBS >
const CostFctWeak_& oops::SaddlePointLMPMatrix< MODEL, OBS >::j_
private

Definition at line 80 of file SaddlePointLMPMatrix.h.

◆ nvec_

template<typename MODEL , typename OBS >
int oops::SaddlePointLMPMatrix< MODEL, OBS >::nvec_
private

Definition at line 87 of file SaddlePointLMPMatrix.h.

◆ pqVEC_

template<typename MODEL , typename OBS >
std::vector<SPVector_> oops::SaddlePointLMPMatrix< MODEL, OBS >::pqVEC_
private

Definition at line 83 of file SaddlePointLMPMatrix.h.

◆ RpVEC_

template<typename MODEL , typename OBS >
std::vector<LagVector_> oops::SaddlePointLMPMatrix< MODEL, OBS >::RpVEC_
private

Definition at line 84 of file SaddlePointLMPMatrix.h.

◆ RqVEC_

template<typename MODEL , typename OBS >
std::vector<CtrlInc_> oops::SaddlePointLMPMatrix< MODEL, OBS >::RqVEC_
private

Definition at line 85 of file SaddlePointLMPMatrix.h.

◆ spvecinit_

template<typename MODEL , typename OBS >
std::unique_ptr<SPVector_> oops::SaddlePointLMPMatrix< MODEL, OBS >::spvecinit_
private

Definition at line 86 of file SaddlePointLMPMatrix.h.

◆ xyVEC_

template<typename MODEL , typename OBS >
std::vector<SPVector_> oops::SaddlePointLMPMatrix< MODEL, OBS >::xyVEC_
private

Definition at line 82 of file SaddlePointLMPMatrix.h.

◆ ZMat_

template<typename MODEL , typename OBS >
Eigen::MatrixXd oops::SaddlePointLMPMatrix< MODEL, OBS >::ZMat_
private

Definition at line 88 of file SaddlePointLMPMatrix.h.

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