vpRowVector.h

/*
 * ViSP, open source Visual Servoing Platform software.
 * Copyright (C) 2005 - 2025 by Inria. All rights reserved.
 *
 * This software is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 * See the file LICENSE.txt at the root directory of this source
 * distribution for additional information about the GNU GPL.
 *
 * For using ViSP with software that can not be combined with the GNU
 * GPL, please contact Inria about acquiring a ViSP Professional
 * Edition License.
 *
 * See https://visp.inria.fr for more information.
 *
 * This software was developed at:
 * Inria Rennes - Bretagne Atlantique
 * Campus Universitaire de Beaulieu
 * 35042 Rennes Cedex
 * France
 *
 * If you have questions regarding the use of this file, please contact
 * Inria at visp@inria.fr
 *
 * This file is provided AS IS with NO WARRANTY OF ANY KIND, INCLUDING THE
 * WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE.
 *
 * Description:
 * Operation on row vectors.
 */
 
#ifndef VP_ROW_VECTOR_H
#define VP_ROW_VECTOR_H
 
#include <vector>
 
#include <visp3/core/vpConfig.h>
#include <visp3/core/vpArray2D.h>
#include <visp3/core/vpColVector.h>
#include <visp3/core/vpMath.h>
#include <visp3/core/vpMatrix.h>
 
BEGIN_VISP_NAMESPACE
 
class vpMatrix;
class vpColVector;


class VISP_EXPORT vpRowVector : public vpArray2D<double>
{
public:
  vpRowVector() : vpArray2D<double>() { }
  VP_EXPLICIT vpRowVector(unsigned int n) : vpArray2D<double>(1, n) { }
  vpRowVector(unsigned int n, double val) : vpArray2D<double>(1, n, val) { }
  vpRowVector(const vpRowVector &v) : vpArray2D<double>(v) { }
  vpRowVector(const vpRowVector &v, unsigned int c, unsigned int ncols);
  VP_EXPLICIT vpRowVector(const vpMatrix &M);
  vpRowVector(const vpMatrix &M, unsigned int i);
  VP_EXPLICIT vpRowVector(const std::vector<double> &v);
  VP_EXPLICIT vpRowVector(const std::vector<float> &v);
#if (VISP_CXX_STANDARD >= VISP_CXX_STANDARD_11)
  vpRowVector(vpRowVector &&v);
  VP_EXPLICIT vpRowVector(const std::initializer_list<double> &list) : vpArray2D<double>(list) { }
#endif
 
  static vpRowVector view(double *raw_data, unsigned int ncols);

  void clear()
  {
    if (!isMemoryOwner) {
      throw vpException(vpException::fatalError, "Cannot clear a vector view");
    }
    if (data != nullptr) {
      free(data);
      data = nullptr;
    }
 
    if (rowPtrs != nullptr) {
      free(rowPtrs);
      rowPtrs = nullptr;
    }
    rowNum = 0;
    colNum = 0;
    dsize = 0;
  }
 
  std::ostream &cppPrint(std::ostream &os, const std::string &matrixName = "A", bool octet = false) const;
  std::ostream &csvPrint(std::ostream &os) const;

  inline void deg2rad()
  {
    double d2r = M_PI / 180.0;
 
    (*this) *= d2r;
  }

  vpRowVector extract(unsigned int c, unsigned int rowsize) const
  {
    if ((c >= colNum) || ((c + rowsize) > colNum)) {
      throw(vpException(vpException::fatalError,
                        "Cannot extract a (1x%d) row vector from a (1x%d) "
                        "row vector starting at index %d",
                        rowsize, colNum, c));
    }
 
    return vpRowVector(*this, c, rowsize);
  }
 
  double frobeniusNorm() const;

  vpRowVector hadamard(const vpRowVector &v) const;
 
  void init(const vpRowVector &v, unsigned int c, unsigned int ncols);
  void insert(unsigned int i, const vpRowVector &v);
 
