#include <iostream>

#include "xtensor/xarray.hpp"
#include "xtensor/xio.hpp"
#include "xtensor/xview.hpp"

xt::xarray<double> arr1
  {{1.0, 2.0, 3.0},
   {2.0, 5.0, 7.0},
   {2.0, 5.0, 7.0}};

xt::xarray<double> arr2
  {5.0, 6.0, 7.0};

xt::view(arr1, 1) + arr2

arr1

#include <iostream>
#include "xtensor/xarray.hpp"
#include "xtensor/xio.hpp"

xt::xarray<int> arr
  {1, 2, 3, 4, 5, 6, 7, 8, 9};

arr.reshape({3, 3});

arr

#include <iostream>
#include "xtensor/xarray.hpp"
#include "xtensor/xmath.hpp"
#include "xtensor/xio.hpp"

xt::xarray<double> arr3
  {1.0, 2.0, 3.0};

xt::xarray<unsigned int> arr4
  {4, 5, 6, 7};

arr4.reshape({4, 1});

xt::pow(arr3, arr4)

#include <iostream>
#include "xtensor/xrandom.hpp"

xt::xarray<double> arr5 = xt::random::randn<double>({4, 3});
arr5

xt::random::randn<double>({5, 3, 8, 10}) > 0.2

#include "xtensor/xbuilder.hpp"

xt::xarray<double> ar = xt::linspace<double>(0.0, 10.0, 12);
ar.reshape({4, 3});
ar

xt::xarray<double> fones = xt::ones<float>({2, 2});
fones

xt::arange<int>(1569325055)

#include <vector>
#include "xtensor/xbroadcast.hpp"

xt::broadcast(xt::linspace<double>(0.0, 10.0, 4),
              std::vector<std::size_t>({3, 4}))

#include <algorithm>

xt::xarray<double> frand = xt::random::randn<double>({2, 2});

// begin() and end() provide and iterator pair iterating over the xexpression in a row-major fashion
std::cout << std::accumulate(frand.begin(), frand.end(), 0.0);

frand

// begin(shape) and end(shape) provide and iterator pair iterating
// over the xexpression broadcasted to the prescrived shape in a row-major fashion
std::vector<std::size_t> shape = {3, 2, 2};
std::cout << std::accumulate(frand.begin(shape), frand.end(shape), 0.0);

#include "xtensor-blas/xlinalg.hpp"

xt::xtensor<double, 2> m = {{1.5, 0.5}, {0.7, 1.0}};

std::cout << "Matrix rank: " << std::endl << xt::linalg::matrix_rank(m);

std::cout << "Matrix inverse: " << std::endl << xt::linalg::inv(m);

std::cout << "Eigen values: " << std::endl << xt::linalg::eigvals(m);

xt::xarray<double> arg1 = xt::arange<double>(9);
xt::xarray<double> arg2 = xt::arange<double>(18);
arg1.reshape({3, 3});
arg2.reshape({2, 3, 3});
std::cout << xt::linalg::dot(arg1, arg2) << std::endl;