Re: [Pkg-octave-devel] Re: Debian packages for Octave 2.9.10 available in experimental

[David Bateman]

Rafael Laboissiere wrote:
> * David Bateman <address@hidden> [2007-03-29 21:01]:
>
>   
>> Hum, I get Nan +Nani for (1+1i)/0 with g++ 4.0.1...
>>     
>
> I have here:
>
> $ g++ --version | fgrep g++
> g++ (GCC) 4.1.2 20061115 (prerelease) (Debian 4.1.1-21)
>   
This might play a part in the issue.

>   
>> Ok I think this
>>
>> http://www.nabble.com/inconsistent--a-b-.%5Ex-vs--a%5Ex-b%5Ex--tf3320209.html#a9230186
>>
>> thread is related to this issue. Given that I understand that your
>> result isn't the same as mine, However that "sparse(1+1i)./sparse(0)" is
>> not the same as "(1+1i)/0" I don't understand. They both use the C++
>> complex class underneath for the division itself...
>>     
>
> Do you think that this may be related to the specific libraries we use?
>   
No I don't think this matters, as the "./" operator directly calls the
libc++ std::<complex>(double)  division operator. However it should do
this for both "./" on full and sparse matrices, therefore "(1+1i)/0" and
"sparse(1+1i)/sparse(0)" should both give the same result. Can you
confirm what "sparse(1+1i)/sparse(0)" gives on your system? However,
given the "make check" failure you sent it appears that the full and
sparse divisions don't agree. Looking at the code is Sparse-op-defs.h
for quotient, I see no reason for this.

D.

-- 
David Bateman                                address@hidden
Motorola Labs - Paris                        +33 1 69 35 48 04 (Ph) 
Parc Les Algorithmes, Commune de St Aubin    +33 6 72 01 06 33 (Mob) 
91193 Gif-Sur-Yvette FRANCE                  +33 1 69 35 77 01 (Fax) 

The information contained in this communication has been classified as: 

[x] General Business Information 
[ ] Motorola Internal Use Only 
[ ] Motorola Confidential Proprietary