CPPScaLapack: _drovector-_dcovector.hpp Source File

00001 //=============================================================================
00002 /*! _drovector*_dcovector operator */
00003 inline double operator*(const _drovector& rovec, const _dcovector& covec)
00004 {
00005 #ifdef  CPPSL_DEBUG
00006   if(rovec.L!=covec.L){
00007     std::cerr << "[ERROR] operator*(const _drovector&, const _dcovector&)"
00008               << std::endl
00009               << "These two vectors can not make a product." << std::endl
00010               << "Your input was (" << rovec.L << ") * (" << covec.L << ")."
00011               << std::endl;
00012     exit(1);
00013   }
00014 #endif//CPPSL_DEBUG
00015   
00016   double val;
00017   pddot_( rovec.L, val,
00018           rovec.Array, 1, 1, rovec.Desc, 1, covec.Array, 1, 1, covec.Desc, 1 );
00019   
00020   rovec.destroy();
00021   covec.destroy();
00022   return val;
00023 }