Skip to content
GitLab
Projects Groups Snippets
  • /
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Sign in / Register
  • openfoam openfoam
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 415
    • Issues 415
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 7
    • Merge requests 7
  • Deployments
    • Deployments
    • Releases
  • Wiki
    • Wiki
  • Activity
  • Graph
  • Create a new issue
  • Commits
  • Issue Boards
Collapse sidebar
  • Development
  • openfoamopenfoam
  • Issues
  • #2212
Closed
Open
Issue created Sep 17, 2021 by Kutalmış Berçin@kutiMaintainer

BUG: QRMatrix: Prepended zero-valued rows lead to wrong QR decompositions

Summary

Matrices prepended with zero-valued rows lead to wrong QR decompositions. (Appended zero rows are fine). e.g.:

A = 
4 3
(
  (0 0 0)
  (1 2 3)
  (5 6 7)
  (9 10 11)
)

Reduced QR decomposition produces the following wrong R matrix:

R =
3 3
(
  (0 0 0)
  (-1 11.8322 13.3534)
  (-5 -8.88178e-16 0.828079)
)

rather than the correct R matrix:

R =
[[-1.03440804e+01 -1.17941852e+01 -1.32442899e+01]
 [ 0.00000000e+00 -9.47204446e-01 -1.89440889e+00]
 [ 0.00000000e+00  0.00000000e+00  7.10058168e-15]]
)

Steps to reproduce

Compile and execute: Test-A.zip

Environment information

    api      = 2107
    patch    = 0
    HEAD     = 72e67c7d1b
    version  = com
    compiler = Clang (system)
             = clang version 9.0.1 
    mpi      = SYSTEMOPENMPI
             = mpirun (Open MPI) 1.10.7.0.5e373bf1fd
    OS       = Description:     openSUSE Leap 15.3
    opts     = linux64ClangDPInt32Opt
Edited Sep 17, 2021 by Kutalmış Berçin
Assignee
Assign to
Time tracking