# \#svd

**URL:** https://discourse.julialang.org/tag/svd/1123.md

[Latest](https://discourse.julialang.org/latest.md) · [Categories](https://discourse.julialang.org/categories.md) · [Tags](https://discourse.julialang.org/tags.md)

---

## [AOCL (not MKL) acceleration on AMD Ryzen CPU's](https://discourse.julialang.org/t/aocl-not-mkl-acceleration-on-amd-ryzen-cpus/112890)

<div class="topic-metadata">

**Author:** [@aasdelat](https://discourse.julialang.org/u/aasdelat)\
**Replies:** 28\
**Last updated:** [August 10, 2026, 10:54am UTC](https://discourse.julialang.org/t/aocl-not-mkl-acceleration-on-amd-ryzen-cpus/112890 "2026-08-10T10:54:52Z")

</div>

Hi: I am interested in calling a lapack routine, dgesvd, that makes a singular value decomposition of a matrix. The library LinearAlgebra.jl already does it, but it is a generic library that does not take full advantag…

---

## [Truncation error in ITensors SVD](https://discourse.julialang.org/t/truncation-error-in-itensors-svd/134559)

<div class="topic-metadata">

**Author:** [@diffeqslvr](https://discourse.julialang.org/u/diffeqslvr)\
**Replies:** 2\
**Last updated:** [December 19, 2025, 10:14pm UTC](https://discourse.julialang.org/t/truncation-error-in-itensors-svd/134559 "2025-12-19T22:14:05Z")

</div>

I am trying to understand the truncation error when performing an SVD using ITensors. According to ITensor Examples · ITensors.jl, the truncation error is defined as sum\_{discarded n} lambda\_n^2 / sum\_{n} lambda\_n^2, whe…

---

## [Which package is suitable for matrix-free partial SVD?](https://discourse.julialang.org/t/which-package-is-suitable-for-matrix-free-partial-svd/118419)

<div class="topic-metadata">

**Author:** [@MAChiyu](https://discourse.julialang.org/u/MAChiyu)\
**Replies:** 2\
**Last updated:** [August 20, 2024, 8:15pm UTC](https://discourse.julialang.org/t/which-package-is-suitable-for-matrix-free-partial-svd/118419 "2024-08-20T20:15:57Z")

</div>

Hi! I am new to Julia. I used to use MATLAB. Within MATLAB, I use lansvd in PROPACK to do partial SVD of some linear operator (needs Afunpar to input other parameters of the matrix free operator). I notice there is PROPA…

---

## [Warning/Error handling from complex, sparse array on SVD (PROPACK)](https://discourse.julialang.org/t/warning-error-handling-from-complex-sparse-array-on-svd-propack/110779)

<div class="topic-metadata">

**Author:** [@KSfairopoulos](https://discourse.julialang.org/u/KSfairopoulos)\
**Replies:** 0\
**Last updated:** [February 26, 2024, 2:16pm UTC](https://discourse.julialang.org/t/warning-error-handling-from-complex-sparse-array-on-svd-propack/110779 "2024-02-26T14:16:45Z")

</div>

Let’s assume I have a sparse matrix with complex entries. I want to perform svd on it. So, I use the Propack package (do you have any better suggestions?). I do this with the command: tsvdvals(gs, k=p), where I retain …

---

## [Incorporating SVD factorization into GPUs with CUDA](https://discourse.julialang.org/t/incorporating-svd-factorization-into-gpus-with-cuda/99252)

<div class="topic-metadata">

**Author:** [@jarroyoe](https://discourse.julialang.org/u/jarroyoe)\
**Replies:** 7\
**Last updated:** [May 23, 2023, 2:31am UTC](https://discourse.julialang.org/t/incorporating-svd-factorization-into-gpus-with-cuda/99252 "2023-05-23T02:31:21Z")

</div>

I’m trying to run the following code to optimize the spectral radius of the weight matrices of a neural network: using CUDA, Random, Lux, ComponentArrays, Optimization, OptimizationOptimisers, OptimizationOptimJL, Linea…

---

## [Compressed svd algorithm implementation](https://discourse.julialang.org/t/compressed-svd-algorithm-implementation/97837)

<div class="topic-metadata">

**Author:** [@Desmond](https://discourse.julialang.org/u/Desmond)\
**Replies:** 4\
**Last updated:** [April 25, 2023, 11:34pm UTC](https://discourse.julialang.org/t/compressed-svd-algorithm-implementation/97837 "2023-04-25T23:34:32Z")

</div>

I am trying to implement the compressed SVD in this \[https://openaccess.thecvf.com/content\_ICCV\_2017\_workshops/papers/w25/Erichson\_Compressed\_Singular\_Value\_ICCV\_2017\_paper.pdf\]. Attached is the algorithm. I have proble…

---

## [NaNs In The Output of SVD](https://discourse.julialang.org/t/nans-in-the-output-of-svd/93712)

<div class="topic-metadata">

**Author:** [@robin-armstrong](https://discourse.julialang.org/u/robin-armstrong)\
**Replies:** 1\
**Last updated:** [January 28, 2023, 10:22pm UTC](https://discourse.julialang.org/t/nans-in-the-output-of-svd/93712 "2023-01-28T22:22:47Z")

</div>

I am currently running numerical experiments which involve computing the SVD of a very large number of random matrices. Occasionally, Julia’s SVD function will return a Vt factor that contains a row of all NaNs. Here is …

---

## [SVD on a matrix that does not fit into memory](https://discourse.julialang.org/t/svd-on-a-matrix-that-does-not-fit-into-memory/92802)

<div class="topic-metadata">

**Author:** [@aasdelat](https://discourse.julialang.org/u/aasdelat)\
**Replies:** 20\
**Last updated:** [January 16, 2023, 3:44pm UTC](https://discourse.julialang.org/t/svd-on-a-matrix-that-does-not-fit-into-memory/92802 "2023-01-16T15:44:17Z")

</div>

I am looking for a function that makes a SVD on a matrix that does not fit into memory. I am interested in using just one computer, although an answer for a cluster can be useful for other users. The python dask package…

---

## [Does the singular value degeneracy cause the gradient to be NAN？](https://discourse.julialang.org/t/does-the-singular-value-degeneracy-cause-the-gradient-to-be-nan/86412)

<div class="topic-metadata">

**Author:** [@F-YF](https://discourse.julialang.org/u/F-YF)\
**Replies:** 11\
**Last updated:** [August 29, 2022, 5:37pm UTC](https://discourse.julialang.org/t/does-the-singular-value-degeneracy-cause-the-gradient-to-be-nan/86412 "2022-08-29T17:37:07Z")

</div>

I recently encountered a situation where the result of a computation was of type NAN. I got the following example after a step-by-step debug: using Zygote using OMEinsum using LinearAlgebra function f(T) @ein M\[ b…