按fortran顺序填充BLAS dger填充矩阵？
BLAS（Basic Linear Algebra Subprograms）是一组用于执行基本线性代数运算的子程序集合。其中，dger是BLAS中的一个子程序，用于计算矩阵的外积。
按fortran顺序填充BLAS dger填充矩阵意味着按照Fortran语言中的内存布局规则来填充矩阵。Fortran语言中，多维数组在内存中是以列优先（column-major）的方式存储的。
对于一个m×n的矩阵A，按fortran顺序填充的BLAS dger填充矩阵可以表示为一个(m+1)×(n+1)的矩阵B。其中，B的第一列和第一行分别用于存储向量x和向量y，而B的剩余部分则用于存储矩阵A。
填充矩阵B的规则如下：
- B的第一列（除第一个元素）存储向量x的元素，即B(i+1,1) = x(i)，其中i为0到m-1。
- B的第一行（除第一个元素）存储向量y的元素，即B(1,j+1) = y(j)，其中j为0到n-1。
- B的剩余部分存储矩阵A的元素，即B(i+1,j+1) = A(i,j)，其中i为0到m-1，j为0到n-1。
BLAS dger填充矩阵的优势在于利用了Fortran语言中的内存布局规则，可以提高矩阵运算的效率。它适用于需要频繁进行矩阵外积计算的场景，如线性回归、矩阵分解等。
腾讯云提供了适用于云计算的各种产品和服务，包括计算、存储、数据库、人工智能等。在BLAS dger填充矩阵的应用场景中，可以使用腾讯云的弹性计算服务（Elastic Compute Service，ECS）来部署和管理计算资源，使用腾讯云对象存储（Tencent Cloud Object Storage，COS）来存储矩阵数据，使用腾讯云人工智能平台（Tencent Cloud AI）来进行矩阵运算和分析。
以下是腾讯云相关产品和产品介绍链接地址：
- 弹性计算服务（ECS）：https://cloud.tencent.com/product/cvm
- 对象存储（COS）：https://cloud.tencent.com/product/cos
- 人工智能平台（AI）：https://cloud.tencent.com/product/ai
相关·内容
1回答
我正在寻找一个用fortran 90/95处理稀疏矩阵的库。我只需要非常基本的运算，比如矩阵向量乘法。你建议我用什么？
我已经搜索过了，一个分机号(？)BLAS称为“稀疏blas"，记录在blast技术论坛规范的第3章中：
，
看起来很理想。根据该文档，库应该有一个fortran 95接口。但是，我在任何我能够下载的地方都找不到这个规范的实际实现。我得到的印象是fortran 95稀疏blas规范并没有真正在任何地方实现？不管怎样，它不是开源的。
我还发现了一个名为sparsekit的库：
有没有人用过这两个库中的任何一个，或者fortran 90/95的其他稀疏矩阵库？我想要关于使
浏览 0提问于2012-09-08得票数 16
回答已采纳
就速度和易用性而言，在分布式矩阵上执行BLAS或稀疏BLAS操作时，您推荐使用哪些c++库(或Fortran中具有c++接口的库)？
浏览 0提问于2011-09-11得票数 0
我正在考虑用Intel进行一些计算，特别是Fortran中程序的矩阵向量稀疏BLAS函数。
我可以用矩阵来表示我的计算，这些矩阵恰好是稀疏的和斜对称的。
据我所见，稀疏BLAS对于一般矩阵和对称矩阵都有稀疏函数，所以我想知道是否有一种方法来处理稀疏的斜对称矩阵，因为我认为它会减少内存占用。
浏览 16提问于2022-07-01得票数 1
回答已采纳
2回答
我在找一个Fortran编译器。我听说过gfortran和Intel Fortran。我听说情报更快了。
必须支持Fortran 90
Fortran 95/2003/2008的更多特性支持得越好。
必须支持OpenMP
必须支持MPI
必须支持某种形式的BLAS
应该输出良好的错误信息。我刚从Fortran开始，所以我会犯很多错误。
重要的是，编译器应该产生高效的、高度优化的可执行程序集，特别是对于矩阵操作(因此，如果它正在对其最喜欢的BLAS库进行严格的优化，这将是很好的)。
编译器可以是免费的，也可以是商业的。
OS: Linux
浏览 0提问于2014-02-12得票数 11
回答已采纳
1回答
我正在寻找一个关于如何将BLAS或LAPACK函数合并到我通过Xcode开发的Objective C程序中的好指南。我在网上能找到的BLAS/LAPACK程序的唯一来源都是用Fortran写的--这很让人困惑。如何在Objective C中使用BLAS/LAPACK和静态代码？有没有在线的例子/指南？我希望能够在矩阵元素具有复数值的情况下执行矩阵运算。
干杯,
阴影
浏览 1提问于2011-10-08得票数 5
所以我试图帮助另一个So用户，在这个过程中，我不能创建一个Cython程序来在NumPy之外做一些简单的事情，这迫使我使用GIL。因此，这使得使用OpenMP (多核)变得不可能。然后我发现了一个有趣的帖子，而你可以从SciPy直接将Fortran库导入到Cython代码(BLAS，LAPACK)中，这些代码与NumPy一起安装，在我的例子中是Intel MKL的等效函数。我所要做的就是简单地将2个向量乘以1000x1维，再将另一个向量转置，得到一个1000x1000矩阵。但是我找不到相关的Fortran例程(相当于NumPy乘法)来完成这项工作。所有的例程似乎都是做矩阵乘法。所以现在SciP
浏览 1提问于2017-06-20得票数 2
1回答
我想使用Fortran95中的稀疏Blas来创建矩阵，并且我正在使用点条目构造。使用以下命令创建矩阵后
调用duscr_begin(n，n，a，istat)
这里a是n乘以n的矩阵的句柄。在向其中插入值之后，我如何使用它的句柄a来查看最终的矩阵？因为我想使用矩阵进行一些其他操作，所以我想看到三个向量(稀疏)形式的矩阵(row_index，Col_index，Value)。
有关此稀疏Blas的详细信息，请参阅第3章，可在此处查看
浏览 2提问于2012-05-30得票数 0
回答已采纳
1回答
我对BLAS很陌生，并试图编译我的Fortran程序，该程序使用BLAS进行矩阵矩阵乘法。我有两个矩阵A和B。
integer, parameter :: RK = selected_real_kind(10,2)
real(kind=RK),pointer :: A(:,:), B(:,:), C(:,:)
然后，我把一些随机的值放进去，然后把它乘以：
call dgemm('N', 'N', N, N, N, 1.0_RK,
A, N, B, N, 0.0_RK, C, N)
N是两个矩阵的大小。所有参数都正确吗？
我无法编译并获得以下
浏览 3提问于2014-10-05得票数 1
回答已采纳
好的，我想比较一下关于稠密矩形矩阵不对称的矩阵逆的结果。通常使用DGETRF和DGETRI Blas来求矩阵逆。假设2000×2000双精度矩阵A，我想求解求矩阵-逆.
然后我得到了SuperLU解算器。用稀疏求解器求逆矩阵的思想是求解A*X=I，其中I是恒等矩阵。如果有解，X就是逆矩阵。所以如果A2000x2000，Ainvers aka X2000x2000，identity I2000x2000
而在fortran中调用superlu的输入
call c_fortran_dgssv ( iopt, n, ncc, nrhs, acc, icc,ccc, b, ldb, factors,
浏览 7提问于2022-03-20得票数 -1
1回答
Numpy Toeplitz矩阵
、、、、
我尝试使用scipy.linalg.toeplitz(c, r=None)在Python语言中构建toeplitz矩阵。虽然我成功了，但我不能保持Fortran的顺序。我需要确保正在构造的toeplitz数组保持Fortran顺序，因为我调用的是BLAS函数。
有什么办法可以做到这一点吗？
浏览 5提问于2014-01-18得票数 2
是否有一种方法可以自定义CMake，以便make命令不会显示每个百分比的进度？例如，现在我的构建显示
[ 1%] Building Fortran object dep/CMakeFiles/depuqtk.dir/blas/cgerc.f.o
[ 2%] Building Fortran object dep/CMakeFiles/depuqtk.dir/blas/cgeru.f.o
[ 3%] Building Fortran object dep/CMakeFiles/depuqtk.dir/blas/cscal.f.o
[ 4%] Building Fortran object
浏览 3提问于2015-01-08得票数 1
1回答
我正在尝试使用CMake来构建一个依赖于blas的程序，我正在使用以下命令检测blas：
include (${CMAKE_ROOT}/Modules/FindBLAS.cmake)
问题是，FindBLAS需要一个fortran编译器，并抱怨
-- Looking for BLAS... - NOT found (Fortran not enabled)
由于我的机器上已经安装了blas (ATLAS Blas)，并且也安装了gfortran，我如何启用Fortran，或者是否有解决方法可以找到C语言的blas库？
浏览 2提问于2010-04-17得票数 4
回答已采纳
我有一个R包，它使用了很多Fortran子例程来进行递归线性代数计算的嵌套循环(很大程度上依赖于BLAS和LAPACK例程)。作为Fortran的接口，我使用了.Fortran函数。我刚刚读到关于在C/C++中编写的子例程中使用.Call而不是.C，这让我思考，当使用Fortran子例程时，通过用C编写一个简单的包装器，然后调用Fortran子例程，使用.Call接口是不是更好？
如上所述，我的Fortran代码在某种意义上非常简单，我只处理双精度或整数类型的多维数组。但我知道，我必须在R端写相当多的检查，以确保一切都不会因为我不小心忘记将某些矩阵的存储模式更改为整数或某些矩阵的维数更改等而崩
浏览 35提问于2013-03-02得票数 19
1回答
我有一个c++代码(用g++编译)，它创建并填充了一个由复杂的双值填充的blitz++矩阵(我称之为lseMatrix)。然后将矩阵与zgetrf fortran例程倒置为
zgetrf_(&n，&n，&((*lseMatrix)(0，0))，&n，&(iPiv(0)，&info)
我有我的旧机器(Ubuntu16.04下的单处理器)，它安装了g++和gfortran (版本5.4)和笔记本/blas(v.3.6.0-2ubuntu2)。代码在那里运行得很好，用zgetrf反演相对较大的矩阵需要大约10分钟。但是，当我使用我的新机器(在Ubu
浏览 2提问于2021-06-12得票数 0
1回答
分析器中的BLAS例程？
、、、
我正在尝试分析一些调用GEMM进行矩阵乘法的Fortran代码。
我曾经使用双精度，正如我所期望的那样看到了dgemm调用，但现在我已经切换到单精度，我在profiler中的任何地方都看不到sgemm。
老实说，我甚至不明白dgemm是怎么出现的。不需要使用-pg编译BLAS库来创建必要的符号吗？
这是为BLAS库使用intel MKL，为分析器使用gprof。
浏览 4提问于2018-02-28得票数 0
我正在尝试在Fortran中计算类似于加权矩阵内积的东西。我用来计算内积的当前脚本如下
! --> In
real(kind=8), intent(in), dimension(ni, nj, nk, nVar) :: U1, U2
real(kind=8), intent(in), dimension(ni, nj, nk) :: intW
! --> Out
real(kind=8), intent(out) :: innerProd
! --> Local
integer :: ni, nj, nk, nVar, iVar
! --> Computing i
浏览 1提问于2018-04-08得票数 2
2回答
需要稀疏矩阵库
、、
是否有复制密集BLAS功能的稀疏矩阵库？
我至少想要：
有效的SYR和SYRK (秩-k更新)，具有稀疏的输入(和可能的密集输出)，
稀疏对称存储的选项，
点，NRM2，转置和其他标准的好东西。
但是，命名和API不需要以BLAS为中心，只要它完成了工作。C或Fortran优先，但C++也可以。
我查看了、、和，但是他们要么不支持我所需要的，要么我就不知道他们是否支持。
浏览 0提问于2011-08-04得票数 2
18回答
我们正在计算其运行时受矩阵操作约束的东西。(如果有兴趣，请在下面提供一些详细信息。)这一经验引发了以下问题：
人们是否对矩阵数学Java库的性能有经验(例如，乘法、逆等)？例如：
我搜索了什么也没找到。
我们的速度比较的细节：
我们正在使用英特尔FORTRAN (IFORT) 10.1 20070913)。我们已经在Java (1.6)中使用Apache 1.2矩阵操作重新实现了它，并且它同意它的所有数字的准确性。(我们有理由在Java中使用它。)(Java双倍，Fortran real*8)Fortran: 6分钟，Java 33分钟，同一台机器。jvisualm
浏览 36提问于2009-02-09得票数 155
回答已采纳
我试图在Anaconda的Numba包中使用cuBLAS函数，并且有问题。我需要输入矩阵是C级的。输出可以按Fortran顺序进行。
我可以运行包提供的示例脚本。脚本有两个函数，gemm_v1和gemm_v2。在gemm_v1中，用户必须按照Fortran顺序创建输入矩阵。在gemm_v2中，它们可以传递给GEMM的cuda实现，并转移到设备上。我可以让这些例子与方阵一起工作。然而，我不知道如何使gemm_v2与非平方输入矩阵一起工作。有没有一种方法可以处理非正方形的C阶输入矩阵？
注：
理想情况下，在调用GEMM用于其他计算之后，输入和输出矩阵都将留在设备上(这是迭代方法的一部分)。
浏览 1提问于2017-07-25得票数 17
1回答
我想用一个下三角矩阵来计算A*x，x是向量。例如： 1 0 0
A = 2 4 0
3 5 6 使用压缩存储 A = (/ 1, 2, 3, 4, 5, 6/) 和 X = (/1, 1, 1/) 现在我想用BLAS函数做A*x，我应该把A转换回3x3矩阵吗？如果没有，你能给我一些提示吗？(我知道在fortran数组的内存中，A是连续存储的)
浏览 14提问于2020-06-18得票数 0
回答已采纳
2回答
我想问你一个问题，我经常在管理CSR/CSC*矩阵使用Fortran。假设我们有一个有N个实值的向量V。向量以前是按一定的大小分配的。现在，我们必须在索引P中添加一个值。一个蛮力代码是：
allocate(tempV(N))
tempV=V
deallocate(V)
allocate(V(N+1))
V=(/tempV(1:P-1), newValue, tempV(P:N)/)
deallocate(tempV)
显然，如果做到一次，这不是一个问题，但重复数千次将不那么有效。内存将填充和空4倍，我想要插入的每个值。
我想知道哪一个程序更适合处理这个问题。您可以提出普通的Fortran (首
浏览 2提问于2013-01-30得票数 1
回答已采纳
1回答
稀疏Blas示例
、、、
我需要使用稀疏算法来存储和乘以Fortran 90中的许多矩阵。我已经看到可以使用稀疏的Blas库来完成它。我想知道在哪里可以找到一些编译和使用这个库的例子(我对Fortran非常陌生)。我试着自己实现，我甚至不可能编写第一行。这是我的节目：
program main
use sparse_blas
IMPLICIT NONE
integer::istat
integer::i,j,NonZA,NonZB,nonZmul
double precision::A(4,4),B(4,4),mul(4,4)
!initialization of matrices
call DUSCR_BEGIN(
浏览 0提问于2014-03-20得票数 0
回答已采纳
我在FORTRAN例程之前已经开始使用LAPACK的LU分解例程来处理这个问题，我已经看到它在其他相关应用程序中被非常有效地使用。然而，我读到过，NumPy和SciPy线性系统求解器大多基于对C中相同的LAPACK/BLAS函数的直接调用，并且想知道切换到FORTRAN语言是否真的可以将计算时间减少到值得放弃更容易的、更高级的语言的水平。
如果有python求解器可以保证这种大小的矩阵(1,000到10000，平方)具有类似的性能，那么它们是什么呢？
我真的需要矩阵求逆，所以切换到迭代Ax=b解决方案不是一个选择。
浏览 1提问于2019-12-24得票数 4
1回答
BLAS定义了GEMV (矩阵向量乘法)2级运算。如何使用BLAS库执行向量矩阵乘法？
这可能是显而易见的，但我不知道如何使用BLAS运算来实现这个乘法。我本以为会有一个GEVM操作。
浏览 0提问于2015-04-06得票数 7
回答已采纳
用于BLAS ?gemm函数的各种?gemm参数可以在较大数组的切片上操作。例如，这个小型C程序对(200,200)矩阵的左上角和右上(100,100)个子矩阵进行矩阵乘法，并将结果存储在左下角(100,100)子矩阵中。
#include <stdlib.h>
#include <stdio.h>
#include <cblas.h>
int main()
{
double * a = malloc(200*200*sizeof(double));
for(int i = 0; i < 200*200; i++) a[i] = 1;
浏览 12提问于2015-06-12得票数 2
回答已采纳
我有一个Fortran子例程，它使用BLAS的子例程dgemm，dgemv和ddot，它们计算矩阵*矩阵，矩阵*向量和向量*向量。我有m*m个矩阵和m*1个向量。在某些情况下，m=1。似乎这些子例程在这些情况下不能很好地工作。它们不会给出误差，但结果似乎存在一些数值上的不稳定性。所以我必须写一些类似这样的东西：
if(m>1) then
vtuni(i,t) = yt(i,t) - ct(i,t) - ddot(m, zt(i,1:m,(t-1)*tvar(3)+1), 1, arec, 1)
else
vtuni(i,t) = yt(i,t) - ct(i,t) - zt
浏览 0提问于2009-06-17得票数 1
回答已采纳
我的问题是关于Java与编译代码的性能，例如C++/fortran/程序集在高性能的数值应用程序中。我知道这是一个有争议的话题，但我正在寻找具体的答案/例子。还有社区维基。我以前也问过类似的问题，但我想我说得很宽泛，没有得到我想要的答案。
双精度矩阵乘法(在blas库中通常称为dgemm )能够达到几乎100 %的峰值CPU性能(以浮点运算每秒计算)。
有几个因素使这一业绩得以实现：
缓存阻塞，以实现最大的内存局部性
循环展开，以最小化控制开销。
向量指令，如SSE
记忆预取
保证没有内存混叠
我已经看到了许多使用程序集、C++、Fortran、Atlas、供应商B
浏览 7提问于2010-02-27得票数 11
朱莉娅()主页上的基准测试结果显示，在"rand_mat_mul“基准测试中，Fortran比Julia/Numpy慢约4倍。
我不明白为什么fortran在从同一个fortran库(BLAS)??调用时速度更慢。
我还对矩阵乘法进化的fortran、julia和numpy进行了简单的测试，得到了类似的结果：
朱莉娅
n = 1000; A = rand(n,n); B = rand(n,n);
@time C = A*B;
运行时间: 0.069577896秒(分配了7MB)
IPython中的Numpy
from numpy import *
n = 1000;
浏览 1提问于2015-02-11得票数 2
1回答
计算v⊗v的特征值
、、、、
我需要计算对称矩阵的特征值。
M=v⊗v，
其中v是已知的一维数组(向量)，⊗表示外部乘积。
目前，我通过一个自定义函数显式计算desired，然后调用⊗例程dsyev获得所需的特征值。
是否有更有效的方法来计算这些特征值(可能直接从v?)在fortran 2003中使用lapack、blas还是自写算法？
浏览 8提问于2022-11-25得票数 1
1回答
我正在用Fortran 77编写一些科学代码，我正在讨论哪种方法更快。
基本上，我有一个MxN矩阵，我们称它为A。M大于N。在后面的代码中，我需要将转置(A)乘以一堆向量。
我的问题是，接受A，我自己转置它并存储它会更快，还是当我调用BLAS时，只给它转置标志？
谢谢！-Patrick
浏览 1提问于2010-09-14得票数 4
我希望执行矩阵指数，这显然是Expokit库是合适的。可悲的是，不像Lapack或OpenMP，这是不容易安装从Cygwin或明格Windows。因此，我已经从下载了这个库，但是一旦解压缩，它就完全由.f文件组成，而对如何使用它们几乎没有指导。我在网上发现的唯一网站是使用，因为它没有给出Expokit库如何链接的任何指示。
对于如何在Windows上安装Expokit，或者如果Windows不适合在Ubuntu上安装Expokit，我将非常感激。
进行ripero建议的更改，并在Ubuntu上运行“makesample_d”，我将得到如下日志。我假设这意味着示例已经成功编译，但我不知道这如何使
浏览 1提问于2018-07-03得票数 0
回答已采纳
我习惯于在cython中使用BLAS例程(来自scipy.linalg.cython_blas)，在那里输入经常被修改到位(例如，在例程中)。我也在尝试对scipy.linalg.blas.dger做同样的事情，但是尽管使用了overwrite_a=True，a还是没有受到影响
(dger执行一级更新，意思是a + np.outer(x, y))
In [29]: x = np.array([1, 0, 2])
In [30]: y = np.array([-1, 1, 0, 2])
In [31]: a = np.arange(12).reshape(3, 4)
In [32]:
浏览 3提问于2017-10-06得票数 2
回答已采纳
3回答
我想使用BLAS的SGEMV子程序在fortran中执行矩阵向量积。我有一个类似于以下代码的代码：
program test
integer, parameter :: DP = selected_real_kind(15)
real(kind=DP), dimension (3,3) :: A
real(kind=DP), dimension (3) :: XP,YP
call sgemv(A,XP,YP)
A是一个3x3矩阵，XP和YP是向量。在包含的模块中，可以看到以下代码：
PURE SUBROUTINE SGEMV_F95(A,X,Y,ALPHA,BETA,TRANS)
!
浏览 0提问于2019-01-14得票数 1
回答已采纳
我正在使用cython对一个矩形矩阵A进行秩1更新。我不能让dger按我想要的方式进行更新，所以我将它隔离在一个函数中：
from scipy.linalg.cython_blas cimport dger
cimport cython
def test_dger(double[:, :] A, double[:] x, double[:] y):
cdef int inc = 1
cdef double one = 1
cdef int n_= A.shape[0]
cdef int m = A.shape[1]
dger(&n, &
浏览 14提问于2017-02-24得票数 0
1回答
我有一个线性方程组，其中as =b，A是一个大矩阵，b也是已知的。矩阵A是用python建立的。现在我想把矩阵A反求得到x。
A和b通过共享对象传递给Fortran 90程序。我使用numpy.f2py编译了Fortran程序：
import numpy.f2py.f2py2e as f2py2e
import sys, os
sys.argv += "-lmkl_rt -c -m MKL_MODULE MKL_WRAPPER.f90".split()
f2py2e.main()
最后，我调用f90子例程：
MKL_MODULE.mkl_wrapper.call_dgelsd
浏览 3提问于2016-10-30得票数 0
回答已采纳
3回答
连接LAPACK/BLAS库
、、、、
背景：
我正在编写一个用C和Fortran 77混合编写的项目，现在需要将LAPACK/BLAS库链接到该项目(全部在Linux环境中)。所讨论的LAPACK是来自netlib.org的3.2.1版本(包括BLAS)。库是使用顶级Makefile编译的(make lapacklib和make blaslib)。
问题：
在链接期间，错误消息声称从LAPACK例程调用的某些(并非所有)BLAS例程是未定义的。这让我有些头疼，但当(在Makefile中)改变了要链接的库的出现顺序时，这个问题最终得到了解决。
代码：
在下面，(a)给出错误，而(b)没有。链接由(c)执行。
(a) LIBS =
浏览 5提问于2010-05-10得票数 1
1回答
在使用BLAS 矩阵乘法函数时，我注意到对于一个未初始化的结果矩阵C，当我像这样调用它时，我会在结果中得到NaNs：
DGEMM('N', 'N', M, N, K, 1.0, A, LDA, B, LDB, 0.0, C, LDC)
但是，如果我事先声明了ALPHA和BETA：
REAL*8 ALPHA, BETA
PARAMETER (ALPHA=1.0)
PARAMETER (BETA=0.0)
DGEMM('N', 'N', M, N, K, ALPHA, A, LDA, B, LDB, BETA, C, LDC)
那么
浏览 2提问于2012-10-14得票数 1
回答已采纳
我正在使用VisualStudio2019中的Armadillo (10.4.1)来完成一些矩阵操作。我使用了来自OpenBlas管理器的NuGet，但是一切都很慢。现在我想切换到最新版本的OpenBlas。我使用了最后一个(0.3.15)，并按照下面的教程: minGW：编译了它。
编译工作正常，但当我尝试矩阵乘法时，控制台中会抛出一个错误：在输入到DSPR2参数1时有一个非法值。
我确实定义了ARMA_USE_BLAS，ARMA_DONT_USE_WRAPPER。我试着和ARMA_BLAS_LONG，ARMA_BLAS_UNDERSCORE，ARMA_USE_FORTRAN_HIDDEN_
浏览 14提问于2021-05-18得票数 0
回答已采纳
我在看。
import numbapro.cudalib.cublas as cublas
blas = cublas.Blas()
n =100
A = np.random.random((n, n)).astype(np.float32)
B = np.random.random((n, n)).astype(np.float32)
C = np.zeros_like(A, order='F')
blas.gemm('T', 'T', n, n, n, 1.0, A, B, 1.0, C)
assert(np.allclose(np.d
浏览 0提问于2019-06-28得票数 0
我想确定任何BLAS Fortran库中整数函数IDAMAX的大小。它可以是4字节(i32lp64型号)或8字节(ilp64)。
知道这个大小可以确定预构建BLAS库中的整体整数声明( integer *4或integer*8)。
问题是下面程序的sizeof(*idamax_(&n，dx，&incx))总是返回4，尽管我希望MKL-integer*8blas返回8。
有什么意见吗？
#include <stdio.h>
void main() {
extern * idamax_(); // external Fortran BLAS function IDAMA
浏览 0提问于2016-06-11得票数 1
首先，我知道这个问题以前有过问题，我的情况不一样。
我将矩阵库用于Xcode中的一个研究项目。我对这个项目很感兴趣，到目前为止，我还没有遇到任何问题。然而，我最近做了一系列的更改，当我尝试重新编译该项目时，我得到了以下错误，尽管没有对Armadillo头文件做任何更改。
In file included from /Users/willandrew/Work/1-4th-Year/1Project/code/TPS-RPM(LAPACK)/TPS-RPM(LAPACK)/main.cpp:16:
In file included from /Users/willandrew/Work/1-4t
浏览 1提问于2014-03-13得票数 1
回答已采纳
2回答
我需要在我的c++代码中使用blas和blas，我想链接MKL、ACML或默认的linux系统中的lapack和blas。不幸的是，它们在c中有不同的约定。例如，MKL ( mkl_blas.h )中的zdotc是
zdotc(&result, &n, x, &incx, y, &incy);
并且调用zdotc的默认call和blas (fortran版本)是：
result = zdotc_(&n, x, &incx, y, &incy);
如果我希望我的代码使用MKL，ACML，或者默认的漏洞blas。我需要为这个写一个包装：
#i
浏览 6提问于2016-10-11得票数 2
回答已采纳
3回答
矩阵/张量三积？
、、、、
我正在研究的一个算法需要在几个地方计算一种矩阵三重乘积。
该运算采用三个维数相同的方阵，并产生一个3-指数张量.标记操作数A、B和C，结果的(i,j,k)-th元素是
X[i,j,k] = \sum_a A[i,a] B[a,j] C[k,a]
在numpy中，您可以使用einsum('ia,aj,ka->ijk', A, B, C)计算这个值。
问题：
这个操作有标准名称吗？
我可以用一个BLAS调用来计算这个值吗？
是否还有其他重优化的数值C/Fortran库可以计算这种类型的表达式？
浏览 2提问于2015-05-13得票数 10
回答已采纳
我正在尝试从转向fortran和c++。我主要处理的是巨大的图像，坐标方向是x:从左到右，y:自上而下，像素被存储在行中，这是图像格式中常见的。
Numpy说，它也按行存储矩阵。就好而言，它与fortran/c++保存图像是一样的。但是，请参阅下面的示例，灰色从黑色到白色依次排列。让我们创建第一行为123等的图像3x3。
mat = np.array([1,2,3,4,5,6,7,8,9], 'i1').reshape(3,3)
接下来，我使用下面的代码获取有关矩阵的信息。例程打印:在位置x=1, y=0的值，在内存中的位置，看数组是否被复制，大步和迭代器。
def pr(a
浏览 4提问于2022-01-17得票数 1
我目前正在翻译一些遗留的fortran代码，我很难理解代码中的某一行代码。编译器似乎也发现这一行奇怪，并抛出一个错误。据我所知，它试图通过1到9的递增顺序初始化数组，并以列的主要形式填充数组矩阵。
program arrayProg
integer :: matrix(3,3), i , j !two dimensional real array
matrix = reshape((/1:9:1/), (/3,3/))
end program arrayProg
这种语法在fortran中可以接受吗？(必须是因为它来自遗留代码)我是否误解了行的功能？
浏览 0提问于2017-10-11得票数 4
回答已采纳
2回答
无法从源头上构建八度
、、、、
当我尝试从源代码构建Octave (使用OpenBLAS优化的BLAS和LAPACK库)时，我会遇到以下错误：
checking for sgemm_ in /usr/lib/libopenblas_sandybridgep-r0.2.8.so... yes
checking whether LSAME is called correctly from Fortran... yes
checking whether ISAMAX is called correctly from Fortran... yes
checking whether SDOT is called correctly f
浏览 1提问于2013-11-06得票数 3
回答已采纳
我不理解Fortran中未格式化文件的格式。
例如：
open (3,file=filename,form="unformatted",access="sequential")
write(3) matrix(i,:)
将矩阵的列输出到文件中。我发现它在文件的两端都填充了4个字节，但是我真的不明白为什么，或者如何控制这种行为。有没有办法去掉填充物？
浏览 1提问于2012-01-06得票数 20
回答已采纳
我在用Armadillo。我有这些变量：
arma::Mat<float> m_matrix;
arma::cx_vec m_eigenvalues;
arma::cx_mat m_eigenvectors;
我想这么做：
void calculate_eigens ()
{
arma :: eig_gen (m_eigenvalues, m_eigenvectors, m_matrix);
}
但是eig_gen的函数原型不允许这些参数类型，它希望第三个参数是一个双精度矩阵。
gcc的错误如下：
错误:调用“eig_gen(arma：：cx_vec&am
浏览 3提问于2017-10-24得票数 2
回答已采纳
我刚开始使用数值线性代数，我刚刚开始使用LAPACK和BLAS。
是否有一个例程可以复制/转换打包和完整存储之间的对称矩阵？
我找到了dtrttp，我可以用它将双精度的全对称矩阵转换为填充存储。然而，这些例程是针对三角形矩阵的，因此相应的dtpttr只填充整个矩阵的一个三角形。另一半我怎么填呢？
浏览 0提问于2015-02-11得票数 4
回答已采纳
如果我有一个张量收缩A[a,b] * B[b,c,d] = C[a,c,d]，它具有B[b,c,d] = B[b,d,c]和C[a,c,d] = C[a,d,c]的性质，那么如何建立BLAS来利用这个对称性呢？
这里假设爱因斯坦求和表示法，即重复指数表示和。
sgemm 似乎是关于矩阵的对称性，而不是秩-3张量.
我可以尝试将张量B平面/重塑成一个更低维的数组，但是看起来平坦/重塑张量也需要时间，至少在Fortran中是这样。
浏览 5提问于2021-03-15得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
云直播活动推荐
腾讯云开发者
