Рассмотрите возможность абстрагирования доступа к данным с помощью Proxy (аналогично итераторам / умным указателям в C ++). К сожалению, синтаксис не так чист, как в C ++, так как operator () недоступен для перегрузки, а operator [] имеет одно аргумент, но все еще закрывается.
Конечно, этот дополнительный уровень абстракции добавляет сложности и самостоятельной работы, но он позволит вам внести минимальные изменения в существующий код, который использует объекты double [,,], в то же время позволяя использовать один double [] массив как для взаимодействия, так и для ваших вычислений в C #.
class Matrix3
{
// referece-to-element object
public struct Matrix3Elem{
private Matrix3Impl impl;
private uint dim0, dim1, dim2;
// other constructors
Matrix3Elem(Matrix3Impl impl_, uint dim0_, uint dim1_, uint dim2_) {
impl = impl_; dim0 = dim0_; dim1 = dim1_; dim2 = dim2_;
}
public double Value{
get { return impl.GetAt(dim0,dim1,dim2); }
set { impl.SetAt(dim0, dim1, dim2, value); }
}
}
// implementation object
internal class Matrix3Impl
{
private double[] data;
uint dsize0, dsize1, dsize2; // dimension sizes
// .. Resize()
public double GetAt(uint dim0, uint dim1, uint dim2) {
// .. check bounds
return data[ (dim2 * dsize1 + dim1) * dsize0 + dim0 ];
}
public void SetAt(uint dim0, uint dim1, uint dim2, double value) {
// .. check bounds
data[ (dim2 * dsize1 + dim1) * dsize0 + dim0 ] = value;
}
}
private Matrix3Impl impl;
public Matrix3Elem Elem(uint dim0, uint dim1, uint dim2){
return new Matrix2Elem(dim0, dim1, dim2);
}
// .. Resize
// .. GetLength0(), GetLength1(), GetLength1()
}
И затем использование этого типа для чтения и записи - 'foo [1,2,3]' теперь записывается как 'foo.Elem (1,2,3) .Value', как для чтения значений, так и для записи. значения в левой части выражений присваивания и значения.
void normalize(Matrix3 m){
double s = 0;
for (i = 0; i < input.GetLength0; i++)
for (j = 0; j < input.GetLength(1); j++)
for (k = 0; k < input.GetLength(2); k++)
{
s += m.Elem(i,j,k).Value;
}
for (i = 0; i < input.GetLength0; i++)
for (j = 0; j < input.GetLength(1); j++)
for (k = 0; k < input.GetLength(2); k++)
{
m.Elem(i,j,k).Value /= s;
}
}
Опять же, добавлены затраты на разработку, но обмениваются данными, устраняя накладные расходы на копирование и связанные с этим расходы на разработку. Это компромисс.