From: SourceForge.net <no...@so...> - 2009-12-14 22:47:55
|
Bugs item #2913614, was opened at 2009-12-13 05:40 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2913614&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Share Libraries Group: None >Status: Pending >Resolution: Fixed Priority: 5 Private: No Submitted By: Andrej Vodopivec (andrejv) Assigned to: Nobody/Anonymous (nobody) Summary: LAPACK: dgesvd is broken Initial Comment: (%i4) A:genmatrix(lambda([i,j], random(2.0)), 4, 4)$ (%i5) dgesvd(A); Maxima encountered a Lisp error: Error during processing of --eval option "(cl-user::run)": The value "G" is not of type (SIMPLE-ARRAY DOUBLE-FLOAT (*)). Automatically continuing. To enable the Lisp debugger set *debugger-hook* to nil. (%i6) build_info(); Maxima version: 5.20post Maxima build date: 8:37 12/12/2009 Host type: i386-apple-darwin9.8.0 Lisp implementation type: SBCL Lisp implementation version: 1.0.19 ---------------------------------------------------------------------- >Comment By: Raymond Toy (rtoy) Date: 2009-12-14 17:47 Message: Fixed in CVS. I'd appreciate if you could test this out. Just need the updated version of dlascl.lisp. ---------------------------------------------------------------------- Comment By: Raymond Toy (rtoy) Date: 2009-12-14 12:56 Message: Thank you for the backtrace. I looks like a bug in the translation. I'll look into it. ---------------------------------------------------------------------- Comment By: Andrej Vodopivec (andrejv) Date: 2009-12-14 09:57 Message: The best I could come up with is attached. (%i3) dgesvd(A); debugger invoked on a SIMPLE-ERROR: Error during processing of --eval option "(cl-user::run)": The value "G" is not of type (SIMPLE-ARRAY DOUBLE-FLOAT (*)). Type HELP for debugger help, or (SB-EXT:QUIT) to exit from SBCL. restarts (invokable by number or by possibly-abbreviated name): 0: [MACSYMA-QUIT] Maxima top-level 1: [CONTINUE ] Ignore and continue with next --eval option. 2: [ABORT ] Skip rest of --eval options. 3: Skip to toplevel READ/EVAL/PRINT loop. 4: [QUIT ] Quit SBCL (calling #'QUIT, killing the process). ((FLET #:LAMBDA139) #<TYPE-ERROR {132837C1}>) 0] backtrace 20 0: ((FLET #:LAMBDA139) #<TYPE-ERROR {132837C1}>) 1: ((FLET #:LAMBDA139) #<TYPE-ERROR {132837C1}>)[:EXTERNAL] 2: (SIGNAL #<TYPE-ERROR {132837C1}>)[:EXTERNAL] 3: (ERROR TYPE-ERROR)[:EXTERNAL] 4: (SB-KERNEL::OBJECT-NOT-TYPE-ERROR-HANDLER #<unavailable argument> #.(SB-SYS:INT-SAP #X011FF7D8) #<SB-ALIEN-INTERNALS:ALIEN-VALUE :SAP #X0006A000 :TYPE (* (SB-ALIEN:STRUCT SB-VM::OS-CONTEXT-T-STRUCT))> (142 14)) 5: (SB-KERNEL::OBJECT-NOT-TYPE-ERROR-HANDLER #<unavailable argument> #.(SB-SYS:INT-SAP #X011FF7D8) #<SB-ALIEN-INTERNALS:ALIEN-VALUE :SAP #X0006A000 :TYPE (* (SB-ALIEN:STRUCT SB-VM::OS-CONTEXT-T-STRUCT))> (142 14))[:EXTERNAL] 6: (SB-KERNEL:INTERNAL-ERROR #.(SB-SYS:INT-SAP #X0006A000) #<unavailable argument>) 7: ("foreign function: call_into_lisp") 8: ("foreign function: funcall2") 9: ("foreign function: interrupt_internal_error") 10: ("foreign function: signal_emulation_wrapper") 11: ("foreign function: os_get_runtime_executable_path") 12: ("foreign function: os_get_runtime_executable_path") 13: (LAPACK::DLASCL "G" 0 0 2.926974774165864d0 9.989595361011175d145 7 1 #(2.926974774165864d0 2.757202750978322d0 1.8889138843296875d0 -0.5641183706403035d0 1.1901286134034068d0 0.2776359985260924d0 0.24215003810892827d0 0.0d0 4.0d0 0.3701265560273339d0 -0.35596118043153524d0 0.0d0 ...) 7 #<unavailable argument>) 14: (LAPACK::DLASQ1 4 #(2.926974774165864d0 1.8889138843296875d0 1.1901286134034068d0 0.24215003810892827d0) #(2.757202750978322d0 -0.5641183706403035d0 0.2776359985260924d0 0.0d0 2.926974774165864d0 2.757202750978322d0 1.8889138843296875d0 -0.5641183706403035d0 1.1901286134034068d0 0.2776359985260924d0 0.24215003810892827d0 0.0d0 ...) #(2.926974774165864d0 2.757202750978322d0 1.8889138843296875d0 -0.5641183706403035d0 1.1901286134034068d0 0.2776359985260924d0 0.24215003810892827d0 0.0d0 4.0d0 0.3701265560273339d0 -0.35596118043153524d0 0.0d0 ...) #<unavailable argument>) 15: (LAPACK::DBDSQR "U" 4 0 0 0 #(2.926974774165864d0 1.8889138843296875d0 1.1901286134034068d0 0.24215003810892827d0) #(2.757202750978322d0 -0.5641183706403035d0 0.2776359985260924d0 0.0d0 2.926974774165864d0 2.757202750978322d0 1.8889138843296875d0 -0.5641183706403035d0 1.1901286134034068d0 0.2776359985260924d0 0.24215003810892827d0 0.0d0 ...) #(0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 ...) 4 #(0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 ...) 4 #(0.0d0) 1 #(2.926974774165864d0 2.757202750978322d0 1.8889138843296875d0 -0.5641183706403035d0 1.1901286134034068d0 0.2776359985260924d0 0.24215003810892827d0 0.0d0 4.0d0 0.3701265560273339d0 -0.35596118043153524d0 0.0d0 ...) #<unavailable argument>) 16: (LAPACK:DGESVD "No columns of U" "No columns of V^T" 4 4 #(-2.926974774165864d0 0.10610322483770288d0 0.21267407794367488d0 0.41800896708772406d0 2.757202750978322d0 1.8889138843296875d0 0.14591288303296437d0 -0.0033539118023266363d0 0.46269703294122566d0 -0.5641183706403035d0 -1.1901286134034068d0 -0.16003447244391097d0 ...) 4 #(2.926974774165864d0 1.8889138843296875d0 1.1901286134034068d0 0.24215003810892827d0) #(0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 ...) 4 #(0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 0.0d0 ...) 4 #(2.757202750978322d0 -0.5641183706403035d0 0.2776359985260924d0 0.0d0 2.926974774165864d0 2.757202750978322d0 1.8889138843296875d0 -0.5641183706403035d0 1.1901286134034068d0 0.2776359985260924d0 0.24215003810892827d0 0.0d0 ...) 268 #<unavailable argument>) 17: ($DGESVD (($MATRIX SIMP) ((MLIST SIMP) 1.8276191992257917d0 0.39036923559557746d0 0.9647810497032392d0 0.7108801143185723d0) ((MLIST SIMP) 0.5044777533707618d0 0.8311254454429657d0 1.2251477785076492d0 1.6836253106718426d0) ((MLIST SIMP) 1.0111788892876237d0 1.356923990037029d0 0.07278663179984379d0 1.2505950366836145d0) ((MLIST SIMP) 1.9874629157389636d0 0.08317849230889962d0 1.2652507787332024d0 1.186304363033953d0)) NIL NIL) 18: (MEVAL1 #<unavailable argument>) 19: (MEVAL (($DGESVD) |$a|)) ---------------------------------------------------------------------- Comment By: Raymond Toy (rtoy) Date: 2009-12-14 08:57 Message: I cannot reproduce this either with cmucl on mac os x. Running this with maxima -g and looking at the Lisp error will help a lot in figuring this problem out. The boolean results are documented; it means the left and right singular vectors were not requested. Perhaps that should be changed, but that's a different issue. ---------------------------------------------------------------------- Comment By: Dieter Kaiser (crategus) Date: 2009-12-13 12:39 Message: On my system I get: (%i2) build_info(); Maxima version: 5.20post Maxima build date: 17:36 12/13/2009 Host type: i686-pc-linux-gnu Lisp implementation type: CLISP Lisp implementation version: 2.44.1 (2008-02-23) (built 3436700604) (memory 3469710995) (%o2) (%i3) A:genmatrix(lambda([i,j], random(2.0)), 4, 4)$ (%i4) dgesvd(A); (%o4) [[3.755933922081446, 1.830035687063915, .7142611635301948, 0.127875116930834], false, false] I do not get a Lisp error. But the result contains boolean values. Dieter Kaiser ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2913614&group_id=4933 |