1 files changed, 60 insertions, 0 deletions
diff --git a/doc/source/reference/c-api/iterator.rst b/doc/source/reference/c-api/iterator.rst
index b4adaef9b..07187e7f1 100644
--- a/doc/source/reference/c-api/iterator.rst
+++ b/doc/source/reference/c-api/iterator.rst
@@ -203,6 +203,66 @@ is used to control the memory layout of the allocated result, typically
     }
 
 
+Multi Index Tracking Example
+----------------------------
+
+This example shows you how to work with the :c:data:`NPY_ITER_MULTI_INDEX` flag. For simplicity, we assume the argument is a two-dimensional array.
+
+.. code-block:: c
+
+   int PrintMultiIndex(PyArrayObject *arr) {
+       NpyIter *iter;
+       NpyIter_IterNextFunc *iternext;
+       npy_intp multi_index[2];
+
+       iter = NpyIter_New(
+           arr, NPY_ITER_READONLY | NPY_ITER_MULTI_INDEX | NPY_ITER_REFS_OK,
+           NPY_KEEPORDER, NPY_NO_CASTING, NULL);
+       if (iter == NULL) {
+           return -1;
+       }
+       if (NpyIter_GetNDim(iter) != 2) {
+           NpyIter_Deallocate(iter);
+           PyErr_SetString(PyExc_ValueError, "Array must be 2-D");
+           return -1;
+       }
+       if (NpyIter_GetIterSize(iter) != 0) {
+           iternext = NpyIter_GetIterNext(iter, NULL);
+           if (iternext == NULL) {
+               NpyIter_Deallocate(iter);
+               return -1;
+           }
+           NpyIter_GetMultiIndexFunc *get_multi_index =
+               NpyIter_GetGetMultiIndex(iter, NULL);
+           if (get_multi_index == NULL) {
+               NpyIter_Deallocate(iter);
+               return -1;
+           }
+
+           do {
+               get_multi_index(iter, multi_index);
+               printf("multi_index is [%" NPY_INTP_FMT ", %" NPY_INTP_FMT "]\n",
+                      multi_index[0], multi_index[1]);
+           } while (iternext(iter));
+       }
+       if (!NpyIter_Deallocate(iter)) {
+           return -1;
+       }
+       return 0;
+   }
+
+When called with a 2x3 array, the above example prints:
+
+.. code-block:: sh
+
+   multi_index is [0, 0]
+   multi_index is [0, 1]
+   multi_index is [0, 2]
+   multi_index is [1, 0]
+   multi_index is [1, 1]
+   multi_index is [1, 2]
+
+
 Iterator Data Types
 ---------------------