  std::ostream &maplePrint(std::ostream &os) const;
  std::ostream &matlabPrint(std::ostream &os) const;
 
  vpRowVector &normalize();
  vpRowVector &normalize(vpRowVector &x) const;

  inline double &operator[](unsigned int n) { return *(data + n); }
  inline const double &operator[](unsigned int n) const { return *(data + n); }

  vpRowVector &operator=(const vpRowVector &v);
  vpRowVector &operator=(const vpMatrix &M);
  vpRowVector &operator=(const std::vector<double> &v);
  vpRowVector &operator=(const std::vector<float> &v);
  vpRowVector &operator=(double x);
#if (VISP_CXX_STANDARD >= VISP_CXX_STANDARD_11)
  vpRowVector &operator=(vpRowVector &&v);
  vpRowVector &operator=(const std::initializer_list<double> &list);
#endif
  bool operator==(const vpRowVector &v) const;
  bool operator!=(const vpRowVector &v) const;
 
  double operator*(const vpColVector &x) const;
  vpRowVector operator*(const vpMatrix &M) const;
  vpRowVector operator*(double x) const;
  vpRowVector &operator*=(double x);
 
  vpRowVector operator/(double x) const;
  vpRowVector &operator/=(double x);
 
  vpRowVector operator+(const vpRowVector &v) const;
  vpRowVector &operator+=(vpRowVector v);
 
  vpRowVector operator-(const vpRowVector &v) const;
  vpRowVector &operator-=(vpRowVector v);
  vpRowVector operator-() const;
 
  vpRowVector &operator<<(const vpRowVector &v);
  vpRowVector &operator<<(double val);
  vpRowVector &operator,(double val);
 
  int print(std::ostream &s, unsigned int length, char const *intro = nullptr) const;
  inline void rad2deg()
  {
    double r2d = 180.0 / M_PI;
 
    (*this) *= r2d;
  }
 
  void reshape(vpMatrix &M, const unsigned int &nrows, const unsigned int &ncols);
  vpMatrix reshape(unsigned int nrows, unsigned int ncols);

  inline void resize(unsigned int i, bool flagNullify = true) { vpArray2D<double>::resize(1, i, flagNullify); }

  void resize(unsigned int nrows, unsigned int ncols, bool flagNullify)
  {
    if (nrows != 1) {
      throw(vpException(vpException::fatalError,
                        "Cannot resize a row vector to a (%dx%d) dimension "
                        "vector that has more than one row",
                        nrows, ncols));
    }
    vpArray2D<double>::resize(nrows, ncols, flagNullify);
  }
 
  void stack(double d);
  void stack(const vpRowVector &v);
 
  double sum() const;
  double sumSquare() const;
  vpColVector t() const;
  std::vector<double> toStdVector() const;
  vpColVector transpose() const;
  void transpose(vpColVector &v) const;
 
  static double mean(const vpRowVector &v);
  static double median(const vpRowVector &v);
  static vpRowVector stack(const vpRowVector &A, const vpRowVector &B);
  static void stack(const vpRowVector &A, const vpRowVector &B, vpRowVector &C);
  static double stdev(const vpRowVector &v, bool useBesselCorrection = false);
 
#if defined(VISP_BUILD_DEPRECATED_FUNCTIONS)

  VP_DEPRECATED void init() { }
  VP_DEPRECATED void stackMatrices(const vpRowVector &r) { stack(r); }
  VP_DEPRECATED static vpRowVector stackMatrices(const vpRowVector &A, const vpRowVector &B) { return stack(A, B); }
  VP_DEPRECATED static void stackMatrices(const vpRowVector &A, const vpRowVector &B, vpRowVector &C)
  {
    stack(A, B, C);
  }

  VP_DEPRECATED void setIdentity(const double &val = 1.0);
  VP_DEPRECATED double euclideanNorm() const;
#endif
};
 
VISP_EXPORT vpRowVector operator*(const double &x, const vpRowVector &v);
 
END_VISP_NAMESPACE
#endif