On Thursday, 2 April 2015 at 09:09:04 UTC, Martin Nowak wrote:
I'm not asking for a linear algebra library in phobos, although
we need one in dub and should consider having one in Phobos at
some point too.
But it would be nice if std.numeric came with a
multiply(T)(T[][] a, T[][] b, T[][] result) to complement
dotProduct and the built-in vector ops.
As you've seen dotProduct set us far ahead of C.
Sounds good:
- a matrix multiplication in phobos that is faster than a naive
implementation but does not try to beat BLAS gemm.
- a dub LA library around BLAS. There is
http://code.dlang.org/packages/cblas but it seems to be not much
more than dstep + cblas.h
How hard would it be to turn that into a lib.
Don't know, I perhaps try.
I have a working prototype of a BFGS non-linear optimizer
currently calling BLAS via array.ptr and the vague plan to make
it opensource. If that ever happens, a LA library around BLAS
would be a great co-product.
