问题标签 [lapacke]

根据官方用户指南，sgelsd 用于解决最小二乘问题

并允许矩阵 A 是矩形且秩亏的。并且根据sgelsd源码中的接口描述，b作为输入输出参数。当 sgelsd 完成时， b 存储解决方案。所以 b 占用 m*sizeof(float) 个字节。而解决方案 x 需要 n*sizeof(float) 个字节（假设 A 是 am*n 矩阵，b 是 am*1 向量）。

但是，当 n>m 时，b 的内存太小，无法存储解 x。如何处理这种情况？我没有从 sgelsd 源代码的注释中得到它。我可以只为 b 分配 n*sizeof(float) 字节并使用第一个 m*sizeof(float) 来存储 b 向量吗？

谢谢。

我试图发现一些库链接的含义，并在标题中找到了这个：

我不明白这些是做什么的，特别是在返回后缀的地方##_

谢谢

从 cython 调用 lapack 函数的最简单方法似乎是使用 lapacke。

我想在 libmkl_* 动态库中使用 Enthought Canopy 附带的 lapack 函数。

但是，我在 Canopy 的任何地方都找不到 lapacke.h。有什么方法可以在不单独安装 lapack 的情况下使用 lapacke？

我是 C++ 新手，我正在使用 Eclipse 编写脚本。我的操作系统是 Ubuntu。我需要为我的代码部分使用 LAPACKE 包。然而，我无法设法链接 Eclipse 和 LAPACKE。我正在尝试编译以下示例代码：

我无法编译代码，因为我的 Eclipse 抛出错误：“Udefined reference to LAPACKE_dgels”。我试图将 Eclipse 链接到 LAPACKE，为此我在 Eclipse 的“路径和符号”选项卡中添加了 LAPACKE 头文件的路径。任何人都可以帮助我解决这个问题需要做些什么吗？我应该错过一些东西......

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

我的文件my_example.c很简单：

当我编译它时，我得到了这个：

这是调用 clapack 函数的正确方法吗？数组传递是否正确？

（PS 我不想动态链接到现有的 clapack 安装）。

我想解决形式的线性系统A*X=B'，其中B'是 的转置B。A是一个方阵N-by-N，B 是N-by-M。在 lapack/lapacke 中，函数 LAPACKE_dgesv（参见此处的示例）用于求解形式为 的系统A*X=B，其中B将其视为多个右侧向量。是否可以通过重新排序其值然后求解来解决表单系统A*X=B'而不必创建Bas的副本？Z=B'A*X=Z

我正在尝试使用 lapacke 的 dgels 功能：当我将它与 malloc 功能一起使用时。它没有给出正确的值。当我使用 malloc 并创建矩阵时，谁能告诉我有什么错误？谢谢你

返回（信息）；}

请帮我！我花了几个小时在互联网上查找，但我还没有找到解决方案....

我正在尝试从C++函数中使用 call lapack 函数，但一开始我失败了。这是我的代码：

我知道“lapacke.h”是一个 C 头文件，所以我使用该extern "C"子句。但是，当我尝试编译这个微不足道的函数时，出现以下错误：

有谁知道是什么导致了这些错误？

非常感谢！

由于这是我第一次在这里发帖，我希望我不会犯太多错误......

我面临以下问题。当我运行这段代码时（很抱歉，我很抱歉）在 LAPACKE 调用期间，我在“解决”函数中遇到了分段错误。

问题是当我不尝试在主函数中打开文件“dataf”时，此代码运行时没有任何（可见）问题。

我想我搞砸了我的类型 Bob 和指针，但我想听听你的观点。

此致

蒂博

使用 Makefile：

编辑：我在原来的帖子中犯了一个错误。A[0] -> A[3] 的分配是双精度数，而不是指向双精度数的指针。

解决后更改 fclose 的位置似乎可以解决问题。但在我的真实代码中它不会改变任何东西......

EDIT2：问题似乎与 LAPACK 功能有关。我改变了这个功能（更适合我的问题）并且内存问题现在已经消失了。

这不是我原来问题的解决方案，但我以某种方式解决了它。感谢所有非常有用的回复和评论:)

EDIT3：我找到了问题的根源。它与LAPACK有关。在函数求解中，我为 ipiv 分配空间。但是 LAPACKE_dgttrs 将 init_A 中的函数 LAPACKE_dgttrf 给出的 ipiv 作为输入。

所以我在我的结构中为 ipiv 分配了一个指针，并在程序的其余部分中使用它。

我正在尝试使用矩形 A 矩阵解决 Xcode 中的 Ax = b，其中 A 和 b 的大小取决于我存储的不同数组的大小。我曾尝试使用 dgels，但由于某种原因，我收到此错误

这对应于工作双矩阵输入，我不知道为什么。

int numberOfPlays = (int)[playArray count]; int numberOfSignals = (int)[signalMappings 计数]；