codac2_SlicedTube_operations.h

Go to the documentation of this file.

 
#pragma once
 
#include <type_traits>
 
#include "codac2_SlicedTube.h"
#include "codac2_math.h"
 
namespace codac2
{
  // These macros stand in for the body of the functions declared below. Doxygen
  // parses this header to produce the Python docstrings, and it cannot name a
  // function whose body it does not recognise, so every macro used that way is
  // listed in EXPAND_AS_DEFINED in doc/api/Doxyfile.in. Add any new one there
  // too, otherwise the bindings lose the docstring identifiers they expect and
  // stop compiling.
  #define macro_unary_tube(f) \
  { \
    auto y = x1; \
    for(auto it = y.tdomain()->begin() ; it != y.tdomain()->end() ; it++) \
    { \
      auto sy = y.slice(it); \
      sy->codomain() = f(sy->codomain()); \
    } \
    return y; \
  } \
 
  #define macro_binary_tube_tube(f,output_type) \
  { \
    assert_release(x1.tdomain() == x2.tdomain()); \
    auto y = output_type; \
    for(auto it = y.tdomain()->begin() ; it != y.tdomain()->end() ; it++) \
    { \
      auto sy = y.slice(it); \
      sy->codomain() = f(x1.slice(it)->codomain(), x2.slice(it)->codomain()); \
    } \
    return y; \
  } \
 
  #define macro_binary_real_tube(f) \
  { \
    auto y = x2; \
    for(auto it = y.tdomain()->begin() ; it != y.tdomain()->end() ; it++) \
    { \
      auto sy = y.slice(it); \
      sy->codomain() = f(x1, sy->codomain()); \
    } \
    return y; \
  } \
 
  #define macro_binary_tube_real(f) \
  { \
    auto y = x1; \
    for(auto it = y.tdomain()->begin() ; it != y.tdomain()->end() ; it++) \
    { \
      auto sy = y.slice(it); \
      sy->codomain() = f(sy->codomain(), x2); \
    } \
    return y; \
  } \
 
  #define macro_member_binary_tube_tube(f) \
  { \
    assert_release(x1.tdomain() == x2.tdomain()); \
    for(auto it = x1.tdomain()->begin() ; it != x1.tdomain()->end() ; it++) \
    { \
      auto sx1 = x1.slice(it); \
      sx1->codomain() = f(sx1->codomain(), x2.slice(it)->codomain()); \
    } \
    return x1; \
  } \
 
  #define macro_member_binary_tube_real(f) \
  { \
    for(auto it = x1.tdomain()->begin() ; it != x1.tdomain()->end() ; it++) \
    { \
      auto sx1 = x1.slice(it); \
      sx1->codomain() = f(sx1->codomain(), x2); \
    } \
    return x1; \
  } \
 
  template<typename U>
  struct is_slicedtube : std::false_type {};
 
  template<typename U>
  struct is_slicedtube<SlicedTube<U>> : std::true_type {};
 
  template<typename U>
  inline constexpr bool is_slicedtube_v = is_slicedtube<std::remove_cvref_t<U>>::value;
 
  template<typename Q>
  concept NonSlicedTube = !is_slicedtube_v<Q>;
 
  template<typename T>
  concept NonScalarTubeCodomain =
    !std::is_same_v<std::remove_cvref_t<T>, double> &&
    !std::is_same_v<std::remove_cvref_t<T>, Interval>;
 
  template<typename T>
  inline T operator_tube_union(const T& x1, const T& x2) { return x1 | x2; }
 
  template<typename T>
  inline T operator_tube_intersection(const T& x1, const T& x2) { return x1 & x2; }
 
  template<typename T, typename X1, typename X2>
  inline T operator_tube_add(const X1& x1, const X2& x2) { return x1 + x2; }
 
  template<typename T, typename X1, typename X2>
  inline T operator_tube_sub(const X1& x1, const X2& x2) { return x1 - x2; }
 
  template<typename T, typename X1, typename X2>
  inline T operator_tube_mul(const X1& x1, const X2& x2) { return x1 * x2; }
 
  template<typename T>
  inline T operator_tube_scal_mul(const Interval& x1, const T& x2) { return x1 * x2; }
 
  template<typename T>
  inline T operator_tube_mul_scal(const T& x1, const Interval& x2) { return x1 * x2; }
 
  inline IntervalVector operator_tube_mul_vec(const IntervalMatrix& x1, const IntervalVector& x2) { return x1 * x2; }
 
  template<typename T, typename X1, typename X2>
  inline T operator_tube_div(const X1& x1, const X2& x2) { return x1 / x2; }
 
  template<typename T, typename X1>
  inline T operator_tube_div_scal(const X1& x1, const Interval& x2) { return x1 / x2; }

