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

Submitted by:  Tatsuro MATSUOKA <tmacchant>
Submitted on:  Mon 03 Jan 2011 08:26:40 AM UTC  
 
Category: NoneSeverity: 3 - Normal
Priority: 5 - NormalItem Group: Incorrect Result
Status: NoneAssigned to: None
Originator Name: Open/Closed: Closed
Release: devOperating System: Microsoft Windows

Add a New Comment(Rich Markup)
   

You are not logged in

Please log in, so followups can be emailed to you.

 

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>
Project 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) < 1e1eps)
assert(norm(vec(triu(R)-R),Inf) == 0)
assert(norm(vec(QR - [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(QR - [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 File(s):
   
   
Comment:
   

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -unavailable- added by rik5 (Posted a comment)
  • -unavailable- added by tmacchant (Submitted the item)
  •  

    Do you think this task is very important?
    If so, you can click here to add your encouragement to it.
    This task has 0 encouragements so far.

    Only project members can vote.

     

    Please enter the title of George Orwell's famous dystopian book (it's a date):

     

     

    Follows 1 latest change.

    Date Changed By Updated Field Previous Value => Replaced By
    Mon 21 Feb 2011 02:10:53 AM UTCrik5Open/ClosedOpen=>Closed

    Back to the top


    Powered by Savane 3.1-cleanup1