[lisplab-cvs] r128 - src/matlisp

[Jørn Inge Vestgården]

Author: jivestgarden
Date: Fri Jan 15 05:53:20 2010
New Revision: 128

Log:
bugfix in eigenvalues

Modified:
   src/matlisp/geev.lisp

Modified: src/matlisp/geev.lisp
==============================================================================
--- src/matlisp/geev.lisp	(original)
+++ src/matlisp/geev.lisp	Fri Jan 15 05:53:20 2010
@@ -30,6 +30,10 @@
 ;;; 
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 
+;;; An issue with this code is the following: it creates output matrices 
+;;; of default kind, ie. matrix-dge and matrix-zge. It should create based 
+;;; on input in stead.
+
 (in-package :lisplab)
 
 (defmethod eigenvectors ((a matrix-blas-dge))
@@ -70,8 +74,8 @@
 
 (defun combine-ri-vectors (wr wi)
   (let* ((n (size wr))
-	 (wr2 (make-instance 'matrix-dge :rows n :cols 1 :store wr))
-	 (wi2 (make-instance 'matrix-dge :rows n :cols 1 :store wi)))
+	 (wr2 (make-instance 'matrix-dge :dim (list n 1) :store wr))
+	 (wi2 (make-instance 'matrix-dge :dim (list n 1) :store wi)))
     (if (.= wi2 0)
 	wr2
 	(.+ wr2 (.* %i (convert wi2 'matrix-zge))))))