codac2_Matrix_addons_Base.h

Go to the documentation of this file.

template<int R=RowsAtCompileTime,int C=ColsAtCompileTime>
inline Scalar& operator()(Index i, Index j)
{
  return const_cast<Scalar&>(static_cast<const Matrix<Scalar,R,C>&>(*this).operator()(i,j));
}

template<int R=RowsAtCompileTime,int C=ColsAtCompileTime>
inline const Scalar& operator()(Index i, Index j) const
{
  assert_release(i >= 0 && i < this->rows() && j >= 0 && j < this->cols());
 
  if constexpr(!IsVectorOrRow<R,C>)
    return this->PlainObjectBase<Matrix<Scalar,R,C>>::operator()(i,j);
 
  else
  {
    if constexpr(R == 1)
      return (*this)[j];
    return (*this)[i];
  }
}

inline auto& init(const Scalar& x)
{
  for(Index i = 0 ; i < this->size() ; i++)
    *(this->data()+i) = x;
  return *this;
}

template<typename U_,int R_,int C_>
inline bool operator==(const Matrix<U_,R_,C_>& x) const
{
  // This operator should be related to Eigen::MatrixBase,
  // for allowing a comparison between two matrix expressions.
  // However, it is not possible to overwrite Eigen::MatrixBase::operator==.
  // Therefore, only comparisons between Eigen::Matrix types is possible.
  // This must be corrected...
 
  if(this->size() != x.size())
    return false;
 
  if constexpr(std::is_same_v<Scalar,codac2::Interval> && std::is_same_v<U_,codac2::Interval>)
  {
    if(this->is_empty() || x.is_empty())
      return this->is_empty() && x.is_empty();
  }
 
  if(this->rows() != x.rows() || this->cols() != x.cols())
    return false;
 
  for(Index i = 0 ; i < this->size() ; i++)
    if(*(this->data()+i) != *(x.data()+i))
      return false;
 
  return true;
}