[Octave-bug-tracker] [bug #32035] qr.cc test error on octave-dev source (3.3.54) on MInGW

[Tatsuro MATSUOKA]

URL:
  <http://savannah.gnu.org/bugs/?32035>

                 Summary: qr.cc test error on octave-dev source (3.3.54) on
MInGW 
                 Project: GNU Octave
            Submitted by: tmacchant
            Submitted on: 2011年01月03日 08時26分40秒
                Category: None
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: Incorrect Result
                  Status: None
             Assigned to: None
         Originator Name: 
        Originator Email: 
             Open/Closed: Open
         Discussion Lock: Any
                 Release: dev
        Operating System: Microsoft Windows

    _______________________________________________________

Details:

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" ???

 







    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?32035>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/