bugGNU Octave - Bugs: bug #32035, qr.cc test error on octave-dev...

 
 

bug #32035: qr.cc test error on octave-dev source (3.3.54) on MInGW

Submitter:  Tatsuro MATSUOKA <tmacchant>
Submitted:  Mon 03 Jan 2011 08:26:40 AM UTC
   
 
Category:  None Severity:  3 - Normal
Priority:  5 - Normal Item Group:  Incorrect Result
Status:  None Assigned to:  None
Originator Name:  Open/Closed:  * Closed
Release:  * dev Operating System:  * Microsoft Windows
Fixed Release:  None Planned Release:  None
* Mandatory Fields

Add a New Comment Rich Markup
   

Mon 21 Feb 2011 07:44:18 AM UTC, comment #3: 

Currently this error does not occur.  Therefore this post should be closed.

Tatsuro MATSUOKA <tmacchant>
Mon 21 Feb 2011 02:10:53 AM UTC, comment #2: 

Closed bug, per your request as this no longer seems to be an issue.

Rik <rik5>
Group administrator
Sat 08 Jan 2011 06:19:44 PM UTC, comment #1: 

I have back the blas runtime from GotoBLAS to the ATLAS and Lapack is updated to 3.3.0.  The error reported does not appear so that I would like to close this item

Tatsuro MATSUOKA <tmacchant>
Mon 03 Jan 2011 08:26:40 AM UTC, original submission:  

Hello

This seem not to be reported.
In the fntests

c:\usr\Tatsu\mingwhome\octaves\hg\octave-work\src/DLD-FUNCTIONS\qr.cc  PASS47/48   FAIL 1

From fntests.log

  *** test
 [Q,R] = qr(Ac);
 [Q,R] = qrinsert(Q,R,3,uc);
 assert(norm(vec(Q'*Q - eye(5)),Inf) < 1e1*eps)
 assert(norm(vec(triu(R)-R),Inf) == 0)
 assert(norm(vec(Q*R - [Ac(:,1:2) uc Ac(:,3)]),Inf) < norm(Ac)*1e1*eps)
!!!!! test failed
assert (norm (vec (Q R - [Ac(:, 1:2), uc, Ac(:, 3)]), Inf) < norm (Ac) 1e1 * eps) failed
shared variables {
  A =

     0.091364   0.613038   0.999083
     0.594638   0.425302   0.603537
     0.383594   0.291238   0.085574
     0.265712   0.268003   0.238409
     0.669966   0.743851   0.445057

  u =

     0.85082
     0.76426
     0.42883
     0.53010
     0.80683

  v =

     0.98810
     0.24295
     0.43167

  Ac =

     0.620405 + 0.956953i   0.480013 + 0.048806i   0.402627 + 0.338171i
     0.589077 + 0.658457i   0.013205 + 0.279323i   0.229284 + 0.721929i
     0.092758 + 0.345687i   0.928679 + 0.241052i   0.764536 + 0.832406i
     0.912098 + 0.721024i   0.049018 + 0.269452i   0.730029 + 0.796517i
     0.112849 + 0.603871i   0.486352 + 0.142337i   0.355646 + 0.151496i

  uc =

     0.20351 + 0.05401i
     0.13141 + 0.43708i
     0.29808 + 0.08789i
     0.69821 + 0.38844i
     0.74871 + 0.25821i

  vc =

     0.85839 + 0.29468i
     0.20820 + 0.93090i
     0.86184 + 0.34689i
#**********************
I have tried the following , after Ac and uc are set properly.

>> [Q,R] = qr(Ac);
>> [Q,R] = qrinsert(Q,R,3,uc);
>>  norm(vec(Q*R - [Ac(:,1:2) uc Ac(:,3)]),Inf)

ans = 4.7103e-016

>> norm (Ac) 1e1 eps

ans = 6.0510e-015

>> norm (vec (Q R - [Ac(:, 1:2), uc, Ac(:, 3)]), Inf) < norm (Ac) 1e1 * eps

ans =  1

>> assert(norm(vec(Q*R - [Ac(:,1:2) uc Ac(:,3)]),Inf) < norm(Ac)*1e1*eps)
>>

???????
Why fntests did report
"assert (norm (vec (Q R - [Ac(:, 1:2), uc, Ac(:, 3)]), Inf) < norm (Ac) 1e1 * eps) failed" ???

 



Tatsuro MATSUOKA <tmacchant>

 

(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 rik5 (Posted a comment)
  • -email is unavailable- added by tmacchant (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.

     

    Follows 1 latest change.

    Date Changed by Updated Field Previous Value => Replaced by
    2011-02-21 rik5 Open/ClosedOpen Closed

    Back to the top

    Powered by Savane 3.13-d3ae.
    Corresponding source code