// fig06_12.cpp
// Compute the product of an array's elements using accumulate.
#include <array>
#include <iostream>
#include <numeric>

template <typename T>
T multiply(T x, T y) {
   return x * y;
}


auto multiply_auto(auto x, auto y) {
   return x * y;
}

int multiply_int(int x, int y) {
   return x * y;
}


double f(double& x){
	return x*x;
}

template <typename T, typename F, long unsigned int N>
T integrate (std::array<T, N>& x, F f){
	T tmp{T(0)};
	for (auto s: x){
            tmp+= f(s);
	}
	return tmp;
}


auto integrate_auto (auto x, auto f){
        auto tmp{x[0]};
	
        for (auto s: x){
            tmp+= f(s);
        }
        tmp-=x[0];
        return tmp;
}


typedef double (*FUNC)(double&);


template <typename T, long unsigned int N>
T integrate_legacy (std::array<T, N> x, FUNC f){
        T tmp{T(0)};
        for (auto s: x){
            tmp+= f(s);
         }
         return tmp;
}


int main() {
   constexpr std::array integers{1, 2, 3, 4, 5};
   std::array floats {1.0, 2.0, 3.0, 4.0, 5.0};
   auto multiply_lambda = [](const auto& x, const auto& y) {return x * y; };

   auto f_lambda = [](auto& x) {return x*x;};

   auto g_lambda = [&integers](auto& x) {
	   auto tmp{0};
	   for (int i{0};auto s: x){
               tmp += s*integers[i];
	   }
	   return tmp;
   };

   std::cout << "Product of integers: " <<  std::accumulate( std::begin(integers), std::end(integers), 1, multiply<int>) << "\n";

  //std::cout << "Product of integers: " <<  std::accumulate( std::begin(integers), std::end(integers), 1, multiply_auto) << "\n";

   std::cout << "Product of integers: " <<  std::accumulate( integers.cbegin(), integers.cend(), 1, multiply<int>) << "\n";

   std::cout << "Product of integers with a lambda: " << std::accumulate(std::begin(integers), std::end(integers), 1, multiply_lambda) << "\n";

   std::cout << "Product of integers with a lambda: " << std::accumulate(std::begin(floats), std::end(floats), 1, multiply_int) << "\n";

   std::cout << "Product of integers with a lambda: " << std::accumulate(std::begin(floats), std::end(floats), 1, multiply_lambda) << "\n";

   std::cout << "Integration2  template  " << integrate_legacy(floats, f_lambda) << "\n";

   std::cout << "Integration  full  template  " << integrate(floats, f) << "\n";

   std::cout << "Integration  auto          " << integrate_auto(floats, f_lambda) << "\n";

   std::cout << "Lambda test   " << g_lambda(floats) << "\n";
}


/**************************************************************************
 * (C) Copyright 1992-2022 by Deitel & Associates, Inc. and               *
 * Pearson Education, Inc. All Rights Reserved.                           *
 *                                                                        *
 * DISCLAIMER: The authors and publisher of this book have used their     *
 * best efforts in preparing the book. These efforts include the          *
 * development, research, and testing of the theories and programs        *
 * to determine their effectiveness. The authors and publisher make       *
 * no warranty of any kind, expressed or implied, with regard to these    *
 * programs or to the documentation contained in these books. The authors *
 * and publisher shall not be liable in any event for incidental or       *
 * consequential damages in connection with, or arising out of, the       *
 * furnishing, performance, or use of these programs.                     *
 **************************************************************************/