Re: Re: (DOT REAL-MAT COMPLEX-MAT) ERROR

Raymond Toy <[email protected]> 26 Jul 2002 17:39:34 -0400
Newsgroups gmane.lisp.matlisp.devel
Message-ID <[email protected]>
>>>>> "Raymond" == Raymond Toy <[email protected]> writes:

    Raymond> In any case, I think the following solves the immediate problem:

    Raymond> (defmethod dot ((x complex-matrix) (y complex-matrix) &optional (conjugate-p t))
    Raymond>   (let* ((nxm (number-of-elements x))
    Raymond> 	 (store-x (store x))
    Raymond> 	 (store-y (store y))
    Raymond> 	 (dot (if conjugate-p
    Raymond> 		  (zdotc nxm store-x 1 store-y 1)
    Raymond> 		  (zdotu nxm store-x 1 store-y 1))))
    Raymond>     (if (zerop (imagpart dot))
    Raymond> 	(realpart dot)
    Raymond> 	(complex (realpart dot) (imagpart dot)))))

    Raymond> I think this need to be fixed in ffi-cmu.lisp in some way.  There's no
    Raymond> telling what other bugs will be found because of this, if I'm right.

I've committed a fix to ffi-cmu.lisp.  I'd appreciate it if you can
try it out.  It works for me. :-)

I hope Tunc will try your example with ACL to see if the same bug
exists there.

Ray