On Monday, 26 September 2016 at 12:05:46 UTC, Martin Tschierschke
On Friday, 23 September 2016 at 13:25:30 UTC, Ilya Yaroshenko
During my time studying mathematics - in the last century, one
Dr. of the mathematical department has developed a C++ library
as glue to some fast C libs using the possibility of operator
overloading, to make the formulation of matrix algorithms quite
So would it be possible to get:
// Performs: c := alpha a x b + beta c
// glas is a pointer to a GlasContext
glas.gemm(alpha, a, b, beta, c);
transformed to something like:
c = alpha * A * B + beta * c;
Maybe not for performance but for readability purpose?
There was a recent thread about operator overloading but I
didn't understood the problem.
(sorry, still a D newbie).
Yes, this is possible. But, first we need a more low level
interface. It will help to extend the Mir library and to write
bindings for other languages (Julia, Python, R).