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