Re: [Matlisp-users] Re: Matrix column access

[Raymond T. <ray...@er...>]

>>>>> "Nicolas" == Nicolas Neuss <Nic...@iw...> writes:

    Nicolas> This makes it possible to have also matrices of single-floats, etc.

Oh, matlisp consciously supports only double-floats, but single-floats
would be possible.  Some kind of macro like you have would help.

    Nicolas> 2. One should automate Matlisp's BLAS interface.  Something like (get-blas
    Nicolas> 'daxpy) should make the Fortran daxpy routine available.
    >> 
    >> I don't follow what you mean by this.  They're always available in
    >> Matlisp.

    Nicolas> But not every routine as much as I remember.  And for
    Nicolas> complete generality, also the single-float versions
    Nicolas> should be available.  I think that it should be possible
    Nicolas> to automate the process of making such a routine
    Nicolas> available.

I think all the BLAS routines are there.  (Well, maybe not the level-3
blas routines).  

Not quite sure what you mean by automating the process of making the
routines available.  Someone has to figure out the function signature
to tell Lisp how to call it.  And some of the functions take
characters or strings, and that is very dependent on the Fortran
compiler.  Same with returning complex values.

I did toy around with the idea of letting f2cl grovel over the
functions and extracting the interfaces, but I didn't get very far
with that.

    Nicolas> interface.  What I don't like is the duplication of code
    Nicolas> for the different matrix classes.

Yes, that's not good.  It's been cleaned up some, but I think there's
still quite a bit of duplication.

Ray