Add CMake examples

Author: tdegeus

.azure-pipelines/unix-build.yml

@@ -11,15 +11,16 @@ steps:
       mkdir build
       cd build
       cmake -DCMAKE_INSTALL_PREFIX=$CONDA_PREFIX -DPYTHON_EXECUTABLE=`which python` -DDOWNLOAD_GTEST=ON $(Build.SourcesDirectory)
+      make install
     displayName: Configure xtensor-python
     workingDirectory: $(Build.BinariesDirectory)
-      
+
   - script: |
       source activate xtensor-python
       make -j2 test_xtensor_python
     displayName: Build xtensor-python
     workingDirectory: $(Build.BinariesDirectory)/build
-          
+
   - script: |
       source activate xtensor-python
       cd test
@@ -32,3 +33,19 @@ steps:
       py.test -s
     displayName: Test xtensor-python (Python)
     workingDirectory: $(Build.SourcesDirectory)
+
+  - script: |
+      source activate xtensor-python
+      cmake . -DPYTHON_EXECUTABLE=`which python`
+      cmake --build .
+      python example.py
+    displayName: Example - readme 1
+    workingDirectory: $(Build.SourcesDirectory)/docs/source/examples/readme_example_1
+
+  - script: |
+      source activate xtensor-python
+      cmake . -DPYTHON_EXECUTABLE=`which python`
+      cmake --build .
+      python example.py
+    displayName: Example - SFINAE
+    workingDirectory: $(Build.SourcesDirectory)/docs/source/examples/sfinae

README.md

@@ -77,7 +77,7 @@ import xtensor_python_test as xt
 
 v = np.arange(15).reshape(3, 5)
 s = xt.sum_of_sines(v)
-s
+print(s)
 ```
 
 **Outputs**
@@ -86,6 +86,14 @@ s
 1.2853996391883833
 ```
 
+**Working example**
+
+Get the working example here:
+
+*   [`CMakeLists.txt`](docs/source/examples/readme_example_1/CMakeLists.txt)
+*   [`main.cpp`](docs/source/examples/readme_example_1/main.cpp)
+*   [`example.py`](docs/source/examples/readme_example_1/example.py)
+
 ### Example 2: Create a universal function from a C++ scalar function
 
 **C++ code**
@@ -122,7 +130,7 @@ import xtensor_python_test as xt
 x = np.arange(15).reshape(3, 5)
 y = [1, 2, 3, 4, 5]
 z = xt.vectorized_func(x, y)
