;;; Lisplab, vector1-d.lisp
;;; Level1, double float vectors

;;; Copyright (C) 2012 Joern Inge Vestgaarden
;;;
;;; This program is free software; you can redistribute it and/or modify
;;; it under the terms of the GNU General Public License as published by
;;; the Free Software Foundation; either version 2 of the License, or
;;; (at your option) any later version.
;;;
;;; This program is distributed in the hope that it will be useful,
;;; but WITHOUT ANY WARRANTY; without even the implied warranty of
;;; MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
;;; GNU General Public License for more details.
;;;
;;; You should have received a copy of the GNU General Public License along
;;; with this program; if not, write to the Free Software Foundation, Inc.,
;;; 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.

(in-package :lisplab)

;; (ll-def-vector1-class-and-vref vector-d element-double-float type-blas-store double-float)

(defclass vector-d (vector-base element-double-float)
   ((store :initarg :store
	   :initform nil
	   :reader vector-store
	   :type type-blas-store)))

(defmethod vref ((vector  vector-d) idx)
  (aref (the type-blas-store (slot-value vector 'store)) idx))

(defmethod (setf vref) (value (vector vector-d) idx)
  (let ((val2 (coerce value 'double-float)))
    (declare (type double-float val2))
    (setf (aref  (the type-blas-store (slot-value vector 'store)) idx)
	  val2)
    val2))