matlisp-users

[Matlisp-users] Problem with function geev

[Thibault L. <tl...@di...>]

Hi,
Ray, I saw your memo on #lisp. I downloaded the csv version and
installed it without problem. But I continue to get the same behavior :-
( : when I use (geev ..) it starts computing and never ends. 
I am still using an "old" version of cmucl (19apre3) could it be the
cause of my problem ? 
What can I do to track the problem?

wok:/usr/local/src/matlisp$ cmucl
; Loading #p"/home/tl/init.lisp".
;; Loading #p"/home/tl/comp/lang/lisp/code/asdf/asdf.x86f".
CMU Common Lisp 19a-pre3, running on wok
With core: /usr/local/cmucl-19a/lib/cmucl/lib/lisp.core
Dumped on: Wed, 2004-06-30 18:13:25+01:00 on lorien
See <http://www.cons.org/cmucl/> for support information.
Loaded subsystems:
    Python 1.1, target Intel x86
    CLOS based on Gerd's PCL 2004/04/14 03:32:47
* (load "/usr/local/src/matlisp/start.lisp")

; Loading #p"/usr/local/src/matlisp/start.lisp".
; Converted SETLOGICALROOT.
; Converted GETLOGICALROOT.
; Converted DEFLOGICALPATH.
;; Loading #p"/usr/local/src/matlisp/system.dcl".
;;; Loading #p"/usr/local/src/matlisp/defsystem.lisp".

[...]

 ** MATLISP is loaded.  Type (HELP MATLISP)
    to see a list of available symbols.
    To use matlisp:

          (use-package "MATLISP")

    or

          (in-package "MATLISP-USER")

 ** The logical pathname matlisp has been
    set to:

          /usr/local/src/matlisp/

; T
* (in-package :matlisp-user)

#<The MATLISP-USER package, 1/9 internal, 0/2 external>
* (setq m (make-real-matrix '((1 1 2) (2 2 4) (3 3 6))))

Warning:  Declaring M special.

#<REAL-MATRIX  3 x 3
    1.0000      1.0000      2.0000     
    2.0000      2.0000      4.0000     
    3.0000      3.0000      6.0000     >
* (geev m)
Terminated
wok:/usr/local/src/matlisp$ 

I am not using the last version of cmucl:
 
CL-USER> (lisp-implementation-version)
"19a-pre3"


-- 
Thibault Langlois <tl...@di...>
FCUL / DI

RE: [Matlisp-users] Problem with function geev

[Ernst v. W. <ev...@in...>]

Hey, I have had the same problem with geev for years.  I use Allegro.  SVD
has always been a perfect replacement for me: it works and is more than fast
enough for my purposes.  Maybe wrong of me not to heve mentioned this
before.

Kind regards,

Ernst van Waning

> -----Oorspronkelijk bericht-----
> Van: mat...@li... 
> [mailto:mat...@li...] Namens 
> Thibault Langlois
> Verzonden: 27 July 2005 12:49
> Aan: mat...@li...
> Onderwerp: [Matlisp-users] Problem with function geev
> 
> Hi,
> Ray, I saw your memo on #lisp. I downloaded the csv version 
> and installed it without problem. But I continue to get the 
> same behavior :- ( : when I use (geev ..) it starts computing 
> and never ends. 
> I am still using an "old" version of cmucl (19apre3) could it 
> be the cause of my problem ? 
> What can I do to track the problem?
> 
> wok:/usr/local/src/matlisp$ cmucl
> ; Loading #p"/home/tl/init.lisp".
> ;; Loading #p"/home/tl/comp/lang/lisp/code/asdf/asdf.x86f".
> CMU Common Lisp 19a-pre3, running on wok With core: 
> /usr/local/cmucl-19a/lib/cmucl/lib/lisp.core
> Dumped on: Wed, 2004-06-30 18:13:25+01:00 on lorien See 
> <http://www.cons.org/cmucl/> for support information.
> Loaded subsystems:
>     Python 1.1, target Intel x86
>     CLOS based on Gerd's PCL 2004/04/14 03:32:47
> * (load "/usr/local/src/matlisp/start.lisp")
> 
> ; Loading #p"/usr/local/src/matlisp/start.lisp".
> ; Converted SETLOGICALROOT.
> ; Converted GETLOGICALROOT.
> ; Converted DEFLOGICALPATH.
> ;; Loading #p"/usr/local/src/matlisp/system.dcl".
> ;;; Loading #p"/usr/local/src/matlisp/defsystem.lisp".
> 
> [...]
> 
>  ** MATLISP is loaded.  Type (HELP MATLISP)
>     to see a list of available symbols.
>     To use matlisp:
> 
>           (use-package "MATLISP")
> 
>     or
> 
>           (in-package "MATLISP-USER")
> 
>  ** The logical pathname matlisp has been
>     set to:
> 
>           /usr/local/src/matlisp/
> 
> ; T
> * (in-package :matlisp-user)
> 
> #<The MATLISP-USER package, 1/9 internal, 0/2 external>
> * (setq m (make-real-matrix '((1 1 2) (2 2 4) (3 3 6))))
> 
> Warning:  Declaring M special.
> 
> #<REAL-MATRIX  3 x 3
>     1.0000      1.0000      2.0000     
>     2.0000      2.0000      4.0000     
>     3.0000      3.0000      6.0000     >
> * (geev m)
> Terminated
> wok:/usr/local/src/matlisp$ 
> 
> I am not using the last version of cmucl:
>  
> CL-USER> (lisp-implementation-version)
> "19a-pre3"
> 
> 
> --
> Thibault Langlois <tl...@di...>
> FCUL / DI
> 
> 
> 
> -------------------------------------------------------
> SF.Net email is sponsored by: Discover Easy Linux Migration Strategies
> from IBM. Find simple to follow Roadmaps, straightforward articles,
> informative Webcasts and more! Get everything you need to get up to
> speed, fast. http://ads.osdn.com/?ad_id=7477&alloc_id=16492&op=click
> _______________________________________________
> Matlisp-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matlisp-users
> 

Re: [Matlisp-users] Problem with function geev

[Raymond T. <ray...@er...>]

>>>>> "Thibault" == Thibault Langlois <tl...@di...> writes:

    Thibault> Ray, I saw your memo on #lisp. I downloaded the csv version and
    Thibault> installed it without problem. But I continue to get the same behavior :-
    Thibault> ( : when I use (geev ..) it starts computing and never ends. 
    Thibault> I am still using an "old" version of cmucl (19apre3) could it be the
    Thibault> cause of my problem ? 

Hmm.  I guess I'll have to try on x86.  19apre3 (why not just 19a?)
shouldn't really matter.

    Thibault> What can I do to track the problem?

    Thibault> * (geev m)
    Thibault> Terminated

Why does it say "Terminated"?  It dies by itself?  You killed it?
Linux's out-of-memory killer killed it?  Is there a core dump?

Ray

Re: Re: [Matlisp-users] Problem with function geev

[Thibault L. <tl...@di...>]

On Wed, 2005-07-27 at 08:48 -0400, Raymond Toy wrote:
> >>>>> "Thibault" == Thibault Langlois <tl...@di...> writes:
> 
>     Thibault> Ray, I saw your memo on #lisp. I downloaded the csv version and
>     Thibault> installed it without problem. But I continue to get the same behavior :-
>     Thibault> ( : when I use (geev ..) it starts computing and never ends. 
>     Thibault> I am still using an "old" version of cmucl (19apre3) could it be the
>     Thibault> cause of my problem ? 
> 
> Hmm.  I guess I'll have to try on x86.  19apre3 (why not just 19a?)

No particular reason except that I have been lazy ...

> shouldn't really matter.
> 
>     Thibault> What can I do to track the problem?
> 
>     Thibault> * (geev m)
>     Thibault> Terminated
> 
> Why does it say "Terminated"?  It dies by itself?  You killed it?
> Linux's out-of-memory killer killed it?  Is there a core dump?
> 

Sorry, I should have explained that. I killed the process after it ran
for several minutes using 100% of cpu.
 
> Ray
> 

Re: [Matlisp-users] Problem with function geev

[Raymond T. <ray...@er...>]

>>>>> "Thibault" == Thibault Langlois <tl...@di...> writes:

    Thibault> #<The MATLISP-USER package, 1/9 internal, 0/2 external>
    Thibault> * (setq m (make-real-matrix '((1 1 2) (2 2 4) (3 3 6))))

    Thibault> Warning:  Declaring M special.

    Thibault> #<REAL-MATRIX  3 x 3
    Thibault>     1.0000      1.0000      2.0000     
    Thibault>     2.0000      2.0000      4.0000     
    Thibault>     3.0000      3.0000      6.0000     >
    Thibault> * (geev m)

Just tried this (finally) using cvs matlisp and cmucl 2005-07.
(geev m) returns

#<REAL-MATRIX  3 x 1
    0.0000     
    9.0000     
   -5.24537E-16 >

which is pretty close to what matlab returns.

I'll have to look at the code some more.  Perhaps it's doing something
wrong.  (There used to be a memory corruption error in some eigenvalue
routine, but that's been fixed.)

Ray

Re: [Matlisp-users] Problem with function geev

[Thibault L. <tl...@di...>]

On Wed, 2005-07-27 at 16:29 -0400, Raymond Toy wrote: 
> >>>>> "Thibault" == Thibault Langlois <tl...@di...> writes:
> 
>     Thibault> #<The MATLISP-USER package, 1/9 internal, 0/2 external>
>     Thibault> * (setq m (make-real-matrix '((1 1 2) (2 2 4) (3 3 6))))
> 
>     Thibault> Warning:  Declaring M special.
> 
>     Thibault> #<REAL-MATRIX  3 x 3
>     Thibault>     1.0000      1.0000      2.0000     
>     Thibault>     2.0000      2.0000      4.0000     
>     Thibault>     3.0000      3.0000      6.0000     >
>     Thibault> * (geev m)
> 
> Just tried this (finally) using cvs matlisp and cmucl 2005-07.
> (geev m) returns
> 
> #<REAL-MATRIX  3 x 1
>     0.0000     
>     9.0000     
>    -5.24537E-16 >
> 
> which is pretty close to what matlab returns.
> 
> I'll have to look at the code some more.  Perhaps it's doing something
> wrong.  (There used to be a memory corruption error in some eigenvalue
> routine, but that's been fixed.)
> 
> Ray
> 

I have been trying to find where this bug occurs. First surprise, it
seems to be in fortran code, somewhere in the DLAMCH routine. I tried to
put some "PRINT *, 'some string'" in the code in order to localize the
point. Doing so I ended with a version of the file that make matlisp
work correctly ! 

Then I went to #lisp to see if someone involved was online. Consulting
the log at http://meme.b9.com I noticed that someone (aeacides) was also
looking at the same file and at the same point in the file ...

So, if I add the prints arround the loop, it works:

         PRINT *, 'before loop 2'
*+       WHILE( C.EQ.A )LOOP
   20    CONTINUE
         IF( C.EQ.A ) THEN
            B = 2*B
            C = DLAMC3( A, B )
            GO TO 20
         END IF
*+       END WHILE
         PRINT *, 'after loop 2'

if I comment the first and last lines it does not work. 

It looks like a nasty memory corruption bug. 

Thibault

PS: the problem is not related to the geev function and should affect
other matlisp functions.