codac2_cart_prod.h

Go to the documentation of this file.

 
#pragma once
 
#include "codac2_IntervalVector.h"
#include "codac2_Matrix.h"
#include "codac2_matrices.h"
 
namespace codac2
{
  inline IntervalVector to_IntervalVector(const Interval& x)
  {
    return { x };
  }
 
  inline IntervalVector to_IntervalVector(const IntervalVector& x)
  {
    return x;
  }
 
  inline IntervalVector to_IntervalVector(const Vector& x)
  {
    return x.template cast<Interval>();
  }
 
  // Flatten some matrix into an interval vector
  template<typename T,int R,int C>
    requires (!Eigen::IsVectorOrRow<R,C>)
  inline IntervalVector to_IntervalVector(const Eigen::Matrix<T,R,C>& x)
  {
    return Eigen::Map<const Eigen::Matrix<T,R,1>>(x.data(), x.size());
  }
 
  inline IntervalVector cart_prod()
  {
    return IntervalVector(0);
  }
 
  template<typename... X>
    requires ((is_interval_based_v<X>) || ...)
      && ((!is_ctc_v<X>) && ...) && ((!is_sep_v<X>) && ...)
  inline IntervalVector cart_prod(const X&... x)
  {
    std::vector<IntervalVector> v_x;
    ((v_x.push_back(to_IntervalVector(x))), ...);
 
    Index n = 0;
    for(const auto& xi : v_x)
      n += xi.size();
    IntervalVector x_(n);
 
    Index i = 0;
    for(const auto& xi : v_x)
    {
      x_.put(i, xi);
      i += xi.size();
    }
    return x_;
  }
 
  template<typename... X>
    requires ((!is_interval_based_v<X>) && ...)
      && ((!is_ctc_v<X>) && ...) && ((!is_sep_v<X>) && ...)
  inline Vector cart_prod(const X&... x)
  {
    return cart_prod(to_IntervalVector(x)...).mid();
  }
}