product.h

#ifndef UTIL_PRODUCT_H
#define UTIL_PRODUCT_H
 
/*
* Util Package - C++ Utilities for Scientific Computation
*
* Copyright 2010 - 2017, The Regents of the University of Minnesota
* Distributed under the terms of the GNU General Public License.
*/
 
#include <util/space/Vector.h>
#include <util/space/Tensor.h>
#include <complex>
using std::complex;
 
namespace Util
{
 
   inline float product(float a, float b)
   {  return a*b; }
 
   inline double product(double a, double b)
   {  return a*b; }
 
   inline double product(const Vector& a, const Vector& b)
   {  return a.dot(b); }
 
   inline double product(const Tensor& a, const Tensor& b)
   {
      double sum = 0.0;
      int i, j;
      for (i = 0; i < Dimension; ++i) {
         for (j = 0; j < Dimension; ++j) {
            sum += a(i, j)*b(i, j);
         }
      }
      return sum;
   }
 
   inline complex<float> product(complex<float> a, complex<float> b)
   {  return conj(a)*b; }
 
   inline complex<double> product(complex<double> a, complex<double> b)
   {  return conj(a)*b; }
 
}
#endif