fixed_row_vector.h

/* Copyright (C) 2003-2005  Martin Förg
 
 * This library 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 2.1 of the License, or (at your option) any later version.
 * 
 * This library 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 this library; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
 *
 * Contact:
 *   martin.o.foerg@googlemail.com
 *
 */
 
#ifndef fixed_row_vector_h
#define fixed_row_vector_h
 
#include "fixed_general_matrix.h"
#include <stdexcept>
 
namespace fmatvec {
 
  template <int N, class AT> class RowVector<Fixed<N>,AT> : public Matrix<General,Fixed<1>,Fixed<N>,AT> {
 
    using Matrix<General,Fixed<1>,Fixed<N>,AT>::ele;
 
    public:
      static constexpr bool isVector {true};
 
      using iterator = AT *;
      using const_iterator = const AT *;
 
      using value_type = AT;
 
    
    protected:
 
    template<class Col> inline RowVector<Fixed<N>,AT>& copy(const RowVector<Col,AT> &x);
 
 
    public:
 
//      template<class Ini=All<AT>>
//        RowVector(Ini ini=All<AT>()) : Matrix<General,Fixed<1>,Fixed<N>,AT>(ini) { }
//      template<class Ini=All<AT>>
//        RowVector(int n, Ini ini=All<AT>()) : Matrix<General,Fixed<1>,Fixed<N>,AT>(ini) { }
 
      explicit RowVector(Noinit ini) : Matrix<General,Fixed<1>,Fixed<N>,AT>(ini) { }
      explicit RowVector(Init ini=INIT, const AT &a=AT()) : Matrix<General,Fixed<1>,Fixed<N>,AT>(ini,a) { }
      explicit RowVector(int n, Noinit ini) : Matrix<General,Fixed<1>,Fixed<N>,AT>(ini) { FMATVEC_ASSERT(n==N, AT); }
      explicit RowVector(int n, Init ini=INIT, const AT &a=AT()) : Matrix<General,Fixed<1>,Fixed<N>,AT>(ini,a) { FMATVEC_ASSERT(n==N, AT); }
 
      RowVector(const RowVector<Fixed<N>,AT> &A) = default;
 
      template<class Row>
      RowVector(const RowVector<Row,AT> &A) : Matrix<General,Fixed<1>,Fixed<N>,AT>(A) {
      }
 
      template<class Type, class Row, class Col>
      explicit RowVector(const Matrix<Type,Row,Col,AT> &A) : Matrix<General,Fixed<1>,Fixed<N>,AT>(A) {
      }
 
      RowVector(const char *str) : Matrix<General,Fixed<1>,Fixed<N>,AT>(str) {
      }
 
      inline RowVector<Fixed<N>,AT>& operator=(const RowVector<Fixed<N>,AT> &x) = default;
 
      template <class Row>
      inline RowVector<Fixed<N>,AT>& operator=(const RowVector<Row,AT> &x) {
        FMATVEC_ASSERT(x.size() == N, AT);
        return copy(x);
      }
 
      template <class Row>
      inline RowVector<Fixed<N>,AT>& operator<<=(const RowVector<Row,AT> &x) { return operator=(x); }
 
      AT& operator()(int i) {
    FMATVEC_ASSERT(i>=0, AT);
    FMATVEC_ASSERT(i<N, AT);
    return e(i);
      }
 
      const AT& operator()(int i) const {
    FMATVEC_ASSERT(i>=0, AT);
    FMATVEC_ASSERT(i<N, AT);
    return e(i);
      }
 
      iterator begin() { return &ele[0][0]; }
      iterator end() { return &ele[0][N]; }
      const_iterator begin() const { return &ele[0][0]; }
      const_iterator end() const { return &ele[0][N]; }
      const_iterator cbegin() const noexcept { return &ele[0][0]; }
      const_iterator cend() const noexcept { return &ele[0][N]; }
 
      AT& e(int i) {
    return ele[0][i];
      }
 
      const AT& e(int i) const {
    return ele[0][i];
      }
 
      inline RowVector<Fixed<N>,AT>& init(const AT &val=AT()); 
      inline RowVector<Fixed<N>,AT>& init(Init, const AT &a=AT()) { return init(a); }
      inline RowVector<Fixed<N>,AT>& init(Noinit, const AT &a=AT()) { return *this; }
 
      constexpr int size() const {return N;}
 
      void resize(int n) {
        if(n!=N)
          throw std::runtime_error("A fixed row vector cannot be resized.");
      }
 
      int inc() const {return 1;}
 
      using Matrix<General,Fixed<1>,Fixed<N>,AT>::operator();
 
      explicit inline operator std::vector<AT>() const;
 
      explicit inline RowVector(const std::vector<AT> &v);
 
//      /*! \brief Cast to AT.
//       *
//       * \return The AT representation of the vector
//       * */
//      explicit operator AT() const {
//        FMATVEC_ASSERT(N==1, AT);
//        return ele[0];
//      }
//
//      /*! \brief AT Constructor.
//       * Constructs and initializes a vector with a AT object.
//       * \param x The AT the vector will be initialized with.
//       * */
//      explicit RowVector(const AT &x) : Matrix<General,Fixed<1>,Fixed<N>,AT>(x) { }
 
      inline const Vector<Fixed<N>,AT> T() const;
  };
 
  template <int N, class AT>
    inline RowVector<Fixed<N>,AT>& RowVector<Fixed<N>,AT>::init(const AT &val) {
      for(int i=0; i<N; i++) 
    e(i) = val; 
      return *this;
    }
 
  template <int N, class AT>
    inline const Vector<Fixed<N>,AT> RowVector<Fixed<N>,AT>::T() const {
      Vector<Fixed<N>,AT> x(NONINIT);
      for(int i=0; i<N; i++)
        x.e(i) = e(i);
      return x;
    }
 
  template <int N, class AT>
   inline RowVector<Fixed<N>,AT>::operator std::vector<AT>() const {
      std::vector<AT> ret(size());
      for(int i=0; i<size(); ++i)
        ret[i] = e(i);
      return ret;
    }
 
  template <int N, class AT>
    inline RowVector<Fixed<N>,AT>::RowVector(const std::vector<AT> &v) : Matrix<General,Fixed<1>,Fixed<N>,AT>(NONINIT) {
      for(int i=0; i<size(); ++i)
        e(i) = v[i];
    }
 
 
  template <int N, class AT> template <class Col>
    inline RowVector<Fixed<N>,AT>& RowVector<Fixed<N>,AT>::copy(const RowVector<Col,AT> &x) {
      for(int i=0; i<N; i++)
        e(i) = x.e(i);
      return *this;
    }
 
 
}
#endif