commons-issues mailing list archives

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

     [ https://issues.apache.org/jira/browse/MATH-1417?page=com.atlassian.jira.plugin.system.issuetabpanels:all-tabpanel
]

Lecomte updated MATH-1417:
--------------------------

Hi Gilles,

Why not ? 
Is common math easily compiled once i checked it out  from git ?

Thanks.

Jef

-----Message d'origine-----
De : Gilles (JIRA) [mailto:jira@apache.org] 
Envoyé : mardi 9 mai 2017 12:52
À : LECOMTE Jean-Francois
Objet : [jira] [Commented] (MATH-1417) weird result in RRQR decomposition.


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

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

Thanks for the report!

Are you willing to dig into the code in order to locate the bug?




--
This message was sent by Atlassian JIRA
(v6.3.15#6346)
__________________________
Avant d'imprimer, pensez à l'environnement ! Please consider the environment before printing
! 
Ce message et toutes ses pièces jointes sont confidentiels et établis à l'intention exclusive
de ses destinataires. Toute utilisation non conforme à sa destination, toute diffusion ou
toute publication, totale ou partielle, est interdite, sauf autorisation expresse. IFP Energies
nouvelles décline toute responsabilité au titre de ce message. This message and any attachments
are confidential and intended solely for the addressees. Any unauthorised use or dissemination
is prohibited. IFP Energies nouvelles should not be liable for this message.
__________________________


> 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