Return-Path: X-Original-To: archive-asf-public-internal@cust-asf2.ponee.io Delivered-To: archive-asf-public-internal@cust-asf2.ponee.io Received: from cust-asf.ponee.io (cust-asf.ponee.io [163.172.22.183]) by cust-asf2.ponee.io (Postfix) with ESMTP id 5D092200C73 for ; Wed, 10 May 2017 11:33:14 +0200 (CEST) Received: by cust-asf.ponee.io (Postfix) id 5ADA8160BB4; Wed, 10 May 2017 09:33:14 +0000 (UTC) Delivered-To: archive-asf-public@cust-asf.ponee.io Received: from mail.apache.org (hermes.apache.org [140.211.11.3]) by cust-asf.ponee.io (Postfix) with SMTP id 7AE25160BA8 for ; Wed, 10 May 2017 11:33:13 +0200 (CEST) Received: (qmail 86336 invoked by uid 500); 10 May 2017 09:33:12 -0000 Mailing-List: contact issues-help@commons.apache.org; run by ezmlm Precedence: bulk List-Help: List-Unsubscribe: List-Post: List-Id: Reply-To: issues@commons.apache.org Delivered-To: mailing list issues@commons.apache.org Received: (qmail 86325 invoked by uid 99); 10 May 2017 09:33:12 -0000 Received: from pnap-us-west-generic-nat.apache.org (HELO spamd3-us-west.apache.org) (209.188.14.142) by apache.org (qpsmtpd/0.29) with ESMTP; Wed, 10 May 2017 09:33:12 +0000 Received: from localhost (localhost [127.0.0.1]) by spamd3-us-west.apache.org (ASF Mail Server at spamd3-us-west.apache.org) with ESMTP id 3540418828C for ; Wed, 10 May 2017 09:33:12 +0000 (UTC) X-Virus-Scanned: Debian amavisd-new at spamd3-us-west.apache.org X-Spam-Flag: NO X-Spam-Score: -99.202 X-Spam-Level: X-Spam-Status: No, score=-99.202 tagged_above=-999 required=6.31 tests=[KAM_ASCII_DIVIDERS=0.8, RP_MATCHES_RCVD=-0.001, SPF_PASS=-0.001, USER_IN_WHITELIST=-100] autolearn=disabled Received: from mx1-lw-eu.apache.org ([10.40.0.8]) by localhost (spamd3-us-west.apache.org [10.40.0.10]) (amavisd-new, port 10024) with ESMTP id qz2lyLjI-E6d for ; Wed, 10 May 2017 09:33:10 +0000 (UTC) Received: from mailrelay1-us-west.apache.org (mailrelay1-us-west.apache.org [209.188.14.139]) by mx1-lw-eu.apache.org (ASF Mail Server at mx1-lw-eu.apache.org) with ESMTP id 6873D5F642 for ; Wed, 10 May 2017 09:33:10 +0000 (UTC) Received: from jira-lw-us.apache.org (unknown [207.244.88.139]) by mailrelay1-us-west.apache.org (ASF Mail Server at mailrelay1-us-west.apache.org) with ESMTP id 6F84CE002A for ; Wed, 10 May 2017 09:33:04 +0000 (UTC) Received: from jira-lw-us.apache.org (localhost [127.0.0.1]) by jira-lw-us.apache.org (ASF Mail Server at jira-lw-us.apache.org) with ESMTP id 2418F21E0A for ; Wed, 10 May 2017 09:33:04 +0000 (UTC) Date: Wed, 10 May 2017 09:33:04 +0000 (UTC) From: "Lecomte (JIRA)" To: issues@commons.apache.org Message-ID: In-Reply-To: References: Subject: [jira] [Updated] (MATH-1417) weird result in RRQR decomposition. MIME-Version: 1.0 Content-Type: text/plain; charset=utf-8 Content-Transfer-Encoding: quoted-printable X-JIRA-FingerPrint: 30527f35849b9dde25b450d4833f0394 archived-at: Wed, 10 May 2017 09:33:14 -0000 [ https://issues.apache.org/jira/browse/MATH-1417?page=3Dcom.atlassian= .jira.plugin.system.issuetabpanels:all-tabpanel ] Lecomte updated MATH-1417: -------------------------- Ok, it's ok for me. I use the current stable release 3.6.1 of apache math in an industrial prod= uct. How should I proceed to take advantage of this fix ? -----Message d'origine----- De=C2=A0: Gilles (JIRA) [mailto:jira@apache.org]=20 Envoy=C3=A9=C2=A0: mardi 9 mai 2017 21:36 =C3=80=C2=A0: LECOMTE Jean-Francois Objet=C2=A0: [jira] [Commented] (MATH-1417) weird result in RRQR decomposit= ion. [ https://issues.apache.org/jira/browse/MATH-1417?page=3Dcom.atlassian.= jira.plugin.system.issuetabpanels:comment-tabpanel&focusedCommentId=3D16003= 376#comment-16003376 ]=20 Gilles commented on MATH-1417: ------------------------------ Got it. :) Please check that commit ed1ce82d822ffe185875b7b7d38352f20171c096 fixed the= problem. -- This message was sent by Atlassian JIRA (v6.3.15#6346) __________________________ Avant d'imprimer, pensez =C3=A0 l'environnement ! Please consider the envir= onment before printing !=20 Ce message et toutes ses pi=C3=A8ces jointes sont confidentiels et =C3=A9ta= blis =C3=A0 l'intention exclusive de ses destinataires. Toute utilisation n= on conforme =C3=A0 sa destination, toute diffusion ou toute publication, to= tale ou partielle, est interdite, sauf autorisation expresse. IFP Energies = nouvelles d=C3=A9cline toute responsabilit=C3=A9 au titre de ce message. Th= is message and any attachments are confidential and intended solely for the= addressees. Any unauthorised use or dissemination is prohibited. IFP Energ= ies 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=20 > java 1.8 > Reporter: Lecomte > > When using RRQRDecomposition on rank deficient matrix, results are wrong. > double[][] Xi =3D { > {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:=20 > 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 =3D dense real , numRows =3D 9 , numCols =3D 8 > -1,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 =20 > 0,000 -1,000 0,000 0,000 -1,000 0,000 -1,000 0,000 =20 > 0,000 0,000 -1,000 0,000 0,000 0,000 0,000 0,000 =20 > 0,000 0,000 0,000 -1,000 0,000 0,000 0,000 0,000 =20 > 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 =20 > 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 =20 > 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 =20 > 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 =20 > 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 =20 > Type =3D dense real , numRows =3D 9 , numCols =3D 9 > 0,000 0,000 0,000 0,000 0,000 0,000 1,000 0,000 0,000 = =20 > -1,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 = =20 > 0,000 -1,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 = =20 > 0,000 0,000 -1,000 0,000 0,000 0,000 0,000 0,000 0,000 = =20 > 0,000 0,000 0,000 0,000 1,000 0,000 0,000 0,000 0,000 = =20 > 0,000 0,000 0,000 0,000 0,000 1,000 0,000 0,000 0,000 = =20 > 0,000 0,000 0,000 -1,000 0,000 0,000 0,000 0,000 0,000 = =20 > 0,000 0,000 0,000 0,000 0,000 0,000 0,000 1,000 0,000 = =20 > 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 1,000 = =20 > That are the results i expect. > =20 -- This message was sent by Atlassian JIRA (v6.3.15#6346)