標簽[qr-decomposition] - 堆棧內存溢出

[英]QR factorization in R not giving correct answer with large sparse matrix

我有一個 dgTMatrix 類型的稀疏矩陣，它是用 R 中的矩陣dgTMatrix構建的，我正在嘗試使用 QR 分解求解Ax=b中的向量x ，但它無法正常工作。例如下面我正在求解隨機A, b ，您可以看到該方法適用於Ax-b = 0的范數我A和b是這樣的：當我對A進行 QR 分解時，我得到以 ...

如何獲得 pivot 並像 base::qr() 一樣從 Matrix::qr() 中排名？

[英]How to get the pivot and rank from Matrix::qr() like that of base::qr()?

在 R 中的稀疏矩陣上應用 Matrix::qr() 時，output 與 base::qr 有很大不同。有V、beta、p、R、q但沒有rank和pivot。下面是一個小示例代碼。我想檢測 A 稀疏矩陣的線性相關列，這需要 pivot 和秩。我應該如何獲得這些信息？還有另一個相關但未回 ...

使用 QR 分解的 Python 線性回歸（相關特征）

[英]Python Linear regression using QR decomposition (correlated features)

我需要在某些特征相關的數據集上應用 Python 線性回歸。 Scikit 的線性回歸使用奇異值分解來最小化平方誤差。但是，這種方法沒有考慮變量之間的相關性。 R 使用 QR 分解，這是一種使用“lm”方法對協變量進行去相關的技術。相關變量的系數作為 NaN 值返回，這對於我想要做的事情很 ...

R中的快速QR分解

[英]Fast QR Factorization in R

我有大量矩陣，我需要對其執行 QR 分解並存儲生成的 Q 矩陣（歸一化以使 R 矩陣在其對角線上具有正數）。除了使用qr()函數之外，還有其他方法嗎？這是工作示例：所以計算 62,500 QR 分解需要 3.07 秒。我想知道是否有更快的東西？ ...

對 `qr()` 的信心動搖了

[英]Shaken faith in `qr()`

