JVectorND.hh

Go to the documentation of this file.

#ifndef __JMATH__JVECTORND__
#define __JMATH__JVECTORND__

#include <vector>

namespace JMATH {}
namespace JPP { using namespace JMATH; }

namespace JMATH {

  /**
   * Nx1 matrix.
   */
  struct JVectorND :
    public std::vector<double>
  {
    /**
     * Default constructor.
     */
    JVectorND() :
      std::vector<double>()
    {}


    /**
     * Constructor.
     *
     * \param  size          size
     */
    JVectorND(size_t size) :
      std::vector<double>(size, 0.0)
    {}
 

    /**
     * Get dot product.
     *
     * \param  vector        vector
     * \return               dot product
     */
    double getDot(const JVectorND& vector) const
    {
      double dot = 0.0;

      const double* p = this-> data();
      const double* q = vector.data();

      for (size_t i = this->size(); i != 0; --i, ++p, ++q) {
        dot += (*p) * (*q);
      }

      return dot;
    }
  };
}

#endif