Re: Double-precision floating point arithmetics behaves differently in Guile and in C

Panicz Maciej Godek schreef op ma 27-09-2021 om 16:04 [+0200]: > niedz., 26 wrz 2021 o 21:11 Matt Wette > napisał(a): > > > Your post is a little vague.  Are you comparing singular values > > vectors > > of A vs > > eigenvalues, eigenvectors of A'*A ?  Also, SVD is iterative and > > different alg

Re: Double-precision floating point arithmetics behaves differently in Guile and in C

niedz., 26 wrz 2021 o 21:11 Matt Wette napisał(a): > Your post is a little vague. Are you comparing singular values vectors > of A vs > eigenvalues, eigenvectors of A'*A ? Also, SVD is iterative and > different algorithms may yield > different results (and some of the SVD algorithms out there a

Re: Double-precision floating point arithmetics behaves differently in Guile and in C

Your post is a little vague.  Are you comparing singular values vectors of A vs eigenvalues, eigenvectors of A'*A ?  Also, SVD is iterative and different algorithms may yield different results (and some of the SVD algorithms out there are not great).  Are you using the exact same algorithm but n

Re: Double-precision floating point arithmetics behaves differently in Guile and in C

Thanks for the feedback! It turned out that SVD was fine and I was looking in the wrong place -- in the C counterpart, I was feeding my ellipsoid algorithm with uninitialized junk, and the compiler didn't shout at me. Sorry for the noise! niedz., 26 wrz 2021 o 15:16 Eli Zaretskii napisał(a): >

Re: Double-precision floating point arithmetics behaves differently in Guile and in C

> From: Panicz Maciej Godek > Date: Sun, 26 Sep 2021 13:55:20 +0200 > > I forgot to mention that I run Guile in an Emacs session running in a WSL > console on Windows 10. > The tests of the C code that I've been performing so far were executed in > an MSYS terminal, but I have just tried running

Re: Double-precision floating point arithmetics behaves differently in Guile and in C

On Sun, Sep 26, 2021 at 10:45:51AM +0200, Panicz Maciej Godek wrote: > Hi, > > I recently ran into a very weird problem. I am using code for computing > Singular Value Decomposition that's been ported do Scheme by Gerald Sussman > as a part of the scmutils package (more specifically, I use some bi

Re: Double-precision floating point arithmetics behaves differently in Guile and in C

I forgot to mention that I run Guile in an Emacs session running in a WSL console on Windows 10. The tests of the C code that I've been performing so far were executed in an MSYS terminal, but I have just tried running them in WSL console, and the radii get ridiculous values. While the values used

Double-precision floating point arithmetics behaves differently in Guile and in C

Hi, I recently ran into a very weird problem. I am using code for computing Singular Value Decomposition that's been ported do Scheme by Gerald Sussman as a part of the scmutils package (more specifically, I use some bits of code extracted from guile-scmutils