Iterating over an mdspan

There are a few ways to iterate over all elements in a multidimensional span (std::mdspan). Here are a couple of common approaches:

Using Nested Loops

You can use nested loops to iterate over each dimension of the mdspan. The number of loops required will match the rank (number of dimensions) of the mdspan.In this example, we use two nested loops to iterate over a 2D mdspan.

The outer loop iterates over the rows using matrix.extent(0), and the inner loop iterates over the columns using matrix.extent(1).

1#include <iostream>
2#include <mdspan>
3
4int main() {
5  std::array<int, 12> arr{
6    1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12};
7
8  std::mdspan<int, std::extents<
9    size_t, 3, 4>> M {arr.data()};
10
11  // Iterating using nested loops
12  for (size_t i = 0; i < M.extent(0); ++i) {
13    for (size_t j = 0; j < M.extent(1); ++j) {
14      std::cout << M[i, j] << " ";  
15    }
16    std::cout << "\n";
17  }
18}

11 2 3 4
25 6 7 8
39 10 11 12

Using a Single Loop with Arithmetic

Alternatively, you can use a single loop to iterate over all elements by treating the mdspan as a flattened 1D array and using arithmetic to calculate the corresponding indices.In this approach, we use a single loop that iterates from 0 to matrix.size() - 1.

Inside the loop, we calculate the corresponding row and column indices using integer division and modulo operations.

1#include <iostream>
2#include <mdspan>
3
4int main() {
5  std::array<int, 12> arr{
6    1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12};
7  std::mdspan<int, std::extents<
8    size_t, 3, 4>> matrix{arr.data()};
9
10  // Iterating using a single loop with arithmetic
11  for (size_t i = 0; i < matrix.size(); ++i) {
12    size_t row = i / matrix.extent(1);  
13    size_t col = i % matrix.extent(1);  
14    std::cout << matrix[row, col] << " ";
15  }
16}

11 2 3 4 5 6 7 8 9 10 11 12

Both approaches yield the same result, so you can choose the one that suits your needs and preferences.