codac_CtcPicard.h

#ifndef __CODAC_CTCPICARD_H__
#define __CODAC_CTCPICARD_H__

#include "codac_DynCtc.h"
#include "codac_TFnc.h"
#include "codac_Slice.h"

namespace codac
{
  class CtcPicard : public DynCtc
  {
    public:

      CtcPicard(const Function& f, float delta = 1.1);
      CtcPicard(const TFnc& f, float delta = 1.1);
      ~CtcPicard();

      void contract(std::vector<Domain*>& v_domains);
      void contract(Tube& x, TimePropag t_propa = TimePropag::FORWARD | TimePropag::BACKWARD);
      void contract(TubeVector& x, TimePropag t_propa = TimePropag::FORWARD | TimePropag::BACKWARD);

      int picard_iterations() const;

    protected:

      void contract_kth_slices(TubeVector& x, int k, TimePropag t_propa);
      void guess_kth_slices_envelope(TubeVector& x, int k, TimePropag t_propa);

      const TFunction* m_f_ptr = nullptr;
      const TFnc& m_f;
      const float m_delta;
      int m_picard_iterations = 0;

      static const std::string m_ctc_name; 
      static std::vector<std::string> m_str_expected_doms; 
      friend class ContractorNetwork;
  };
}

#endif