Skip to content

Add underdetermined (m < n) LQ support to zgels_gpu#89

Open
jgillis wants to merge 1 commit into
icl-utk-edu:masterfrom
jgillis:zgels-lq
Open

Add underdetermined (m < n) LQ support to zgels_gpu#89
jgillis wants to merge 1 commit into
icl-utk-edu:masterfrom
jgillis:zgels-lq

Conversation

@jgillis

@jgillis jgillis commented Apr 8, 2026

Copy link
Copy Markdown

Summary

  • Extends zgels_gpu to handle underdetermined systems (m < n) via LQ factorization
  • Adds zgelqs_gpu (LQ solve) and zunmlq_gpu (apply Q from LQ), adapted from their QR counterparts
  • Updates test harness to cover the m < n case

Test plan

  • testing_dgels_gpu passes for default (square/overdetermined) sizes
  • Numerical residuals within machine epsilon tolerance

🤖 Generated with Claude Code

zgels_gpu previously only handled overdetermined (m >= n) systems via QR.
This adds the LQ factorization path for underdetermined systems, computing
the minimum-norm solution via gelqf_gpu + gelqs_gpu.

Co-Authored-By: Claude Opus 4.6 (1M context) <noreply@anthropic.com>
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant