
    using namespace boost::numeric::ublas;
    vector&lt;double&gt; v (3);
    for (unsigned i = 0; i &lt; v.size (); ++ i)
        v (i) = i;

    std::cout &lt;&lt; sum (v) &lt;&lt; std::endl;
    std::cout &lt;&lt; norm_1 (v) &lt;&lt; std::endl;
    std::cout &lt;&lt; norm_2 (v) &lt;&lt; std::endl;
    std::cout &lt;&lt; norm_inf (v) &lt;&lt; std::endl;
    std::cout &lt;&lt; index_norm_inf (v) &lt;&lt; std::endl;
}
</pre>
<h3>Binary Reductions</h3>
<h4>Prototypes</h4>
<pre>
<code>template&lt;class E1, 