[Octave-bug-tracker] [bug #49904] Inversion of a Hermitian matrix produces incorrect result

[Jonathan Lai]

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

                 Summary: Inversion of a Hermitian matrix produces incorrect
result
                 Project: GNU Octave
            Submitted by: jlai319
            Submitted on: Wed 21 Dec 2016 04:06:58 PM GMT
                Category: Libraries
                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: 4.2.0
        Operating System: Microsoft Windows

    _______________________________________________________

Details:

>> Rxx
Rxx =

   1.00003 + 0.00000i  -0.68292 + 0.26406i  -0.57670 + 0.25291i   0.34986 +
0.03335i
  -0.68292 - 0.26406i   0.94539 + 0.00000i   0.69006 - 0.03967i  -0.43847 +
0.38619i
  -0.57670 - 0.25291i   0.69006 + 0.03967i   0.95260 + 0.00000i  -0.34734 +
0.40728i
   0.34986 - 0.03335i  -0.43847 - 0.38619i  -0.34734 - 0.40728i   0.98356 +
0.00000i

>> inv(Rxx)
ans =

   7412.81 +    0.00i   4404.95 - 4826.59i   3353.84 - 1286.54i  -1916.52 -
5975.72i
   4404.95 - 4826.59i   5762.91 +    0.00i   2829.13 + 1419.57i   2752.81 -
4799.16i
   3353.84 - 1286.54i   2829.13 + 1419.57i   1743.11 +    0.00i    170.06 -
3036.83i
  -1916.52 - 5975.72i   2752.81 - 4799.16i    170.06 - 3036.83i   5314.48 +   
0.00i

Comment: Rxx is a Hermitian, but inv(Rxx) is not





    _______________________________________________________

Reply to this item at:

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

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