std::inner_product

template< class InputIt1, class InputIt2, class T >
T inner_product( InputIt1 first1, InputIt1 last1,
                 InputIt2 first2, T init );

template< class InputIt1, class InputIt2, class T,
          class BinaryOp1, class BinaryOp2 >
T inner_product( InputIt1 first1, InputIt1 last1,
                 InputIt2 first2, T init,
                 BinaryOp1 op1, BinaryOp2 op2 );

Calcule le produit scalaire (c'est-à-dire la somme des produits) ou effectue une opération map/reduce ordonnée sur la plage [first1, last1) et la plage de std::distance(first1, last1) éléments commençant à first2.

Étant donné last2 comme le std::distance(first1, last1)ième prochain itérateur de first2, si l'une des conditions suivantes est satisfaite, le comportement est indéfini :

• T n'est pas CopyConstructible.
• T n'est pas CopyAssignable.
• op1 ou op2 modifie tout élément de [first1, last1) ou [first2, last2).
• op1 ou op2 invalide tout itérateur ou sous-plage dans [first1, last1] ou [first2, last2].

Paramètres

Valeur de retour

acc après toutes les modifications.

Implémentation possible

template<class InputIt1, class InputIt2, class T>
constexpr // since C++20
T inner_product(InputIt1 first1, InputIt1 last1, InputIt2 first2, T init)
{
    while (first1 != last1)
    {
        init = std::move(init) + (*first1) * (*first2); // std::move since C++20
        ++first1;
        ++first2;
    }
    
    return init;
}

template<class InputIt1, class InputIt2, class T,
         class BinaryOp1, class BinaryOp2>
constexpr // since C++20
T inner_product(InputIt1 first1, InputIt1 last1, InputIt2 first2, T init,
                BinaryOp1 op1, BinaryOp2 op2)
{
    while (first1 != last1)
    {
        init = op1(std::move(init), op2(*first1, *first2)); // std::move since C++20
        ++first1;
        ++first2;
    }
    
    return init;
}

Notes

La version parallélisable de cet algorithme, std::transform_reduce, nécessite op1 et op2 d'être commutatifs et associatifs, mais std::inner_product n'impose pas une telle exigence, et effectue toujours les opérations dans l'ordre donné.

Exemple

#include <functional>
#include <iostream>
#include <numeric>
#include <vector>

int main()
{
    std::vector<int> a{0, 1, 2, 3, 4};
    std::vector<int> b{5, 4, 2, 3, 1};
    
    int r1 = std::inner_product(a.begin(), a.end(), b.begin(), 0);
    std::cout << "Inner product of a and b: " << r1 << '\n';
    
    int r2 = std::inner_product(a.begin(), a.end(), b.begin(), 0,
                                std::plus<>(), std::equal_to<>());
    std::cout << "Number of pairwise matches between a and b: " <<  r2 << '\n';
}

Inner product of a and b: 21
Number of pairwise matches between a and b: 2

Rapports de défauts

Les rapports de défauts suivants modifiant le comportement ont été appliqués rétroactivement aux normes C++ précédemment publiées.

Voir aussi