C++23 (or later) succinct way to extract elements of an array using an index vector?

Question

I've been reviewing new ways of slicing arrays using the latest standards, but it's been a bit overwhelming and I have a question.

Is there any new, succinct way, with the new C++23 additions to ranges, views, mdspan, or even the ones to come submdspan, mdarray, to extract elements of an mdspan into another container using a vector with the indices to be extracted?

I have to work with std::array that I know their size beforehand, however, the indices will be known only at runtime so I can use temporary std::vector types or other containers.

Here's a godbolt example, reproduced below

#include <iostream>
#include <array>
#include <vector>
#include <https://raw.githubusercontent.com/kokkos/mdspan/single-header/mdspan.hpp>

int main()
{
    constexpr std::size_t pDeg = 16 ;
    std::array<std::size_t,pDeg*pDeg> arr;
    for (std::size_t i=0;i<arr.size();i++)
        arr[i] = i ;

    auto md = std::mdspan<std::size_t,std::extents<std::size_t,pDeg,pDeg>>(arr.data());

    std::vector<std::size_t> ind1 = {1,5,6,7,8,9,4,5};
    std::vector<std::size_t> ind2 = {14,11,12,13,11,1,5,6};

}

I have functions implemented to operate on matrices: multiply(), add(), etc

I would like to write matrix algebra expressions like

auto result = add(multiply(md[ind1,:],md[:,ind2]),md[ind1,ind2]);

where md[ind1,:] is a submatrix containing rows 1,5,6,7,8,9,4,5 of md, md[:,ind2] contains columns 14,11,12,13,11,1,5,6 of md and md[ind1,ind2] is a 8-by-8 submatrix of md containing all combinations of rows and columns given by ind1 and ind2 respectively (notice that indices can be repeated).

or something like

auto result = add(multiply(md.extract<0>(ind1),md.extract<1>(ind2)),md.extract<0,1>(ind1,ind2));

where result is always a std::array that I know the size at compile time.

The methods do not necessarily have to be the ones I propose, but they needs to be in a sort of short notation since there will be many matrix algebra operations.