agenda_wrapper.h

Go to the documentation of this file.

#ifndef agenda_wrapper_h
#define agenda_wrapper_h

#include "oem.h"


class AgendaWrapper
{
public:

    const unsigned int m, n;
    OEMMatrixReference jacobian;
    OEMVector          yi;


    AgendaWrapper( Workspace * ws_,
                   unsigned int m_,
                   unsigned int n_,
                   Matrix & jacobian_,
                   Vector & yi_,
                   const Agenda *inversion_iterate_agenda_ ) :
        m(m_),
        n(n_),
        jacobian(jacobian_), yi(yi_), ws(ws_),
        inversion_iterate_agenda( inversion_iterate_agenda_ ),
        reuse_jacobian((jacobian_.nrows() != 0) &&
                        (jacobian_.ncols() != 0) &&
                       (yi_.nelem() != 0)),
        iteration_counter(0)
        {}

    AgendaWrapper(const AgendaWrapper &) = delete;
    AgendaWrapper(      AgendaWrapper &&) = delete;
    AgendaWrapper& operator=(const AgendaWrapper &)  = delete;
    AgendaWrapper& operator=(      AgendaWrapper &&) = delete;



    OEMMatrixReference Jacobian(const OEMVector & xi,
                                OEMVector & yi_)
    {
        if (!reuse_jacobian)
        {
            inversion_iterate_agendaExecute(
                *ws, yi, jacobian, xi, 1, 0,
                *inversion_iterate_agenda);
            yi_ = yi;
            iteration_counter += 1;
        } else {
            reuse_jacobian = false;
            yi_ = yi;
        }
        return jacobian;
    }


    OEMVector evaluate(const OEMVector &xi)
    {
        if (!reuse_jacobian)
        {
            Matrix dummy;
            inversion_iterate_agendaExecute(
                *ws, yi, dummy, xi, 0, iteration_counter,
                *inversion_iterate_agenda );
        } else {
            reuse_jacobian = false;
        }
        return yi;
    }

private:

    Workspace    * ws;
    const Agenda * inversion_iterate_agenda;
    bool           reuse_jacobian;
    unsigned int   iteration_counter;

};

#endif // agenda_wrappers_h