[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Octave-bug-tracker] [bug #62867] wrong result when solving large-scale

From: anonymous
Subject: [Octave-bug-tracker] [bug #62867] wrong result when solving large-scale linear equation
Date: Fri, 5 Aug 2022 15:48:15 -0400 (EDT)


                 Summary: wrong result when solving large-scale linear
                 Project: GNU Octave
               Submitter: None
               Submitted: Fri 05 Aug 2022 07:48:14 PM UTC
                Category: Octave Function
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: Matlab Compatibility
                  Status: None
             Assigned to: None
         Originator Name: Junjie Yu
        Originator Email: junjieyu@usc.edu
             Open/Closed: Open
                 Release: 7.2.0
         Discussion Lock: Any
        Operating System: Any


Follow-up Comments:

Date: Fri 05 Aug 2022 07:48:14 PM UTC By: Anonymous
I was working on running some Matlab files in Octave (7.2.0) and got an
incorrect answer. After digging for a while, I found Octave meets 'matrix
singular to machine precision' warning when solving large-size linear
equations. Below I provide a simple test case for reproducing the error.
We want to solve a linear system: A*x = rhs(see test.m file)
Where A has dimension 93248*93248 and is very sparse, rhs has the dimension of
we can solve x using x = mldivide(A, rhs);
matlab can solve this correctly, octave lead to precision error and give a
totally wrong result.

warning: matrix singular to machine precision, rcond = 1.31465e-17
warning: called from
    test at line 6 column 3

ans = -0.027239

ans =


Could you check whether you can reproduce this error and give me any
suggestions to handle this?
Thanks so much!


Reply to this item at:


Message sent via Savannah

reply via email to

[Prev in Thread] Current Thread [Next in Thread]