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 |