lapack 相关问题

LAPACK(Linear Algebra PACKage)是一个用于求解线性代数方程的软件库包。 LAPACK是用Fortran 90编写的。

为什么scalapack中没有GEEV的等价物?

我希望使用 ScaLAPACK 找到大型复杂非埃尔米特矩阵的特征值。我似乎找不到这样一个通用特征求解器的任何文档(大概名称是 so...

回答 1 投票 0

从 lapack 库调用 F77_NAME() 函数时我的代码崩溃了

我试图从 lapack lib 调用 F77_NAME(dgeqrf) 函数,但程序崩溃了。 你觉得有什么不对吗? #包括 #包括 #包括 我试图从 lapack lib 调用 F77_NAME(dgeqrf) 函数,但程序崩溃了。 你觉得哪里不对? #include <R.h> #include <R_ext/BLAS.h> #include <R_ext/Lapack.h> #include <stdio.h> double* getcopy(int* nrow, int* ncol, double* a) { double* copy = malloc(*nrow * *ncol * sizeof(double)); int i; for (i = 0; i < *nrow * *ncol; i++) { copy[i] = a[i]; } return copy; } int main() { int m = 3, n = 3; double a[] = {12, -51, 4, 6, 167, -68, -4, 24, -41}; double* acopy = getcopy(&m, &n, a); double tau[3]; int info = 0; int i; int one = 1; double* work = malloc(max(&one, &m) * sizeof(double)); int lwork = 0; int lda = max(&m, &n); F77_NAME(dgeqrf)(&m, &n, acopy, &m, tau, work, &lwork, &info); for (i = 0; i < m * n; i++) printf("%f", acopy[i]); return 0; } 来自 dgeqrf.f 源代码: * LWORK (input) INTEGER * The dimension of the array WORK. LWORK >= max(1,N). * For optimum performance LWORK >= N*NB, where NB is * the optimal blocksize. * * If LWORK = -1, then a workspace query is assumed; the routine * only calculates the optimal size of the WORK array, returns * this value as the first entry of the WORK array, and no error * message related to LWORK is issued by XERBLA. 因此,如果您想要最佳运行,您首先需要调用该函数并将“lwork”设置为-1,以获得“work”的最佳大小。

回答 1 投票 0

使用 CMake 检测 BLAS/LAPACK 供应商

所以我的代码希望在不同的 BLAS/LAPACK 供应商发生时包含不同的头文件。是否有任何预定义的宏或类似的东西让我检查它?

回答 2 投票 0

Quantum Espresso v6.8 无法在 MacBook (Monterey 12.0.1 (21A559)) 上编译

按照此处发布的指南,我尝试在我的MacBook(macOS:Monterey 12.0.1(21A559))上安装最新版本的量子浓缩咖啡(v6.8)。因此,我安装了所有必备软件包...

回答 2 投票 0

如何从C调用CLAPACK?

我正在尝试从 C 调用 CLAPACK 函数。我下载了 CLAPACK-3.2.1(来自 http://www.netlib.org/clapack/clapack.tgz)并按照此处的说明进行操作(http://people.tgz)。 sc.fsu.edu/~%20jburkardt/c_src/

回答 2 投票 0

如何在C中使用lapack构建特征向量

我正在用C代码中的lapacke研究特征向量,我想知道当特征值有虚部时它们是如何构建的?我的意思是我必须在相同的维度上进行 malloc,就好像它们是...

回答 1 投票 0

Ipopt 找不到 Lapack,需要安装它

我正在尝试在 Linux 计算机上安装 Ipopt,但是我没有任何 sudo 权限。我认为 Lapack 安装成功(有什么方法可以检查吗?)。由于没有 sudo 权限...

回答 1 投票 0

在没有 sudo 权限的情况下安装 Lapack

我想在 Linux 机器上安装 Lapack,但是我没有任何 sudo 权限。 我当前所做的是将 git 存储库克隆到我的 $HOME 文件夹中。在 $HOME/lapack/ 中运行“make”。然后cr...

回答 1 投票 0

LAPACK zgemm op(A) 尺寸

在此来自 netlib 的链接中,它将 M 指定为: 输入时,M 指定矩阵的行数 矩阵 C 的 op( A ) 和 M 必须至少为零。

回答 1 投票 0

在 EMR 上使用 OpenBLAS Spark

尝试在 EMR 实例上运行 Spark 2.1.0 中的 MlLib ALS 算法时,我不断收到臭名昭著的警告: 警告 BLAS:无法从以下位置加载实现:com.github.fommil.netlib.NativeSyst...

回答 2 投票 0

sklearn 在新的 conda 环境中找不到 lapack

问题: 为什么 sklearn 无法检测到其 liblapack 依赖项以及如何解决此问题? 背景: 我在以下 conda 环境(torch-nightly)中安装了 scikit-learn(请注意,liblapack

回答 1 投票 0

重复单精度复矩阵向量乘法(速度和精度提高)

我将一个长时间运行的函数简化为一系列“简单”的矩阵向量乘法。矩阵没有变化,但是向量却很多。我做了一个测试

回答 1 投票 0

ATLAS:计算矩阵逆时出错

我在使用 ATLAS 库计算矩阵的逆时遇到错误。 我得到错误输出: 在进入 DTRSM 参数编号 9 时具有非法值 中止(核心转储) 我有一个

回答 1 投票 0

为什么我会得到两个相同类型的主元值进行 LU 分解? LAPACK sgetrf

我正在使用 LAPACK,我得到了两种相同的主值。我正在使用 LAPACK 例程 sgetrf 来计算 LU 分解 A = L*U*P 下面的 C 代码给出了与我的 MA 相同的结果...

回答 1 投票 0

“未定义对‘dgesvx_’的引用”,尽管 LAPACK 已安装并定位

我正在尝试运行一个包含函数 dgesvx() 的 Fortran 脚本。我运行脚本为 gfortran -lblas -llapack -o cft numtype.o cfapprox.o cftest.o 输出是 gfortran -lblas -llapack ...

回答 1 投票 0

LAPACK 的奇异值分解:大矩阵的问题

我正在使用 LAPACK 的 C 接口来计算矩阵的奇异值分解 (SVD)。为此,我使用例程 dgesvd_。 我创建了一个简单的 C++ 脚本,它创建了一个...

回答 2 投票 0

macOS Monterey - 安装 Scipy 错误“未找到 lapack/blas 资源”

我已经使用brew install pyenv在Mac上安装了pyenv。我想运行一个需要安装 scipy==1.5.4 的 python 项目。我尝试使用命令 pip3 install scipy==1.5.4...

回答 2 投票 0

将自定义 Matrix 类传递给 Lapack 子例程

我想要 lapack 包来实现一些非常有用的功能,但我不喜欢/无法自己很好地实现这些功能。问题是我无法将自定义 Matrix 类传递给 lapack

回答 1 投票 0

为什么这个简单的 Fortran 矩阵求逆代码没有返回 LAPACK 的预期值?

我正在尝试使用 LAPACK 包来计算矩阵的逆。为此,我使用了例程 sgetrf 和 sgetri。但是,在测试代码时,它不会返回预期的 v...

回答 0 投票 0

Fortran 中 LAPACK 的错误结果

我通过粘贴此链接中的命令在 Msys2 中安装了 LAPACK:https://packages.msys2.org/package/mingw-w64-x86_64-lapack。然后我将库与 -llapack -lblas 链接起来。一切似乎都有效

回答 1 投票 0

最新问题
© www.soinside.com 2019 - 2024. All rights reserved.