-z
+print(z)
 ```
 
 **Outputs**

docs/source/examples.rst

@@ -0,0 +1,100 @@
+
+****************
+(CMake) Examples
+****************
+
+Basic example (from readme)
+===========================
+
+Consider the following C++ code:
+
+:download:`main.cpp <examples/readme_example_1/main.cpp>`
+
+.. literalinclude:: examples/readme_example_1/main.cpp
+   :language: cpp
+
+There are several options to build the module,
+whereby we will CMake here with the following ``CMakeLists.txt``:
+
+:download:`CMakeLists.txt <examples/readme_example_1/CMakeLists.txt>`
+
+.. literalinclude:: examples/readme_example_1/CMakeLists.txt
+   :language: cmake
+
+Then we can test the module:
+
+:download:`example.py <examples/readme_example_1/example.py>`
+
+.. literalinclude:: examples/readme_example_1/example.py
+   :language: cmake
+
+.. note::
+
+    Since we did not install the module,
+    we should compile and run the example from the same folder.
+    To install, please consult
+    `this pybind11 / CMake example <https://github.com/pybind/cmake_example>`_.
+
+
+Type restriction with SFINAE
+============================
+
+.. seealso::
+
+    `Medium post by Johan Mabille <https://medium.com/@johan.mabille/designing-language-bindings-with-xtensor-f32aa0f20db>`__
+    This example covers "Option 4".
+
+In this example we will design a module with a function that accepts an ``xt::xtensor`` as argument,
+but in such a way that an ``xt::pyxtensor`` can be accepted in the Python module.
+This is done by having a templated function
+
+.. code-block:: cpp
+
+    template <class T>
+    void times_dimension(T& t);
+
+As this might be a bit too permissive for your liking, we will show you how to limit the
+scope to *xtensor* types, and allow other overloads using the principle of SFINAE
+(Substitution Failure Is Not An Error).
+In particular:
+
+:download:`mymodule.hpp <examples/sfinae/mymodule.hpp>`
+
+.. literalinclude:: examples/sfinae/mymodule.hpp
+   :language: cpp
+
+Consequently from C++, the interaction with the module's function is trivial
+
+:download:`main.cpp <examples/sfinae/main.cpp>`
+
+.. literalinclude:: examples/sfinae/main.cpp
+   :language: cpp
+
+For the Python module we just have to specify the template to be
+``xt::pyarray`` or ``xt::pytensor``. E.g.
+
+:download:`src/python.cpp <examples/sfinae/python.cpp>`
+
+.. literalinclude:: examples/sfinae/python.cpp
+   :language: cpp
+
+We will again use CMake to compile, with the following ``CMakeLists.txt``:
+
+:download:`CMakeLists.txt <examples/sfinae/CMakeLists.txt>`
+
+.. literalinclude:: examples/sfinae/CMakeLists.txt
+   :language: cmake
+
+Then we can test the module:
+
+:download:`example.py <examples/readme_example_1/example.py>`
+
+.. literalinclude:: examples/readme_example_1/example.py
+   :language: cmake
+
+.. note::
+
+    Since we did not install the module,
+    we should compile and run the example from the same folder.
+    To install, please consult
+    `this pybind11 / CMake example <https://github.com/pybind/cmake_example>`_.

docs/source/examples/readme_example_1/CMakeLists.txt

@@ -0,0 +1,13 @@
+cmake_minimum_required(VERSION 3.1..3.19)
+
+project(mymodule)
+
+find_package(pybind11 CONFIG REQUIRED)
+find_package(xtensor REQUIRED)
+find_package(xtensor-python REQUIRED)
+find_package(Python REQUIRED COMPONENTS NumPy)
+
+pybind11_add_module(mymodule main.cpp)
+target_link_libraries(mymodule PUBLIC pybind11::module xtensor-python Python::NumPy)
+
+target_compile_definitions(mymodule PRIVATE VERSION_INFO=0.1.0)

docs/source/examples/readme_example_1/example.py

@@ -0,0 +1,6 @@
+import mymodule
+import numpy as np
+
+a = np.array([1, 2, 3])
+assert np.isclose(np.sum(np.sin(a)), mymodule.sum_of_sines(a))
+

docs/source/examples/readme_example_1/main.cpp

@@ -0,0 +1,18 @@
+#include <numeric>
+#include <xtensor.hpp>
+#include <pybind11/pybind11.h>
+#define FORCE_IMPORT_ARRAY
+#include <xtensor-python/pyarray.hpp>
+
+double sum_of_sines(xt::pyarray<double>& m)
+{
+    auto sines = xt::sin(m);  // sines does not actually hold values.
+    return std::accumulate(sines.begin(), sines.end(), 0.0);
+}
+
+PYBIND11_MODULE(mymodule, m)
+{
+    xt::import_numpy();
+    m.doc() = "Test module for xtensor python bindings";
+    m.def("sum_of_sines", sum_of_sines, "Sum the sines of the input values");
+}

docs/source/examples/sfinae/CMakeLists.txt

@@ -0,0 +1,16 @@
+cmake_minimum_required(VERSION 3.1..3.19)
+
+project(mymodule)
+
+find_package(pybind11 CONFIG REQUIRED)
+find_package(xtensor REQUIRED)
+find_package(xtensor-python REQUIRED)
+find_package(Python REQUIRED COMPONENTS NumPy)
+
+pybind11_add_module(mymodule python.cpp)
+target_link_libraries(mymodule PUBLIC pybind11::module xtensor-python Python::NumPy)
+
+target_compile_definitions(mymodule PRIVATE VERSION_INFO=0.1.0)
+
+add_executable(myexec main.cpp)
+target_link_libraries(myexec PUBLIC xtensor)

