{

auto vecf1 = pyvectorize(f1);

shape_type shape = { 3, 2 };

pyarray<double> a(shape, 1.5);

pyarray<double> b(shape, 2.3);

pyarray<double> c = vecf1(a, b);

EXPECT_EQ(a(0, 0) + b(0, 0), c(0, 0));

{

auto vecf1 = pyvectorize([](double a, double b) { return a + b; });

shape_type shape = { 3, 2 };

pyarray<double> a(shape, 1.5);

pyarray<double> b(shape, 2.3);

pyarray<double> c = vecf1(a, b);

EXPECT_EQ(a(0, 0) + b(0, 0), c(0, 0));

{

using complex_t = std::complex<double>;

shape_type shape = { 3, 2 };

pyarray<complex_t> a(shape, complex_t(1.2, 2.5));

auto f = pyvectorize([](complex_t x) { return std::abs(x); });

auto res = f(a);

double exp = std::abs(a(1, 1));

EXPECT_EQ(exp, res(1, 1));