help-octave
[Top][All Lists]
Advanced

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

Re: Octave criterion for matrix singularity


From: Carlo De Falco
Subject: Re: Octave criterion for matrix singularity
Date: Tue, 26 Feb 2019 13:00:25 +0000


> On 25 Feb 2019, at 19:34, mabalenk <address@hidden> wrote:
> 
> Dear all,
> 
> I'm repeatedly solving a linear system with a Vandermonde matrix. I'm using
> a "backslash" operator:
> 
>  y = -W'/b'
> 
> At each iteration I'm increasing the matrix order by one. At some point
> Octave issues a warning:
> 
>  warning: matrix singular to machine precision, rcond = 1.56508e-17
> 
> If possible, would you please direct me to the Octave source code, where
> this decision is made? I would like to understand the criterion, that
> triggers the warning. I found a post on StackOverflow that says MATLAB would
> issue a similar warning, once rcond < 1e-12. I would like to know, what
> threshold does Octave use and how is it defined. Thank you!
> 
> --
> Best wishes,
> Maxim



The code for solving A\b with A real double precision, full matrix is here :

 
http://hg.savannah.gnu.org/hgweb/octave/file/14815cb62fbe/liboctave/array/dMatrix.cc#l1445

the condition being checked is

 rcond + 1.0 == 1.0

which means rcond is less or equal than eps(1.0)/2

HTH,
c.






reply via email to

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