codac2_CtcInter.h

Go to the documentation of this file.

 
#pragma once
 
#include <type_traits>
#include "codac2_CtcWrapper.h"
#include "codac2_Collection.h"
#include "codac2_template_tools.h"
 
namespace codac2
{
  template<typename... X>
  class CtcInter : public Ctc<CtcInter<X...>,X...>
  {
    public:

      explicit CtcInter(Index n)
        : Ctc<CtcInter<X...>,X...>(n)
      {
        if constexpr(std::is_same_v<X...,Interval>)
        {
          assert(n == 1);
        }
      }

      template<typename C>
        requires (IsCtcBaseOrPtr<C,X...> && !std::is_same_v<CtcInter<X...>,C>)
      CtcInter(const C& c)
        : Ctc<CtcInter<X...>,X...>(size_of(c)), _ctcs(c)
      { }

      template<typename... C>
        requires (IsCtcBaseOrPtr<C,X...> && ...)
      CtcInter(const C&... c)
        : Ctc<CtcInter<X...>,X...>(size_first_item(c...)), _ctcs(c...)
      {
        assert_release(all_same_size(c...));
      }

      CtcInter(const Collection<CtcBase<IntervalVector>>& ctcs)
        : Ctc<CtcInter<X...>,X...>(ctcs.front()->size()), _ctcs(ctcs)
      {
        for(const auto& ci : _ctcs)
        {
          assert_release(ci->size() == this->size());
        }
      }

      size_t nb() const
      {
        return _ctcs.size();
      }

      void contract(X&... x) const
      {
        for(const auto& ci : _ctcs)
        {
          ci->contract(x...);
          if((x.is_empty() | ...))
            return;
        }
      }

      template<typename C>
        requires IsCtcBaseOrPtr<C,X...>
      CtcInter<X...>& operator&=(const C& c)
      {
        assert_release(size_of(c) == this->size());
        _ctcs.push_back(c);
        return *this;
      }
 
    protected:

      Collection<CtcBase<X...>> _ctcs;
  };

  template <>
  class CtcInter<> : public CtcInter<IntervalVector>
  { };

  template<typename Tuple>
  struct CtcInterType;

  template<typename... T>
  struct CtcInterType<std::tuple<T...>> {
    using Ctc = CtcInter<T...>;
  };

  template<typename C1,typename C2>
  typename CtcInterType<typename C1::ContractedTypes>::Ctc operator&(const C1& c1, const C2& c2)
  {
    return { c1, c2 };
  }

  template<typename C1,typename C2>
  typename CtcInterType<typename C1::ContractedTypes>::Ctc operator&(const std::shared_ptr<C1>& c1, const std::shared_ptr<C2>& c2)
  {
    return { c1, c2 };
  }

  template<typename C1,typename C2>
  typename CtcInterType<typename C1::ContractedTypes>::Ctc operator&(const std::shared_ptr<C1>& c1, const C2& c2)
  {
    return { c1, c2 };
  }

  template<typename C1,typename C2>
  typename CtcInterType<typename C1::ContractedTypes>::Ctc operator&(const C1& c1, const std::shared_ptr<C2>& c2)
  {
    return { c1, c2 };
  }
  
  template<typename C2>
    requires std::is_base_of_v<CtcBase<IntervalVector>,C2>
  inline CtcInter<IntervalVector> operator&(const IntervalVector& c1, const C2& c2)
  {
    assert_release(c1.size() == c2.size());
    return CtcInter<IntervalVector>(CtcWrapper(c1),c2);
  }

  template<typename C1>
    requires std::is_base_of_v<CtcBase<IntervalVector>,C1>
  inline CtcInter<IntervalVector> operator&(const C1& c1, const IntervalVector& c2)
  {
    assert_release(c1.size() == c2.size());
    return CtcInter<IntervalVector>(c1,CtcWrapper(c2));
  }
 
  // Template deduction guides
 
  CtcInter(Index) -> CtcInter<IntervalVector>;
}