|
From: Akshay S. <ak...@us...> - 2013-01-27 03:14:32
|
This is an automated email from the git hooks/post-receive script. It was
generated because a ref change was pushed to the repository containing
the project "matlisp".
The branch, tensor has been updated
via 7dc82f982eff648874c67454a6f6e5a947104f3d (commit)
from 8230a46f93849cdb60ce09173ab31687f998d07d (commit)
Those revisions listed above that are new to this repository have
not appeared on any other notification email; so we list those
revisions in full, below.
- Log -----------------------------------------------------------------
commit 7dc82f982eff648874c67454a6f6e5a947104f3d
Author: Akshay Srinivasan <aks...@gm...>
Date: Sat Jan 26 19:11:06 2013 -0800
Fixed bug with (:c :h) handling in GEMV.
diff --git a/src/level-2/gemv.lisp b/src/level-2/gemv.lisp
index c7ad77a..babc897 100644
--- a/src/level-2/gemv.lisp
+++ b/src/level-2/gemv.lisp
@@ -91,15 +91,16 @@
(type symbol job))
(if (member job '(:n :t))
(complex-base-typed-gemv! alpha A x beta y job)
- ;;The CBLAS way.
- (let-typed ((cx (let-typed ((ret (apply #'make-real-tensor (lvec->list (dimensions x))) :type complex-vector))
- (complex-typed-axpy! #c(-1d0 0d0) x ret))
- :type complex-vector))
- (complex-typed-num-scal! #c(-1d0 0d0) (tensor-realpart~ y))
+ ;;
+ (let-typed ((cx (let ((ret (complex-typed-copy! x (complex-typed-zeros (dimensions x)))))
+ (real-typed-num-scal! -1d0 (tensor-imagpart~ ret))
+ ret) :type complex-vector)
+ (y.view (tensor-imagpart~ y)))
+ (real-typed-num-scal! -1d0 y.view)
(complex-base-typed-gemv! (cl:conjugate alpha) A cx
(cl:conjugate beta) y (ecase job (:h :t) (:c :n)))
- (complex-typed-num-scal! #c(-1d0 0d0) (tensor-realpart~ y))
- y)))
+ (real-typed-num-scal! -1d0 y.view)))
+ y)
;;Symbolic
#+maxima
-----------------------------------------------------------------------
Summary of changes:
src/level-2/gemv.lisp | 15 ++++++++-------
1 files changed, 8 insertions(+), 7 deletions(-)
hooks/post-receive
--
matlisp
|