bugGNU Octave - Bugs: bug #36368, complex diagonal matrix powers...

 
 

bug #36368: complex diagonal matrix powers read uninit memory

Submitter:  None
Submitted:  Wed 02 May 2012 05:48:06 AM UTC
   
 
Category:  None Severity:  5 - Blocker
Priority:  7 - High Item Group:  Regression
Status:  Fixed Assigned to:  jordigh
Originator Name:  muslim Originator Email:  -email is unavailable-
Open/Closed:  * Closed Release:  * dev
Operating System:  * Any Fixed Release:  None
Planned Release:  None
* Mandatory Fields

Add a New Comment Rich Markup
   

Thu 03 May 2012 04:57:19 AM UTC, comment #4: 

Patched in this changeset:

http://hg.savannah.gnu.org/hgweb/octave/rev/6250e1232c9c

This should be part of the 3.6.2 bugfix release that should happen very soon now.

Jordi Gutiérrez Hermoso <jordigh>
Group Member
Thu 03 May 2012 12:42:48 AM UTC, comment #3: 

Er, sorry, I mean, I'll push this to the stable branch.

Jordi Gutiérrez Hermoso <jordigh>
Group Member
Thu 03 May 2012 12:38:06 AM UTC, comment #2: 

I got a patch for this. DiagArray2<T>::operator() didn't have a non-const version. I'm running some tests, and I'll push this soon to the development branch.

Jordi Gutiérrez Hermoso <jordigh>
Group Member
Wed 02 May 2012 06:10:41 AM UTC, comment #1: 

I also see this behavior with the current development sources (both stable and default branches).

This appears to be a regression in going from 3.2.x to 3.4.x.

John W. Eaton <jwe>
Group administrator
Wed 02 May 2012 05:48:06 AM UTC, original submission:  

Hi,

When I running the code below, I get wrong answer.


A=diag(exp(2*pi*1i*[0:2]/3))
A =

Diagonal Matrix

   1.00000 + 0.00000i                    0                    0
                    0  -0.50000 + 0.86603i                    0
                    0                    0  -0.50000 - 0.86603i

But:

A^2
ans =

Diagonal Matrix

 Columns 1 and 2:

   6.4406e-247 + 5.3118e+006i                            0
                            0   4.9279e+247 + 4.4528e+252i
                            0                            0

 Column 3:

                            0
                            0
   4.4600e-086 + 1.7804e-306i

which is wrong( but A*A gets right)

Anonymous

 

(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

Attach Files:
   
   
Comment:
   

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -email is unavailable- added by jordigh (Posted a comment)
  • -email is unavailable- added by melrobin
  • -email is unavailable- added by jwe (Posted a comment)
  • -email is unavailable- added by None (Submitted the item)
  •  

    There are 0 votes so far. Votes easily highlight which items people would like to see resolved in priority, independently of the priority of the item set by tracker managers.

    Only group members can vote.

     

    Follow 12 latest changes.

    Date Changed by Updated Field Previous Value => Replaced by
    2012-05-03 jordigh StatusIn Progress Fixed
        Open/ClosedOpen Closed
    2012-05-03 jordigh Summarymatrix operation incorrect complex diagonal matrix powers read uninit memory
    2012-05-03 jordigh StatusConfirmed In Progress
        Assigned toNone jordigh
    2012-05-02 melrobin Carbon-Copy- Added melrobin
    2012-05-02 jwe Severity3 - Normal 5 - Blocker
        Priority5 - Normal 7 - High
        Item GroupIncorrect Result Regression
        StatusNone Confirmed
        Release3.6.1 dev
        Operating SystemMicrosoft Windows Any

    Back to the top

    Powered by Savane 3.13-4448.
    Corresponding source code