generic_expression.hpp

/*
 *    This file is part of CasADi.
 *
 *    CasADi -- A symbolic framework for dynamic optimization.
 *    Copyright (C) 2010-2023 Joel Andersson, Joris Gillis, Moritz Diehl,
 *                            KU Leuven. All rights reserved.
 *    Copyright (C) 2011-2014 Greg Horn
 *
 *    CasADi is free software; you can redistribute it and/or
 *    modify it under the terms of the GNU Lesser General Public
 *    License as published by the Free Software Foundation; either
 *    version 3 of the License, or (at your option) any later version.
 *
 *    CasADi is distributed in the hope that it will be useful,
 *    but WITHOUT ANY WARRANTY; without even the implied warranty of
 *    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 *    Lesser General Public License for more details.
 *
 *    You should have received a copy of the GNU Lesser General Public
 *    License along with CasADi; if not, write to the Free Software
 *    Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
 *
 */
 
 
#ifndef CASADI_GENERIC_EXPRESSION_HPP
#define CASADI_GENERIC_EXPRESSION_HPP
 
#include "calculus.hpp"
 
namespace casadi {
 
  struct CASADI_EXPORT GenericExpressionCommon {};
 
#ifndef SWIG
template<typename ExType>
class GenericExpression : public GenericExpressionCommon {
  protected:
    // Helper functions
    inline const ExType& self() const { return static_cast<const ExType&>(*this); }
    inline ExType& self() { return static_cast<ExType&>(*this); }
  public:
 
 
  static ExType plus(const ExType &x, const ExType &y) {
    return ExType::binary(OP_ADD, x, y);
  }
  friend inline ExType plus(const ExType &x, const ExType &y) {
    return ExType::plus(x, y);
  }
  friend inline ExType operator+(const ExType &x, const ExType &y) {
    return plus(x, y);
  }
  inline ExType& operator+=(const ExType &y) { return self() = self() + y; }
 
 
  static ExType minus(const ExType &x, const ExType &y) {
    return ExType::binary(OP_SUB, x, y);
  }
  friend inline ExType minus(const ExType &x, const ExType &y) {
    return ExType::minus(x, y);
  }
  friend inline ExType operator-(const ExType &x, const ExType &y) {
    return minus(x, y);
  }
  inline ExType& operator-=(const ExType &y) { return self() = self() - y; }
 
 
  static ExType times(const ExType &x, const ExType &y) {
    return ExType::binary(OP_MUL, x, y);
  }
  friend inline ExType times(const ExType &x, const ExType &y) {
    return ExType::times(x, y);
  }
  friend inline ExType operator*(const ExType &x, const ExType &y) {
    return times(x, y);
  }
  inline ExType& operator*=(const ExType &y) {return self() = self() * y;}
 
 
  static ExType rdivide(const ExType &x, const ExType &y) {
    return ExType::binary(OP_DIV, x, y);
  }
  friend inline ExType rdivide(const ExType &x, const ExType &y) {
    return ExType::rdivide(x, y);
  }
  friend inline ExType operator/(const ExType &x, const ExType &y) {
    return rdivide(x, y);
  }
  inline ExType& operator/=(const ExType &y) {return self() = self() / y;}
 
 
  static ExType lt(const ExType &x, const ExType &y) {
    return ExType::binary(OP_LT, x, y);
  }
  friend inline ExType lt(const ExType &x, const ExType &y) {
    return ExType::lt(x, y);
  }
  friend inline ExType operator<(const ExType &x, const ExType &y) {
    return lt(x, y);
  }
 
 
  static ExType le(const ExType &x, const ExType &y) {
    return ExType::binary(OP_LE, x, y);
  }
  friend inline ExType le(const ExType &x, const ExType &y) {
    return ExType::le(x, y);
  }
  friend inline ExType operator<=(const ExType &x, const ExType &y) {
    return le(x, y);
  }
 
 
  static ExType gt(const ExType &x, const ExType &y) {
    return ExType::lt(y, x);
  }
  friend inline ExType gt(const ExType &x, const ExType &y) {
    return ExType::gt(x, y);
  }
  friend inline ExType operator>(const ExType &x, const ExType &y) {
    return gt(x, y);
  }
 
 
  static ExType ge(const ExType &x, const ExType &y) {
    return ExType::le(y, x);
  }
  friend inline ExType ge(const ExType &x, const ExType &y) {
    return ExType::ge(x, y);
  }
  friend inline ExType operator>=(const ExType &x, const ExType &y) {
    return ge(x, y);
  }
 
 
  static ExType eq(const ExType &x, const ExType &y) {
    return ExType::binary(OP_EQ, x, y);
  }
  friend inline ExType eq(const ExType &x, const ExType &y) {
    return ExType::eq(x, y);
  }
  friend inline ExType operator==(const ExType &x, const ExType &y) {
    return eq(x, y);
  }
 
 
  static ExType ne(const ExType &x, const ExType &y) {
    return ExType::binary(OP_NE, x, y);
  }
  friend inline ExType ne(const ExType &x, const ExType &y) {
    return ExType::ne(x, y);
  }
  friend inline ExType operator!=(const ExType &x, const ExType &y) {
    return ne(x, y);
  }
 
 
   static ExType logic_and(const ExType &x, const ExType &y) {
     return ExType::binary(OP_AND, x, y);
   }
   friend inline ExType logic_and(const ExType &x, const ExType &y) {
     return ExType::logic_and(x, y);
   }
   friend inline ExType operator&&(const ExType &x, const ExType &y) {
     return logic_and(x, y);
   }
 
 
   static ExType logic_or(const ExType &x, const ExType &y) {
     return ExType::binary(OP_OR, x, y);
   }
   friend inline ExType logic_or(const ExType &x, const ExType &y) {
     return ExType::logic_or(x, y);
   }
   friend inline ExType operator||(const ExType &x, const ExType &y) {
     return logic_or(x, y);
   }
 
 
    static ExType logic_not(const ExType& x) {
      return ExType::unary(OP_NOT, x);
    }
    friend inline ExType logic_not(const ExType& x) {
      return ExType::logic_not(x);
    }
    inline ExType operator!() const {
      return logic_not(self());
    }
 
 
    static ExType abs(const ExType& x) {
      return ExType::unary(OP_FABS, x);
    }
    friend inline ExType abs(const ExType& x) {
      return ExType::abs(x);
    }
    friend inline ExType fabs(const ExType& x) {
      return abs(x);
    }
 
 
    static ExType sqrt(const ExType& x) {
      return ExType::unary(OP_SQRT, x);
    }
    friend inline ExType sqrt(const ExType& x) {
      return ExType::sqrt(x);
    }
 
 
    static ExType sq(const ExType& x) {
      return ExType::unary(OP_SQ, x);
    }
    friend inline ExType sq(const ExType& x) {
      return ExType::sq(x);
    }
 
 
    static ExType sin(const ExType& x) {
      return ExType::unary(OP_SIN, x);
    }
    friend inline ExType sin(const ExType& x) {
      return ExType::sin(x);
    }
 
 
    static ExType cos(const ExType& x) {
      return ExType::unary(OP_COS, x);
    }
    friend inline ExType cos(const ExType& x) {
      return ExType::cos(x);
    }
 
 
    static ExType tan(const ExType& x) {
      return ExType::unary(OP_TAN, x);
    }
    friend inline ExType tan(const ExType& x) {
      return ExType::tan(x);
    }
 
 
    static ExType atan(const ExType& x) {
      return ExType::unary(OP_ATAN, x);
    }
    friend inline ExType atan(const ExType& x) {
      return ExType::atan(x);
    }
 
 
    static ExType asin(const ExType& x) {
      return ExType::unary(OP_ASIN, x);
    }
    friend inline ExType asin(const ExType& x) {
      return ExType::asin(x);
    }
 
 
    static ExType acos(const ExType& x) {
      return ExType::unary(OP_ACOS, x);
    }
    friend inline ExType acos(const ExType& x) {
      return ExType::acos(x);
    }
 
 
    static ExType tanh(const ExType& x) {
      return ExType::unary(OP_TANH, x);
    }
    friend inline ExType tanh(const ExType& x) {
      return ExType::tanh(x);
    }
 
 
    static ExType sinh(const ExType& x) {
      return ExType::unary(OP_SINH, x);
    }
    friend inline ExType sinh(const ExType& x) {
      return ExType::sinh(x);
    }
 
 
    static ExType cosh(const ExType& x) {
      return ExType::unary(OP_COSH, x);
    }
    friend inline ExType cosh(const ExType& x) {
      return ExType::cosh(x);
    }
 
 
    static ExType atanh(const ExType& x) {
      return ExType::unary(OP_ATANH, x);
    }
    friend inline ExType atanh(const ExType& x) {
      return ExType::atanh(x);
    }
 
 
    static ExType asinh(const ExType& x) {
      return ExType::unary(OP_ASINH, x);
    }
    friend inline ExType asinh(const ExType& x) {
      return ExType::asinh(x);
    }
 
 
    static ExType acosh(const ExType& x) {
      return ExType::unary(OP_ACOSH, x);
    }
    friend inline ExType acosh(const ExType& x) {
      return ExType::acosh(x);
    }
 
 
    static ExType exp(const ExType& x) {
      return ExType::unary(OP_EXP, x);
    }
    friend inline ExType exp(const ExType& x) {
      return ExType::exp(x);
    }
 
 
    static ExType log(const ExType& x) {
      return ExType::unary(OP_LOG, x);
    }
    friend inline ExType log(const ExType& x) {
      return ExType::log(x);
    }
 
 
    static ExType log10(const ExType& x) {
      return log(x)*(1/std::log(10.));
    }
    friend inline ExType log10(const ExType& x) {
      return ExType::log10(x);
    }
 
 
    static ExType log1p(const ExType& x) {
      return ExType::unary(OP_LOG1P, x);
    }
    friend inline ExType log1p(const ExType& x) {
      return ExType::log1p(x);
    }
 
 
    static ExType expm1(const ExType& x) {
      return ExType::unary(OP_EXPM1, x);
    }
    friend inline ExType expm1(const ExType& x) {
      return ExType::expm1(x);
    }
 
 
    static ExType floor(const ExType& x) {
      return ExType::unary(OP_FLOOR, x);
    }
    friend inline ExType floor(const ExType& x) {
      return ExType::floor(x);
    }
 
 
    static ExType ceil(const ExType& x) {
      return ExType::unary(OP_CEIL, x);
    }
    friend inline ExType ceil(const ExType& x) {
      return ExType::ceil(x);
    }
 
 
    static ExType erf(const ExType& x) {
      return ExType::unary(OP_ERF, x);
    }
    friend inline ExType erf(const ExType& x) {
      return ExType::erf(x);
    }
 
 
    static ExType erfinv(const ExType& x) {
      return ExType::unary(OP_ERFINV, x);
    }
    friend inline ExType erfinv(const ExType& x) {
      return ExType::erfinv(x);
    }
 
 
    static ExType sign(const ExType& x) {
      return ExType::unary(OP_SIGN, x);
    }
    friend inline ExType sign(const ExType& x) {
      return ExType::sign(x);
    }
 
 
    static ExType pow(const ExType& x, const ExType& y) {
      return ExType::binary(OP_POW, x, y);
    }
    friend inline ExType pow(const ExType& x, const ExType& y) {
      return ExType::pow(x, y);
    }
 
 
    static ExType mod(const ExType& x, const ExType& y) {
      return ExType::binary(OP_FMOD, x, y);
    }
    friend inline ExType mod(const ExType& x, const ExType& y) {
      return ExType::mod(x, y);
    }
    friend inline ExType fmod(const ExType& x, const ExType& y) {
      return mod(x, y);
    }
 
 
    static ExType remainder(const ExType& x, const ExType& y) {
      return ExType::binary(OP_REMAINDER, x, y);
    }
    friend inline ExType remainder(const ExType& x, const ExType& y) {
      return ExType::remainder(x, y);
    }
 
 
    static ExType atan2(const ExType& y, const ExType& x) {
      return ExType::binary(OP_ATAN2, y, x);
    }
    friend inline ExType atan2(const ExType& y, const ExType& x) {
      return ExType::atan2(y, x);
    }
 
 
    static ExType if_else_zero(const ExType& x, const ExType& y) {
      return ExType::binary(OP_IF_ELSE_ZERO, x, y);
    }
    friend inline ExType if_else_zero(const ExType& x, const ExType& y) {
      return ExType::if_else_zero(x, y);
    }
 
 
    static ExType fmin(const ExType& x, const ExType& y) {
      return ExType::binary(OP_FMIN, x, y);
    }
    friend inline ExType fmin(const ExType& x, const ExType& y) {
      return ExType::fmin(x, y);
    }
 
 
    static ExType fmax(const ExType& x, const ExType& y) {
      return ExType::binary(OP_FMAX, x, y);
    }
    friend inline ExType fmax(const ExType& x, const ExType& y) {
      return ExType::fmax(x, y);
    }
 
 
     friend inline bool is_equal(const ExType& x, const ExType& y, casadi_int depth=0) {
       return ExType::is_equal(x, y, depth);
     }
 
 
     static ExType copysign(const ExType& x, const ExType& y) {
       return ExType::binary(OP_COPYSIGN, x, y);
     }
     friend inline ExType copysign(const ExType& x, const ExType& y) {
       return ExType::copysign(x, y);
     }
 
 
     static ExType constpow(const ExType& x, const ExType& y) {
       return ExType::binary(OP_CONSTPOW, x, y);
     }
     friend inline ExType constpow(const ExType& x, const ExType& y) {
       return ExType::constpow(x, y);
     }
 
 
     static ExType printme(const ExType& x, const ExType& y) {
       return ExType::binary(OP_PRINTME, x, y);
     }
     friend inline ExType printme(const ExType& x, const ExType& y) {
       return ExType::binary(OP_PRINTME, x, y);
     }
 
 
    static ExType hypot(const ExType& x, const ExType& y) {
      return ExType::binary(OP_HYPOT, x, y);
    }
    friend inline ExType hypot(const ExType& x, const ExType& y) {
      return ExType::hypot(x, y);
    }
 
 
};
#endif // SWIG
 
} // namespace casadi
 
#endif // CASADI_GENERIC_EXPRESSION_HPP