[ https://issues.apache.org/jira/browse/MATH-1417?page=com.atlassian.jira.plugin.system.issuetabpanels:comment-tabpanel&focusedCommentId=16002571#comment-16002571 ]
Gilles commented on MATH-1417: ------------------------------ bq. Is common math easily compiled once i checked it out from git ? Asumimg that you have "maven" installed, it should just be a matter of running this command: {noformat} mvn compile {noformat} > 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)