_drovector-_dcovector.hpp

Go to the documentation of this file.

//=============================================================================
/*! _drovector*_dcovector operator */
inline double operator*(const _drovector& rovec, const _dcovector& covec)
{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(rovec.l!=covec.l){
    ERROR_REPORT;
    std::cerr << "These two vectors can not make a product." << std::endl
              << "Your input was (" << rovec.l << ") * (" << covec.l << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  CPPL_INT inc =1;
  double val =ddot_( &rovec.l, rovec.array, &inc, covec.array, &inc );
  
  rovec.destroy();
  covec.destroy();
  return val;
}