This bugzilla service is closed. All entries have been migrated to https://gitlab.com/libeigen/eigen
Bug 941 - Lacking accuracy in ColPivHouseholderQR
Summary: Lacking accuracy in ColPivHouseholderQR
Status: RESOLVED FIXED
Alias: None
Product: Eigen
Classification: Unclassified
Component: QR (show other bugs)
Version: 3.2
Hardware: All All
: Normal Accuracy Problem
Assignee: Nobody
URL:
Whiteboard:
Keywords:
Depends on:
Blocks:
 
Reported: 2015-01-29 14:47 UTC by ettersi
Modified: 2019-12-04 14:08 UTC (History)
3 users (show)



Attachments
Matrix causing the problem (236.25 KB, text/plain)
2015-01-29 14:47 UTC, ettersi
no flags Details
MWE (1.15 KB, text/x-c++src)
2015-01-29 14:48 UTC, ettersi
no flags Details

Description ettersi 2015-01-29 14:47:57 UTC
Created attachment 532 [details]
Matrix causing the problem

For the attached matrix, ColPivHouseholderQR yields a decomposition error 

(A*P - Q*R).norm()

of 4.28e-11 (using Scalar = double), i.e. it loses 5 orders of magnitude. 

For comparison, non-pivoting HouseholderQR achieves a decomposition error of 1.98e-16, and also the pivoting QR implementation of Matlab achieves a decomposition error of 4.11e-16.
Comment 1 ettersi 2015-01-29 14:48:35 UTC
Created attachment 533 [details]
MWE
Comment 2 Gael Guennebaud 2015-01-30 19:05:21 UTC
Thank you for the detailed report.

https://bitbucket.org/eigen/eigen/commits/f3ecd92ab19f/
Changeset:   f3ecd92ab19f
User:        ggael
Date:        2015-01-30 18:04:04+00:00
Summary:     Bug 941: fix accuracy issue in ColPivHouseholderQR, do not stop decomposition on a small pivot

3.2:
https://bitbucket.org/eigen/eigen/commits/cb8ac3b2bbd6/
Comment 3 Nobody 2019-12-04 14:08:52 UTC
-- GitLab Migration Automatic Message --

This bug has been migrated to gitlab.com's GitLab instance and has been closed from further activity.

You can subscribe and participate further through the new bug through this link to our GitLab instance: https://gitlab.com/libeigen/eigen/issues/941.

Note You need to log in before you can comment on or make changes to this bug.