https://code.dlang.org/packages/gged/~main

I have been working on extending the functionality of multidimensional arrays using mir.ndslice, and now einstein summation is available.

It may not be practical, but it can be described nicely, like below, so I will introduce this.

    auto A  = iota(9.).array.gged!double(3,3);
    auto x  = iota(3.).gged!double(3);
    auto Ax = Einsum | A.ij*x.j; // [5, 14, 23]

On Monday, 16 January 2023 at 10:00:53 UTC, Zenw wrote:

    auto A  = iota(9.).array.gged!double(3,3);
    auto x  = iota(3.).gged!double(3);
    auto Ax = Einsum | A.ij*x.j; // [5, 14, 23]

Very nice. Do you plan to make it @nogc nothrow like mir-algorithm?

On Wednesday, 18 January 2023 at 21:32:50 UTC, klknn wrote:

    auto A  = iota(9.).array.gged!double(3,3);
    auto x  = iota(3.).gged!double(3);
    auto Ax = Einsum | A.ij*x.j; // [5, 14, 23]

Thank you for your comment. I studied the code of the mir-algorithm and now managed to make it @nogc nothrow. :)