Re: [Gcl-devel] blas and complex

[Camm Maguire]

Greetings!

Vanuxem Grégory <address@hidden> writes:

> Hi,
> 
> Le samedi 23 juillet 2005 à 11:20 -0400, Camm Maguire a écrit :
> > Greetings!
> > 
> > Right now you need to double the array of the same type as the
> > real-part of the comlex, i.e. (make-array 20 :element-type long-float)
> > -> complex double blas vector of dimension 10.  Eventually  we can
> > make use of the gcc complex extensions.
> > 
> > 2d arrays hav the leading dimension stored in x->a.a_dims[1] (or
> > [0]).  This gives the blas 'lda'. 
> 
> I use another interface the user can set lda 
> (since lda can be superior).
> 
> > 
> > Take care,
> > 
> > Vanuxem Grégory <address@hidden> writes:
> > 
> > > Hi,
> > > 
> > > I run into problem with complex version of blas.
> > > How do you think someone can implement an interface of 
> > > complex version of blas.
> > > 
> > > For example:
> > > 
> > >     alpha * X + Y
> > > 
> > > how can i send alpha to c code ?
> > > A vector of two elements ?
> And here ?
> problematic no ?
> 

Yes, unfortunately the wrapper here will have to make the two float
bodies contiguous, which is easy enough but tedious.  Lisp comlex
numbers must have the same real and imag parts, which would seem to
indicate that they could be stored together, except for integer and
ratio which are both permitted in the same complex -- hence GCL uses
object pointers for both.

Take care,

> Cheers,
> 
> Greg
> 
> > > 
> > > And X, 2D-array of long-float ?
> > > 
> > > Cheers,
> > > 
> > > Greg
> > > 
> > > PS: there is an interface for cmucl or allegro named matlisp.
> > > It's use CLOS.
> > > 
> > > 
> > > 
> > > 
> > > 
> > > 
> > > _______________________________________________
> > > Gcl-devel mailing list
> > > address@hidden
> > > http://lists.gnu.org/mailman/listinfo/gcl-devel
> > > 
> > > 
> > > 
> > 
> 
> 
> 
> 
> _______________________________________________
> Gcl-devel mailing list
> address@hidden
> http://lists.gnu.org/mailman/listinfo/gcl-devel
> 
> 
> 

-- 
Camm Maguire                                            address@hidden
==========================================================================
"The earth is but one country, and mankind its citizens."  --  Baha'u'llah