  template<typename T>
  inline SlicedTube<T>& operator|=(SlicedTube<T>& x1, const SlicedTube<T>& x2)
    macro_member_binary_tube_tube(operator_tube_union<T>)
 
  
  template<typename T>
  inline SlicedTube<T> operator|(const SlicedTube<T>& x1, const SlicedTube<T>& x2)
    macro_binary_tube_tube(operator_tube_union<T>,x1)
 
  
  template<typename T>
  inline SlicedTube<T>& operator&=(SlicedTube<T>& x1, const SlicedTube<T>& x2)
    macro_member_binary_tube_tube(operator_tube_intersection<T>)
 
  
  template<typename T>
  inline SlicedTube<T> operator&(const SlicedTube<T>& x1, const SlicedTube<T>& x2)
    macro_binary_tube_tube(operator_tube_intersection<T>,x1)
 
  
  template<typename T>
  inline const SlicedTube<T>& operator+(const SlicedTube<T>& x1)
  {
    return x1;
  }

  template<typename T>
  inline SlicedTube<T> operator+(const SlicedTube<T>& x1, const SlicedTube<T>& x2)
    macro_binary_tube_tube(operator_tube_add<T>,x1)
 
  
  template<typename T, typename Q>
    requires NonSlicedTube<Q>
  inline SlicedTube<T> operator+(const SlicedTube<T>& x1, const Q& x2)
    macro_binary_tube_real(operator_tube_add<T>)
 
  
  template<typename T, typename Q>
    requires NonSlicedTube<Q>
  inline SlicedTube<T> operator+(const Q& x1, const SlicedTube<T>& x2)
    macro_binary_real_tube(operator_tube_add<T>)
 
  
  template<typename T, typename Q>
    requires NonSlicedTube<Q>
  inline SlicedTube<T>& operator+=(SlicedTube<T>& x1, const Q& x2)
    macro_member_binary_tube_real(operator_tube_add<T>)
 
  
  template<typename T>
  inline SlicedTube<T>& operator+=(SlicedTube<T>& x1, const SlicedTube<T>& x2)
    macro_member_binary_tube_tube(operator_tube_add<T>)
 
  
  template<typename T>
  inline SlicedTube<T> operator-(const SlicedTube<T>& x1)
  {
    return -1. * x1;
  }

  template<typename T>
  inline SlicedTube<T> operator-(const SlicedTube<T>& x1, const SlicedTube<T>& x2)
    macro_binary_tube_tube(operator_tube_sub<T>,x1)
 
  
  template<typename T, typename Q>
    requires NonSlicedTube<Q>
  inline SlicedTube<T> operator-(const SlicedTube<T>& x1, const Q& x2)
    macro_binary_tube_real(operator_tube_sub<T>)
 
  
  template<typename T, typename Q>
    requires NonSlicedTube<Q>
  inline SlicedTube<T> operator-(const Q& x1, const SlicedTube<T>& x2)
    macro_binary_real_tube(operator_tube_sub<T>)
 
  
  template<typename T, typename Q>
    requires NonSlicedTube<Q>
  inline SlicedTube<T>& operator-=(SlicedTube<T>& x1, const Q& x2)
    macro_member_binary_tube_real(operator_tube_sub<T>)
 
  
  template<typename T>
  inline SlicedTube<T>& operator-=(SlicedTube<T>& x1, const SlicedTube<T>& x2)
    macro_member_binary_tube_tube(operator_tube_sub<T>)
 
  
  template<typename T>
    requires (!std::is_same_v<T, double>)
  inline SlicedTube<T> operator*(const Interval& x1, const SlicedTube<T>& x2)
    macro_binary_real_tube(operator_tube_scal_mul<T>)
 
  
  template<typename T>
    requires (!std::is_same_v<T, double>)
  inline SlicedTube<T> operator*(const SlicedTube<T>& x1, const Interval& x2)
    macro_binary_tube_real(operator_tube_mul_scal<T>)
 
  
  inline SlicedTube<Interval> operator*(const SlicedTube<Interval>& x1, const SlicedTube<Interval>& x2)
    macro_binary_tube_tube(operator_tube_mul<Interval>,x1)

  inline SlicedTube<IntervalVector> operator*(const SlicedTube<IntervalVector>& x1, const SlicedTube<IntervalVector>& x2)
    macro_binary_tube_tube(operator_tube_mul<IntervalVector>,x1)

  inline SlicedTube<IntervalMatrix> operator*(const SlicedTube<IntervalMatrix>& x1, const SlicedTube<IntervalMatrix>& x2)
    macro_binary_tube_tube(operator_tube_mul<IntervalMatrix>,x1)

  template<typename T>
    requires NonScalarTubeCodomain<T>
  inline SlicedTube<T> operator*(const SlicedTube<Interval>& x1, const SlicedTube<T>& x2)
    macro_binary_tube_tube(operator_tube_mul<T>,x2)

  template<typename T>
    requires NonScalarTubeCodomain<T>
  inline SlicedTube<T> operator*(const SlicedTube<T>& x1, const SlicedTube<Interval>& x2)
    macro_binary_tube_tube(operator_tube_mul<T>,x1)

