codac2_tube_cart_prod.h

Go to the documentation of this file.

 
#pragma once
 
#include <tuple>
#include "codac2_cart_prod.h"
#include "codac2_SlicedTube.h"
 
namespace codac2
{
  template<typename... X>
  SlicedTube<IntervalVector> tube_cart_prod(const X&... x)
  {
    static_assert(sizeof...(X) > 0, "tube_cart_prod requires at least one tube");
    Index s = 0;
    ((s += x.size()), ...);
 
    const auto& x0 = std::get<0>(std::forward_as_tuple(x...));
    auto tdomain = x0.tdomain();
    assert_release(((x.tdomain() == tdomain) && ...));
    SlicedTube v(tdomain, IntervalVector(s));
 
    for(auto it = tdomain->begin() ; it != tdomain->end() ; it++)
      v.slice(it)->codomain() = cart_prod(x.slice(it)->codomain()...);
 
    return v;
  }
}