Учитывая плотный домен dom: domain(n);, где n ‹ 3, объявление sps1: sparse subdomain(dom); дает разреженный поддомен sps1 из dom. С sps1 возможна обычная нарезка массива/матрицы. То есть, учитывая матрицу A: [sps1], можно взять n - 1-мерные срезы A. Однако обычная матричная операция transpose() неприменима.
Определение второй матрицы B:[sps2] над другим разреженным поддоменом sps2 = CSRDomain(dom) позволяет получить transpose()s из B, но возможность разрезать на B утрачивается.
Обе эти способности, казалось бы, должны быть всегда доступны. Есть ли лучший способ объявить разреженные поддомены, которые сохранят два?