我在處理排名不足的情況時非常依賴qr()函數，但最近遇到了一些無法正常工作的示例。考慮下面的矩陣badX ： badX <- structure(c(-1.641906809157e-10, 0, 0, 0, 0, -0.5, 0, 0, -1.10482935525559e-16, ...

為什么 QR 分解不能正常工作？（Lapacke，復雜案例）

[英]Why doesn't QR decomposition work correctly? (Lapacke, complex case)

我用拉帕克。我正在嘗試在 C 中對復雜數據進行 QR 分解。為此，我編寫了 function （基於 Haatschii 代碼How to get the Q from the QR factorization output? ）：值得注意的是，Alireza 的作者也遇到了 functio ...

特征向量的 QR 方法 Python

[英]QR method for eigenvectors Python

我正在嘗試使用 QR 方法找到矩陣 A 的特征向量。我找到了對應於最大特征值的特征值和特征向量。如何在不使用 numpy.linalg.eig 的情況下找到特征向量的 rest？接下來我檢查條件：檢查條件： ...

Xilinx Vivado_HLS 提供的線性代數庫的 QRF_ALT 函數中使用的算法是什么？

[英]What is the algorithm used in QRF_ALT function in linear Algebra library provided with the Xilinx Vivado_HLS?

隨着 Vivado HLS 安裝，提供了一個線性代數庫。在這個庫中有這個函數 QRF_ALT，它應該是 Qrf-basic 的高吞吐量版本：使用給定旋轉的 QR 分解。我的問題是它如何實現高吞吐量？從代碼附帶的簡短評論中，我的猜測是它使用了批量矩陣行的並行處理，但是它是如何調度它的，我無法理 ...

numpy 中的 QuiRk - 所有歸一化元素在使用反射器的 QR 分解中分配為 0

[英]QuiRk in numpy - All normalized elements get assigned to 0 in QR decomposition using reflectors

數組中的所有標准化元素都設置為 0。我正在嘗試為陣列創建反射器。當我通過 Python 中的第一個元素(x[1:] /= x[0], x[0] = 1)對數組中的所有元素進行歸一化時，我只會在第一個元素中得到全零和 1。結果應如下所示(u = [1, 4/9, 1/9, 3/9, 1/9]) ...

用於 QR 分解的戶主反射

[英]House-Holder Reflection for QR Decomposition

我正在嘗試通過家庭反射器實現 QR 分解。在一個非常簡單的數組上嘗試這個時，我得到了奇怪的數字。任何人都可以告訴我，為什么在 function 定義的最后一行在 vec 和 vec.T 之間使用 @ vs * 運算符可以獲得主要獎勵積分。截至今天早上，這已經難倒了兩位數學/計算機科學博士。 ...

使用Numpy優化正交於矢量的矩陣投影

[英]Optimized projection of a matrix orthogonaly to a vector with Numpy

我需要使矩陣A所有其他列正交於其列j 。我使用以下算法：這個想法來自用Gram-Schmidt過程進行QR分解。但是由於Gram-Schmidt過程，該代碼並未經過優化且不穩定。 Numpy是否提供任何方法來計算這些向量的正交投影？帶戶主矩陣我 ...

在 Python 中使用 QR 分解解決超定系統

[英]Solve overdetermined system with QR decomposition in Python

我試圖用 QR 分解和 linalg.solve 解決一個超定系統，但我得到的錯誤是 LinAlgError: 數組的最后 2 個維度必須是 square 。當 R 數組不是正方形時會發生這種情況，對嗎？代碼看起來像這樣有沒有辦法以更有效的方式為任意 A 維度編寫此代碼？如果沒有，我 ...

如何在 Python 中實現基於 Householder 的 QR 分解？

[英]How can you implement Householder based QR decomposition in Python?

我目前正在嘗試為矩形矩陣實現基於 Householder 的 QR 分解，如http://eprints.ma.man.ac.uk/1192/1/qrupdating_12nov08.pdf （第 3、4、5 頁）中所述。顯然我弄錯了一些偽代碼，因為 (1) 我的結果不同於numpy.qr.li ...

QR分解Fortran錯誤

[英]QR decomposition Fortran error

我對QR分解方法有疑問。我使用dgeqrf子例程進行分解，但編譯器中沒有錯誤，但此后出現了問題。我還沒發現錯誤在哪里。另一個問題是，如果A矩陣為零，則A = Q * R =>可以分解為零或失去秩。 ...

MATLAB中具有給定旋轉的GMRES方法

[英]GMRES method with given rotations in MATLAB

我有以下算法實現我的問題是在每次迭代（可能是每次重新啟動）時獲得一個向量（可能是矩陣） vetnorm_r ，其中包含殘差的所有范數（作為輸出）。我不知道如何建立這個向量或矩陣。％輸入A REAL非對稱正定矩陣％x REAL初始猜測向量％b REAL右側向量％M REAL ...

cuSolverDN或其他CUDA庫是否具有針對稠密矩陣的QR分解的批處理版本以解決A * x = b？

[英]Does cuSolverDN or another CUDA library have a batched-version of QR decomposition for dense matrices to solve A*x = b?

我正在嘗試求解A * x = b，其中A具有復雜的值並且很密集。我使用了cuSolverDN庫中的cusolverDnCgeqrf（）方法對一組線性方程組進行QR分解。但是，我想多次執行此操作以加快處理速度。是否有此方法的“批處理”版本？還是我可以使用另一個CUDA庫？ ...

“輸入意外結束”錯誤是什么意思？

[英]What does the “Unexpected end of input” error mean?

我目前正在R中構建一個應用程序，以計算QR矩陣分解，QR非負矩陣分解和計算ICA。目前，我正在完成第一個任務。我收到以下錯誤： source（“ trial.R”）source（“ trial.R”）中的錯誤：trial.R：153：0：輸入151的意外結束：152：^ ...

張量分解能做什么？

[英]What can we do with a tensor factorization?

我有三個關於張量因子分解的問題。張量因子分解（分解）的情況（或應用）是什么？這有可能成為未來的主流技術？你如何使用它？ ...

如何在Julia中最有效地使用QR分解？

[英]How to most efficiently use QR-decomposition in Julia?

避免陣列分配有利於提高性能。但是，我還沒有理解什么是可以執行矩陣A的QR分解的最有效方式。（注意：需要Q和R矩陣）簡單地使用Q, R = qr(A)可能不是最好的想法，因為它分配Q和R兩者，其中兩者都可以重新分配。函數qrfact允許以打包格式存儲分解。但是，我之后 ...

通過Householder成功找到Q＆R之后，無法獲得R中的線性回歸系數

[英]Unable to get Linear Regression Cofficients in R after Successfully finding Q & R via Householder

我正在手動嘗試計算回歸系數，而不是對數據使用任何默認值http://people.sc.fsu.edu/~jburkardt/datasets/regression/x31.txt 這是我的代碼，可以正確產生滿足A = QR的Q＆R。但是我找不到系數作為問與答產生問題的維度。有什么 ...