docs/source/examples/sfinae/example.py

@@ -0,0 +1,8 @@
+import mymodule
+import numpy as np
+
+a = np.array([[1, 2, 3], [4, 5, 6]], dtype=np.float64)
+b = np.array(a, copy=True)
+mymodule.times_dimension(b) # changing in-place!
+assert np.allclose(2 * a, b)
+

docs/source/examples/sfinae/main.cpp

@@ -0,0 +1,10 @@
+#include "mymodule.hpp"
+#include <xtensor/xio.hpp>
+
+int main()
+{
+    xt::xtensor<size_t, 2> a = xt::arange<size_t>(2 * 3).reshape({2, 3});
+    mymodule::times_dimension(a);
+    std::cout << a << std::endl;
+    return 0;
+}

docs/source/examples/sfinae/mymodule.hpp

@@ -0,0 +1,32 @@
+#include <xtensor/xtensor.hpp>
+
+namespace mymodule {
+
+template <class T>
+struct is_std_vector
+{
+    static const bool value = false;
+};
+
+template <class T>
+struct is_std_vector<std::vector<T> >
+{
+    static const bool value = true;
+};
+
+// any xtensor object
+template <class T, std::enable_if_t<xt::is_xexpression<T>::value, bool> = true>
+void times_dimension(T& t)
+{
+    using value_type = typename T::value_type;
+    t *= (value_type)(t.dimension());
+}
+
+// an std::vector
+template <class T, std::enable_if_t<is_std_vector<T>::value, bool> = true>
+void times_dimension(T& t)
+{
+    // do nothing
+}
+
+}

docs/source/examples/sfinae/python.cpp

@@ -0,0 +1,11 @@
+#include "mymodule.hpp"
+#include <pybind11/pybind11.h>
+#define FORCE_IMPORT_ARRAY
+#include <xtensor-python/pyarray.hpp>
+
+PYBIND11_MODULE(mymodule, m)
+{
+    xt::import_numpy();
+    m.doc() = "Test module for xtensor python bindings";
+    m.def("times_dimension", &mymodule::times_dimension<xt::pyarray<double>>);
+}

docs/source/index.rst

@@ -16,7 +16,7 @@ What are ``xtensor`` and ``xtensor-python``?
  - ``xtensor`` is a C++ library for multi-dimensional arrays enabling numpy-style broadcasting and lazy computing.
  - ``xtensor-python`` enables inplace use of numpy arrays with all the benefits from ``xtensor``
 
-     - C++ universal functions and broadcasting 
+     - C++ universal functions and broadcasting
      - STL - compliant APIs.
 
 
@@ -62,6 +62,7 @@ This software is licensed under the BSD-3-Clause license. See the LICENSE file f
    basic_usage
    array_tensor
    numpy_capi
+   examples
    cookiecutter
 
 .. toctree::
@@ -79,7 +80,7 @@ This software is licensed under the BSD-3-Clause license. See the LICENSE file f
 
 .. _NumPy: http://www.numpy.org
 .. _`Buffer Protocol`: https://docs.python.org/3/c-api/buffer.html
-.. _`numpy to xtensor cheat sheet`: http://xtensor.readthedocs.io/en/latest/numpy.html 
+.. _`numpy to xtensor cheat sheet`: http://xtensor.readthedocs.io/en/latest/numpy.html
 .. _xtensor: https://github.com/xtensor-stack/xtensor
-.. _pybind11: https://github.com/pybind/pybind11 
-.. _xtensor-python-cookiecutter: https://github.com/xtensor-stack/xtensor-python-cookiecutter 
+.. _pybind11: https://github.com/pybind/pybind11
+.. _xtensor-python-cookiecutter: https://github.com/xtensor-stack/xtensor-python-cookiecutter