[
https://issues.apache.org/jira/browse/MATH-1417?page=com.atlassian.jira.plugin.system.issuetabpanels:comment-tabpanel&focusedCommentId=16002565#comment-16002565
]
Gilles edited comment on MATH-1417 at 5/9/17 12:04 PM:
-------------------------------------------------------
Hi Gilles,
Why not ?
Is common math easily compiled once i checked it out from git ?
Thanks.
Jef
was (Author: lecomtje):
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:j...@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)
