#ifndef X66AJaG3Qy6m9CIBTEKhiCYWag0
#define X66AJaG3Qy6m9CIBTEKhiCYWag0

#include "pla_evolve.hpp"



namespace EvolveImpl {

  template<typename Quant>
  struct Evolve {
    template<class PC>
    void operator() (PC &pc, const Time ct, const Ptr<Quant> id) {
      typedef typename QuantChild<Quant>::type Child;
      // first evolve childs
      for (Ptr<Child> i : id.childs())
	Evolve<Child>()(pc, ct, i);
      // then self
      pc.cast(PLA_Evolve<Quant>{ContinuousContext<Quant>{id}, ct});
    }
  };

  template<>
  struct Evolve<Synapse> {
    template<class PC>
    void operator() (PC &pc, const Time ct, const Ptr<Synapse> id) {
      // no childs -> evolve only self
      pc.cast(PLA_Evolve<Synapse>{ContinuousContext<Synapse>{id}, ct});
    }
  };
}

using EvolveImpl::Evolve;

#endif // X66AJaG3Qy6m9CIBTEKhiCYWag0