commons-issues mailing list archives

Site index · List index
Message view « Date » · « Thread »
Top « Date » · « Thread »
From "Gilles (JIRA)" <j...@apache.org>
Subject [jira] [Commented] (MATH-1417) weird result in RRQR decomposition.
Date Tue, 09 May 2017 14:42:04 GMT

    [ https://issues.apache.org/jira/browse/MATH-1417?page=com.atlassian.jira.plugin.system.issuetabpanels:comment-tabpanel&focusedCommentId=16002798#comment-16002798
] 

Gilles commented on MATH-1417:
------------------------------

Have you checked out the git "master" branch?
With your example (and no change to the source), I obtain 4 as the rank value.

> weird result in RRQR decomposition.
> -----------------------------------
>
>                 Key: MATH-1417
>                 URL: https://issues.apache.org/jira/browse/MATH-1417
>             Project: Commons Math
>          Issue Type: Bug
>    Affects Versions: 3.6.1
>         Environment: linux RH6
> netbeans 8.2 
> java 1.8
>            Reporter: Lecomte
>
> When using RRQRDecomposition on rank deficient matrix, results are wrong.
> double[][] Xi = {
>             {0, 0, 0, 0, 0, 0, 0, 0, 0},
>             {0, 1, 0, 0, 0, 0, 0, 0, 0},
>             {0, 0, 1, 0, 0, 0, 0, 0, 0},
>             {0, 0, 1, 0, 0, 0, 0, 0, 0},
>             {0, 0, 1, 0, 0, 0, 0, 0, 0},
>             {0, 0, 0, 1, 0, 0, 0, 0, 0},
>             {0, 0, 0, 0, 0, 0, 1, 0, 0},
>             {0, 0, 0, 0, 0, 0, 0, 0, 0},};
> With this matrix, i obtain: 
> rank 6
> R:
> |1,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000|
> |0,000 1,000 1,000 1,000 0,000 0,000 0,000 0,000|
> |0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000|
> |0,000 0,000 0,000 0,000 1,000 0,000 0,000 0,000|
> |0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000|
> |0,000 0,000 0,000 0,000 0,000 1,000 0,000 0,000|
> |0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000|
> |0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000|
> |0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000|
> Q:
> |0,000 0,000 1,000 0,000 0,000 0,000 0,000 0,000 0,000|
> |1,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000|
> |0,000 1,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000|
> |0,000 0,000 0,000 1,000 0,000 0,000 0,000 0,000 0,000|
> |0,000 0,000 0,000 0,000 1,000 0,000 0,000 0,000 0,000|
> |0,000 0,000 0,000 0,000 0,000 0,000 1,000 0,000 0,000|
> |0,000 0,000 0,000 0,000 0,000 1,000 0,000 0,000 0,000|
> |0,000 0,000 0,000 0,000 0,000 0,000 0,000 1,000 0,000|
> |0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 1,000|
> Where Scipy (lapack) or ejml gives me:
> rank 4
> Type = dense real , numRows = 9 , numCols = 8
> -1,000   0,000   0,000   0,000   0,000   0,000   0,000   0,000  
>  0,000  -1,000   0,000   0,000  -1,000   0,000  -1,000   0,000  
>  0,000   0,000  -1,000   0,000   0,000   0,000   0,000   0,000  
>  0,000   0,000   0,000  -1,000   0,000   0,000   0,000   0,000  
>  0,000   0,000   0,000   0,000   0,000   0,000   0,000   0,000  
>  0,000   0,000   0,000   0,000   0,000   0,000   0,000   0,000  
>  0,000   0,000   0,000   0,000   0,000   0,000   0,000   0,000  
>  0,000   0,000   0,000   0,000   0,000   0,000   0,000   0,000  
>  0,000   0,000   0,000   0,000   0,000   0,000   0,000   0,000  
> Type = dense real , numRows = 9 , numCols = 9
>  0,000   0,000   0,000   0,000   0,000   0,000   1,000   0,000   0,000  
> -1,000   0,000   0,000   0,000   0,000   0,000   0,000   0,000   0,000  
>  0,000  -1,000   0,000   0,000   0,000   0,000   0,000   0,000   0,000  
>  0,000   0,000  -1,000   0,000   0,000   0,000   0,000   0,000   0,000  
>  0,000   0,000   0,000   0,000   1,000   0,000   0,000   0,000   0,000  
>  0,000   0,000   0,000   0,000   0,000   1,000   0,000   0,000   0,000  
>  0,000   0,000   0,000  -1,000   0,000   0,000   0,000   0,000   0,000  
>  0,000   0,000   0,000   0,000   0,000   0,000   0,000   1,000   0,000  
>  0,000   0,000   0,000   0,000   0,000   0,000   0,000   0,000   1,000  
> That are the results i expect.
>  



--
This message was sent by Atlassian JIRA
(v6.3.15#6346)

Mime
View raw message