[Octave-bug-tracker] [bug #36612] rref incorrect output

[anonymous]

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

                 Summary: rref incorrect output
                 Project: GNU Octave
            Submitted by: None
            Submitted on: Sat 09 Jun 2012 07:34:25 AM UTC
                Category: None
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: Regression
                  Status: None
             Assigned to: None
         Originator Name: Will Kelly
        Originator Email: 
             Open/Closed: Open
         Discussion Lock: Any
                 Release: 3.4.3
        Operating System: GNU/Linux

    _______________________________________________________

Details:

For a given matrix X (see input.m) if one of the values is small enough
rref(X) will not scale this value to one. I haven't really had a chance to
compare Octave's rref script with my own. My script, however, does produce the
correct results. I've provided it so that if you wish, you can compare and
contrast the two. (Or use it if it meets your standards.) I only briefly took
a look at your script, long enough to copy the header information and the
usage script. I'm not sure your script utilized partial pivoting, but I
haven't really delved too deep into the matter. 
GNU Octave, version 3.4.3
Octave was configured for \"x86_64-redhat-linux-gnu\".
Fedora 16. Running Octave in QtOctave.


rref([1E-30,0;0,1;])
ans = 
      0 1
      0 0
Not equal to expected :
ans = 
      1 0
      0 1





    _______________________________________________________

File Attachments:


-------------------------------------------------------
Date: Sat 09 Jun 2012 07:34:25 AM UTC  Name: RREF.diff  Size: 8kB   By: None

<http://savannah.gnu.org/bugs/download.php?file_id=25998>
-------------------------------------------------------
Date: Sat 09 Jun 2012 07:34:25 AM UTC  Name: input.m  Size: 436B   By: None

<http://savannah.gnu.org/bugs/download.php?file_id=25999>

    _______________________________________________________

Reply to this item at:

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

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