# \#lapack

**URL:** https://discourse.julialang.org/tag/lapack/202.md

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

---

## [AOCL.jl v5.3.1](https://discourse.julialang.org/t/aocl-jl-v5-3-1/138718)

<div class="topic-metadata">

**Author:** [@agodbole](https://discourse.julialang.org/u/agodbole)\
**Replies:** 8\
**Last updated:** [August 14, 2026, 5:48am UTC](https://discourse.julialang.org/t/aocl-jl-v5-3-1/138718 "2026-08-14T05:48:05Z")

</div>

Hello all, AOCL is a set of numerical libraries optimized for AMD processors built on the AMD “Zen” architecture, spanning multiple generations. It supports AMD EPYC™, AMD Ryzen™, and AMD Ryzen™ Threadripper™ processor …

---

## [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…

---

## [Lapack wrappers and final arguments 1](https://discourse.julialang.org/t/lapack-wrappers-and-final-arguments-1/124391)

<div class="topic-metadata">

**Author:** [@juthohaegeman](https://discourse.julialang.org/u/juthohaegeman)\
**Replies:** 2\
**Last updated:** [January 3, 2025, 12:25am UTC](https://discourse.julialang.org/t/lapack-wrappers-and-final-arguments-1/124391 "2025-01-03T00:25:17Z")

</div>

Studying some of the Lapack wrappers over at LinearAlgebra.jl/src/lapack.jl at master · JuliaLang/LinearAlgebra.jl · GitHub, I am confused by the ccall structure. Most relevant lapack methods have info as their final arg…

---

## [\[ANN\] FastLapackInterface.jl v1.0.0: Non-allocating LAPACK factorizations](https://discourse.julialang.org/t/ann-fastlapackinterface-jl-v1-0-0-non-allocating-lapack-factorizations/83354)

<div class="topic-metadata">

**Author:** [@louisponet](https://discourse.julialang.org/u/louisponet)\
**Replies:** 14\
**Last updated:** [January 21, 2024, 6:58pm UTC](https://discourse.julialang.org/t/ann-fastlapackinterface-jl-v1-0-0-non-allocating-lapack-factorizations/83354 "2024-01-21T18:58:44Z")

</div>

The DynareJulia team and I are excited to announce the first major release of FastLapackInterface.jl. The package facilitates the pre-allocation of workspaces to be used with some of the most important LAPACK functions w…

---

## [Inplace version for LU decomposition in Julia](https://discourse.julialang.org/t/inplace-version-for-lu-decomposition-in-julia/108795)

<div class="topic-metadata">

**Author:** [@mrVeng](https://discourse.julialang.org/u/mrVeng)\
**Replies:** 8\
**Last updated:** [January 15, 2024, 2:55pm UTC](https://discourse.julialang.org/t/inplace-version-for-lu-decomposition-in-julia/108795 "2024-01-15T14:55:06Z")

</div>

Hi there! I have to compute the lu decomposition of a matrix that is dependent - I was thinking to create a buffer first and then inplace compute the lu decomposition when needed, but after googling and searching throug…

---

## [ERROR: LAPACKException(4)](https://discourse.julialang.org/t/error-lapackexception-4/104301)

<div class="topic-metadata">

**Author:** [@Sahil\_Khan](https://discourse.julialang.org/u/Sahil_Khan)\
**Replies:** 0\
**Last updated:** [September 27, 2023, 8:13am UTC](https://discourse.julialang.org/t/error-lapackexception-4/104301 "2023-09-27T08:13:36Z")

</div>

Can someone explain what this error is about? and what could be possible solution? I can’t share my model but it is differential equations based and have equations such du\[1\] = a\*(b)^n \* c\*(k)^m + a\*b. So, it contains e…

---

## [Multithreaded LAPACK function in a Threads.@threads loop](https://discourse.julialang.org/t/multithreaded-lapack-function-in-a-threads-threads-loop/45205)

<div class="topic-metadata">

**Author:** [@migarstka](https://discourse.julialang.org/u/migarstka)\
**Replies:** 2\
**Last updated:** [September 24, 2022, 6:42am UTC](https://discourse.julialang.org/t/multithreaded-lapack-function-in-a-threads-threads-loop/45205 "2022-09-24T06:42:31Z")

</div>

Hi everyone, Problem I am trying to understand how to optimally set the number of BLAS threads when I am using a multithreaded LAPACK function in a Threads.@threads loop. More specifically I want to compute the eigende…

---

## [Syevd is absent in LinearAlgebra.LAPACK](https://discourse.julialang.org/t/syevd-is-absent-in-linearalgebra-lapack/51714)

<div class="topic-metadata">

**Author:** [@Jae-Mo\_Lihm](https://discourse.julialang.org/u/Jae-Mo_Lihm)\
**Replies:** 0\
**Last updated:** [December 12, 2020, 12:48pm UTC](https://discourse.julialang.org/t/syevd-is-absent-in-linearalgebra-lapack/51714 "2020-12-12T12:48:35Z")

</div>

In the LinearAlgebra package, functions LAPACK.syev! and LAPACK.syevr! for solving the eigenvalue problem are present. But, LAPACK.syevd! is absent. As I understand, the ?syevd (and its Hermitian counterpart, zheevd) al…

---

## [ArgumentError: invalid argument #4 to LAPACK call](https://discourse.julialang.org/t/argumenterror-invalid-argument-4-to-lapack-call/46455)

<div class="topic-metadata">

**Author:** [@PhillipBC](https://discourse.julialang.org/u/PhillipBC)\
**Replies:** 7\
**Last updated:** [September 15, 2020, 3:51pm UTC](https://discourse.julialang.org/t/argumenterror-invalid-argument-4-to-lapack-call/46455 "2020-09-15T15:51:02Z")

</div>

This problem arose using Julia 1.4.2, and I then got version 1.5.1 and the problem is still there. If I create a matrix, then get its eigenvalues, then plot the eigenvalues, then get the eigenvalues again, I get the err…

---

## [LAPACK exception with gels!](https://discourse.julialang.org/t/lapack-exception-with-gels/43620)

<div class="topic-metadata">

**Author:** [@Joris\_Pinkse](https://discourse.julialang.org/u/Joris_Pinkse)\
**Replies:** 1\
**Last updated:** [July 28, 2020, 7:43pm UTC](https://discourse.julialang.org/t/lapack-exception-with-gels/43620 "2020-07-28T19:43:28Z")

</div>

If I run the following program then Julia quits with an exception. I suspect this has something to do with rounding in the factorization. The matrix B below can be replaced with a matrix of zeros of the same size and t…

---

## [OpenBLAS is faster than Intel MKL on AMD Hardware (Ryzen)](https://discourse.julialang.org/t/openblas-is-faster-than-intel-mkl-on-amd-hardware-ryzen/8033)

<div class="topic-metadata">

**Author:** [@Elrod](https://discourse.julialang.org/u/Elrod)\
**Replies:** 40\
**Last updated:** [June 19, 2020, 1:03am UTC](https://discourse.julialang.org/t/openblas-is-faster-than-intel-mkl-on-amd-hardware-ryzen/8033 "2020-06-19T01:03:05Z")

</div>

Questions about MKL vs OpenBLAS come up a lot, for example in comparisons with Matlab (linked to MKL), and a lot of users have issues building with MKL, eg here. Of course, one can easily download an MKL binary with Juli…

---

## [On solving sets of linear equations in Fortran vs Julia](https://discourse.julialang.org/t/on-solving-sets-of-linear-equations-in-fortran-vs-julia/33729)

<div class="topic-metadata">

**Author:** [@razvangheorghe](https://discourse.julialang.org/u/razvangheorghe)\
**Replies:** 14\
**Last updated:** [February 7, 2020, 12:23pm UTC](https://discourse.julialang.org/t/on-solving-sets-of-linear-equations-in-fortran-vs-julia/33729 "2020-02-07T12:23:27Z")

</div>

Hi, I’d like to translate a Fortran 90 program to Julia and then continue with my development in Julia. The bottleneck right now seem to be the solving of around 100 systems of linear equations where the matrix (~400x400…

---

## [Adding more lapack wrappers to base julia](https://discourse.julialang.org/t/adding-more-lapack-wrappers-to-base-julia/25490)

<div class="topic-metadata">

**Author:** [@gasagna](https://discourse.julialang.org/u/gasagna)\
**Replies:** 1\
**Last updated:** [June 23, 2019, 12:08pm UTC](https://discourse.julialang.org/t/adding-more-lapack-wrappers-to-base-julia/25490 "2019-06-23T12:08:27Z")

</div>

I find myself in need of few additional LAPACK routines that are not currently wrapped in base Julia. What is the policy for adding new wrappers, should I make a PR?

---

## [Eigenvalues of symmetric matrix in packed format (upper triangular entries stacked in vector)](https://discourse.julialang.org/t/eigenvalues-of-symmetric-matrix-in-packed-format-upper-triangular-entries-stacked-in-vector/18323)

<div class="topic-metadata">

**Author:** [@migarstka](https://discourse.julialang.org/u/migarstka)\
**Replies:** 0\
**Last updated:** [December 5, 2018, 2:46pm UTC](https://discourse.julialang.org/t/eigenvalues-of-symmetric-matrix-in-packed-format-upper-triangular-entries-stacked-in-vector/18323 "2018-12-05T14:46:20Z")

</div>

I am trying to compute the eigenvalues of a symmetric matrix in packed format, i.e. you stack the upper triangular entries in a vector (Packed Storage ) It seems like there exists a LAPACK function for that, e.g. spev. …

---

## [Multiprocessing in LAPACK is gone](https://discourse.julialang.org/t/multiprocessing-in-lapack-is-gone/10448)

<div class="topic-metadata">

**Author:** [@Lian\_Yunlong](https://discourse.julialang.org/u/Lian_Yunlong)\
**Replies:** 0\
**Last updated:** [April 20, 2018, 2:50am UTC](https://discourse.julialang.org/t/multiprocessing-in-lapack-is-gone/10448 "2018-04-20T02:50:07Z")

</div>

Hello everyone, I had built my own Julia excutable because I need to install CUDAdrv. But I was probably not aware of the multiprocessing options. Now when I call eigfact(), there is always only one core running, but on…