  template<typename T, typename Q>
    requires NonSlicedTube<Q>
  inline SlicedTube<T> operator*(const SlicedTube<T>& x1, const Q& x2)
    macro_binary_tube_real(operator_tube_mul<T>)

  template<typename T, typename Q>
    requires NonSlicedTube<Q>
  inline SlicedTube<T> operator*(const Q& x1, const SlicedTube<T>& x2)
    macro_binary_real_tube(operator_tube_mul<T>)

  inline SlicedTube<IntervalVector> operator*(const SlicedTube<IntervalMatrix>& x1, const SlicedTube<IntervalVector>& x2)
    macro_binary_tube_tube(operator_tube_mul_vec,x2)

  template<typename T, typename Q>
    requires NonSlicedTube<Q>
  inline SlicedTube<T>& operator*=(SlicedTube<T>& x1, const Q& x2)
    macro_member_binary_tube_real(operator_tube_mul<T>)

  template<typename T>
  inline SlicedTube<T>& operator*=(SlicedTube<T>& x1, const SlicedTube<T>& x2)
    macro_member_binary_tube_tube(operator_tube_mul<T>)

  template<typename T>
    requires NonScalarTubeCodomain<T>
  inline SlicedTube<T>& operator*=(SlicedTube<T>& x1, const SlicedTube<Interval>& x2)
    macro_member_binary_tube_tube(operator_tube_mul<T>)

  template<typename T>
    requires (!std::is_same_v<T, double>)
  inline SlicedTube<T> operator/(const SlicedTube<T>& x1, const Interval& x2)
    macro_binary_tube_real(operator_tube_div_scal<T>)

  inline SlicedTube<Interval> operator/(const SlicedTube<Interval>& x1, const SlicedTube<Interval>& x2)
    macro_binary_tube_tube(operator_tube_div<Interval>,x1)

  template<typename T>
    requires NonScalarTubeCodomain<T>
  inline SlicedTube<T> operator/(const SlicedTube<T>& x1, const SlicedTube<Interval>& x2)
    macro_binary_tube_tube(operator_tube_div<T>,x1)

  template<typename T, typename Q>
    requires NonSlicedTube<Q>
  inline SlicedTube<T> operator/(const SlicedTube<T>& x1, const Q& x2)
    macro_binary_tube_real(operator_tube_div<T>)

  template<typename T, typename Q>
    requires NonSlicedTube<Q>
  inline SlicedTube<T> operator/(const Q& x1, const SlicedTube<T>& x2)
    macro_binary_real_tube(operator_tube_div<T>)

  template<typename T, typename Q>
    requires NonSlicedTube<Q>
  inline SlicedTube<T>& operator/=(SlicedTube<T>& x1, const Q& x2)
    macro_member_binary_tube_real(operator_tube_div<T>)

  template<typename T>
  inline SlicedTube<T>& operator/=(SlicedTube<T>& x1, const SlicedTube<T>& x2)
    macro_member_binary_tube_tube(operator_tube_div<T>)

  template<typename T>
    requires NonScalarTubeCodomain<T>
  inline SlicedTube<T>& operator/=(SlicedTube<T>& x1, const SlicedTube<Interval>& x2)
    macro_member_binary_tube_tube(operator_tube_div<T>)

  SlicedTube<Interval> sqr(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> sqrt(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> pow(const SlicedTube<Interval>& x1, int x2);

  SlicedTube<Interval> pow(const SlicedTube<Interval>& x1, const Interval& x2);

  SlicedTube<Interval> exp(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> log(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> cos(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> sin(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> tan(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> acos(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> asin(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> atan(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> atan2(const SlicedTube<Interval>& x1, const SlicedTube<Interval>& x2);

  SlicedTube<Interval> atan2(const SlicedTube<Interval>& x1, const Interval& x2);

  SlicedTube<Interval> atan2(const Interval& x1, const SlicedTube<Interval>& x2);

  SlicedTube<Interval> cosh(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> sinh(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> tanh(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> acosh(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> asinh(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> atanh(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> abs(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> min(const SlicedTube<Interval>& x1, const SlicedTube<Interval>& x2);

  SlicedTube<Interval> min(const SlicedTube<Interval>& x1, const Interval& x2);

  SlicedTube<Interval> min(const Interval& x1, const SlicedTube<Interval>& x2);

  SlicedTube<Interval> max(const SlicedTube<Interval>& x1, const SlicedTube<Interval>& x2);

  SlicedTube<Interval> max(const SlicedTube<Interval>& x1, const Interval& x2);

  SlicedTube<Interval> max(const Interval& x1, const SlicedTube<Interval>& x2);

  SlicedTube<Interval> sign(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> integer(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> floor(const SlicedTube<Interval>& x1);

  SlicedTube<Interval> ceil(const SlicedTube<Interval>& x1);
}