[
https://issues.apache.org/jira/browse/MATH-1417?page=com.atlassian.jira.plugin.system.issuetabpanels:comment-tabpanel&focusedCommentId=16002933#comment-16002933
]
Gilles edited comment on MATH-1417 at 5/9/17 7:35 PM:
------------------------------------------------------
Sorry, but i am searching a nullspace and so i work with Xi transposed, not Xi
by itself. And with this transposed matrix I should also obtain rank 4.
I didn't test master branch but the code seems the same.
was (Author: lecomtje):
Sorry, but i am searching a nullspace and so i work with Xi transposed, not Xi
by itself. And with this transposed matrix I should also obtain rank 4.
I didn't test master branch but the code seems the same.
-----Message d'origine-----
De : Gilles (JIRA) [mailto:j...@apache.org]
Envoyé : mardi 9 mai 2017 16:42
À : 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=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.
--
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)
