Commit acdc5d17 authored by Liam Healy's avatar Liam Healy

BLAS1 completed (untested)

Change the generic specializer for scalars to :element-type (from
:element-c-type), and restore BLAS1 to gsll.asd.
parent bf94e20d
;; Definition of GSLL system
;; Liam Healy
;; Time-stamp: <2008-07-15 22:24:27EDT gsll.asd>
;; Time-stamp: <2008-08-03 23:27:56EDT gsll.asd>
;; $Id$
(asdf:defsystem "gsll"
......@@ -76,20 +76,20 @@
(:file "trigonometry" :depends-on (return-structures))
(:file "zeta" :depends-on (return-structures))))
(:file "sorting" :depends-on (init data))
#+no
(:module linear-algebra
:depends-on (init data)
:components
((:file "blas1")
(:file "blas2")
(:file "blas3" :depends-on (blas2))
(:file "lu")
(:file "qr")
(:file "qrpt")
(:file "svd")
(:file "cholesky")
(:file "diagonal")
(:file "householder")))
;(:file "blas2")
;(:file "blas3" :depends-on (blas2))
;(:file "lu")
;(:file "qr")
;(:file "qrpt")
;(:file "svd")
; (:file "cholesky")
; (:file "diagonal")
;(:file "householder")
))
#+no
(:file "eigensystems" :depends-on (init data))
;; Skip fft for now, I'm not sure how it works in C
......
;; Macro for defining GSL functions.
;; Liam Healy 2008-04-16 20:49:50EDT defmfun.lisp
;; Time-stamp: <2008-08-02 19:18:11EDT defmfun.lisp>
;; Time-stamp: <2008-08-03 23:27:25EDT defmfun.lisp>
;; $Id$
(in-package :gsl)
......@@ -302,7 +302,7 @@
(if (and replacing (listp arg))
(list (first arg)
(case (second arg)
(:element-c-type (number-class-from-type element-type))
(:element-type (number-class-from-type element-type))
(:component-float-type
(number-class-from-type (component-float-type element-type)))
(otherwise
......
;; BLAS level 1, Vector operations
;; Liam Healy, Wed Apr 26 2006 - 15:23
;; Time-stamp: <2008-08-02 20:05:38EDT blas1.lisp>
;; Time-stamp: <2008-08-03 23:24:43EDT blas1.lisp>
;; $Id$
(in-package :gsl)
......@@ -100,7 +100,7 @@
:documentation ; FDL
"Copy the elements of the vector x into the vector y.")
(defmfun scale ((alpha :element-c-type) (x vector))
(defmfun scale ((alpha :element-type) (x vector))
;; Alpha is the same type as the elements of vector, so for complex
;; vectors it must be complex.
("gsl_blas_" :type "scal")
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment