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