octave-bug-tracker
[Top][All Lists]
Advanced

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

[Octave-bug-tracker] [bug #53440] Wrong factorization with ilutp type


From: Marco Caliari
Subject: [Octave-bug-tracker] [bug #53440] Wrong factorization with ilutp type
Date: Fri, 23 Mar 2018 09:31:49 -0400 (EDT)
User-agent: Mozilla/5.0 (X11; Ubuntu; Linux x86_64; rv:59.0) Gecko/20100101 Firefox/59.0

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

                 Summary: Wrong factorization with ilutp type
                 Project: GNU Octave
            Submitted by: caliari
            Submitted on: Fri 23 Mar 2018 01:31:47 PM UTC
                Category: Octave Function
                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: GNU/Linux

    _______________________________________________________

Details:

According to the documentation (and matlab behavior) the following code


A=sprand(10,10,1/2);opts.type="ilutp";[L,U]=ilu(A,opts);norm(L*U-A,inf)


should give something like eps. In fact, droptol is zero, which means complete
LU factorization. The explicit pivoted version


A=sprand(10,10,1/2);opts.type="ilutp";[L,U,P]=ilu(A,opts);norm(L*U-P*A,inf)


works fine.




    _______________________________________________________

Reply to this item at:

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

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




reply via email to

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