ra-ra/examples/outer.C

// Daniel Llorens - 2015
// Adapted from blitz++/examples/outer.cpp

#include "ra/ra-operators.H"

using std::cout; using std::endl; using std::flush;

int main()
{
    ra::Owned<float,1> x { 1, 2, 3, 4 }, y { 1, 0, 0, 1 };
    ra::Owned<float,2> A({4,4}, 99.);

    x = { 1, 2, 3, 4 };
    y = { 1, 0, 0, 1 };

    ra::TensorIndex<0> i;
    ra::TensorIndex<1> j;

    A = x(i) * y(j);

    cout << A << endl;

// [ra] an alternative

    cout << from(std::multiplies<float>(), x, y) << endl;

    return 0;
}