问题标签 [lapack]

我需要找到 $A \times B$ 的 SVD，其中 $A$ 是上三角矩阵，B 是下三角矩阵。在 LAPACK 中执行此操作的有效方法和例程是什么？

我正在尝试从 lapack lib 调用 F77_NAME(dgeqrf) 函数，但程序崩溃了。你觉得有什么问题吗？

我想在 ac 程序中使用 Lapack lib 中的 F77_NAME(dgeqrf) 函数计算 qr 分解。

对于矩阵 3x3 ：

我得到输出 3x3（R 矩阵和一些用于构造 Q 的向量的组合）（线性形式）：

然后我使用来自 Lapack 的 F77_NAME(dorgqr) 来提取 Q 矩阵，得到输出 3x3（线性形式）：

这是一个取自维基百科的例子，我的 Q 似乎与维基百科 Q 不同： http ://en.wikipedia.org/wiki/QR_decomposition#Example_2

fortran 和 c 数组表示之间的差异可能是原因吗？初始矩阵的转置会解决问题吗？

我想使用Blas和Lapack库来使用一些 rutines，但是我不知道如何在Visual C++ 2010.

在这种情况下如何使用它们？

我想在 Visual C/C++ 2010 中使用 LAPACK 和 BLAS（将 LAPACK 和 BLAS 库添加到 Visual C/C++ 2010），所以我在这里阅读了一个非常相似的问题（相同但使用 Visual C/C++ 2008），他声称他解决了这个问题：

我想出了我的问题的答案。问题是我知道我必须将包含 LAPACK .lib 文件的目录放在库目录列表中（工具 -> 选项 -> 项目和解决方案 -> VC++ 目录 -> 库文件）但忘记了我还必须将其添加到项目->（项目名称）属性->链接器->输入下的“附加依赖项”。

当我尝试遵循他的解决方案时，在得到Tools -> Options -> Projects and Solutions -> VC++ Directories Visual后，C/C++ 2010 说 VC++ 目录未使用...

那么现在如何添加 .libs 以便 Visual C/C++ 2010 识别它们？

假设我得到一个可以看作矩阵的结果表，就像在示例中一样（矩阵是 aStore Procedure或 aUDF应用于表中的结果SQL Server 2008）

所以假设我得到类似的东西：

我想操纵这种结果（从 Store 过程或 UDF 获得）outside the DBMS，我想以 2 种方式进行：

1.-Send the complete matrix作为某种二进制文件，我的意思是avoid passing through a text file

（然后我想使用or MATLAB）。LAPACKread this binary file and manipulate it

2.- 将这个完整的矩阵作为一个direct memory pointer或可能1 MB具有矩阵信息的大小的内存块发送，然后MATLAB或者LAPACK would access这个指针或那些内存块并能够操作矩阵。

• 这些程序是否可行？如果可行
• 步骤是什么？
• 所有这些都可以在一个中完成吗Store procedure，是否有必要使用 Visual C/C++ 或类似的想法？

它是什么LAPACK bindings我如何使用它们来读取那块内存？

• 如何制作包含此矩阵“表”的内存块SQL Server 2008？
• 如果那不可能，包含矩阵的内存指针是否可行？

我正在编写 c 代码以使用 R 中的 lapack lib 获取矩阵 A 的 QR 分解。我的结果与我使用 R 语言命令得到的结果不同。

这是lapack的事情还是我的代码中的错误？

对于矩阵（行专业）：

结果是：

我的结果：

在过去的一周中，我一直在构建一个庞大的静态库（>1000 个 C 文件），其中包含由独立第三方生成的文件。

目前，我正在编写一个将进行一些密集计算的应用程序，看来我需要添加 CLAPACK 库。我相信加速框架是要走的路，但我似乎无法让它表现出来：

如果没有框架，我有一堆链接器错误告诉我我错过了，其中包括 cblas_zgemm 和 dgetrf（总共有十几个分布在大约 30 个文件中）。

在我将框架添加到项目中并且没有更改我的代码一 iota 之后，cblas_zgemm 链接器错误消失了，但 dgetrf 仍然存在（尽管属于框架）。

如果我添加#include (or #import) <Accelerate/Accelerate.h>，我会收到超过 1000 个链接器错误，告诉我我正在重新声明一堆枚举器（例如 CblasTrans）并且 xyz 的类型存在冲突......错误消息似乎在重复，同时 dgetrf 仍未定义.

我做错了什么/我错过了什么？

谢谢

编辑：我添加后的完整错误消息是：

编辑2：添加框架后的原始链接器错误：

我正在使用 Octave 和 R 使用简单的矩阵计算 SVD 并得到两个不同的答案！代码如下：

R

八度

我是 Octave 和 R 的新手，所以我的第一个问题是我这样做对吗？如果是这样，哪一个是“正确的”？他们俩都对吗？我也在 numpy 中尝试过，并直接调用 LAPACK 函数 dgesdd 和 dgesvd。Numpy 给我一个类似于 Octave 的答案，调用 LAPACK 函数给我一个类似于 R 的答案。

谢谢！