On Tuesday, 23 December 2014 at 03:11:20 UTC, Laeeth Isharc wrote:
On Monday, 22 December 2014 at 22:46:57 UTC, aldanor wrote:
On Monday, 22 December 2014 at 22:36:16 UTC, H. S. Teoh via
Digitalmars-d wrote:
FYI, Kenji's merge has since been merged. So now the stage is
set for
somebody to step up and write a nice multidimensional array
implementation.
One important thing to wish for, in my opinion, is that the
design of such implementation would allow for (future
potential) integration with linear algebra libraries like
blas/lapack without having to be rewritten from scratch (e.g.
so it doesn't end up like Python's array module which got
completely superceded by numpy).
You mean especially for sparse matrices ? What is it that
needs to be borne in mind for regular matrices ?
The layout in lapck/blas is column major so it can be handy using
a wrapper around arrays to provide the FORTRAN indexing.
Also you need to pass the .ptr property of the array or &a[0]. D
arrays are fat and include their length.
Cheers,
uri
