MKL矩形矩阵就地转置:不使用多核？

文章/答案/技术大牛

发布

2回答

c++、c、intel、blas、intel-mkl

我想要一个非常大的矩阵的地方内存转置。我正在使用mkl_simatcopy。但我观察到一些性能问题，而转置就地。我目前使用的是具有72个物理核心和RedHat操作系统的英特尔至强处理器E7-8867 v4 @2.40 and。我的观察是，当我执行转置操作时，只使用了一个核心，并且它没有使用所有的核心。我已经尝试了所有的环境变量，如MK_NUM_THREADS、MKL_DYNAMIC="FALSE"

浏览 70提问于2019-05-13得票数 0

1回答

利用MKL库确定特征值的性能

c++、matrix、eigenvalue、intel-mkl、eigenvector

我的目标是找到mxn大小的输入矩阵的特征值和向量。由于它是一个矩形矩阵，所以通过做一个转置，然后用源矩阵进行矩阵相乘，把它转换成平方对称矩阵。在此之后，我使用Intel MKL库来确定平方对称矩阵的特征值和向量。不幸的是，对于1000000 x 100000的平方对称矩阵大小，所需时间为35~40秒，这不适合我的应用。

浏览 4提问于2016-11-29得票数 0

1回答

在Cython中使用Fortran NumPy操作和NOGIL，什么是相当于NumPy NumPy的Fortran库？

python、scipy、fortran、cython、blas

所以我试图帮助另一个So用户，在这个过程中，我不能创建一个Cython程序来在NumPy之外做一些简单的事情，这迫使我使用GIL。因此，这使得使用OpenMP (多核)变得不可能。然后我发现了一个有趣的帖子，而你可以从SciPy直接将Fortran库导入到Cython代码(BLAS，LAPACK)中，这些代码与NumPy一起安装，在我的例子中是Intel MKL的等效函数。我所要做的就是简单地将2个向量乘以1000x1维，再将另一个向量转置，得到一个1000x1000矩阵。但是我找不到相

浏览 1提问于2017-06-20得票数 2

3回答

R的REvolution

r、mathematical-packages

它似乎是一个用于高性能矩阵计算的库集合。显然，它似乎真的很有效。例如，在使用REvolution进行矩阵转置时： user system elapsed

浏览 2提问于2009-11-17得票数 5

回答已采纳

3回答

高效的矩阵转置算法

arrays、algorithm、matrix、transpose

我需要转置一个由char数组表示的方阵。有没有一种方法可以以低于o(n^2)的复杂度执行它？不管怎样，最高效的缓存方式是什么呢？

浏览 0提问于2015-01-19得票数 3

1回答

多核硬件上的numpy

python、performance、numpy、multicore、linear-algebra

关于让numpy使用多核(在英特尔硬件上)来处理内部和外部向量积、向量矩阵乘法等方面的最新技术是什么？如果有必要，我很乐意重新构建numpy，但在这一点上，我正在寻找在不更改代码的情况下加快速度的方法。作为参考，我的show_config()如下所示，我从未观察到numpy使用多个内核： libraries = ['lapack', 'ptf77blas"\\"3.9.16\\""

浏览 1提问于2011-05-13得票数 17

回答已采纳

1回答

变换在路径和矩形中的工作方式不同

c++、gdi+

我正在使用GDI+和C++在自定义控件上绘制一个矩形和一条路径。这是我的代码。myRectF);gdiGraphics.FillPath(new SolidBrush(Color.Green), myPath);在路径中绘制的矩形是它应该在的位置。

浏览 0提问于2012-06-04得票数 0

回答已采纳

5回答

Matlab中的快速矩阵乘法

performance、matlab、matrix、matrix-multiplication

我需要在Matlab中进行一个非常大的矩阵/向量乘法："A“是一个655360乘以5的实数矩阵，它不一定是稀疏的；"B”是一个655360乘以1的实数向量。我的问题是如何有效地计算: B'*A。通过一点搜索，我找到了James Tursa的一个有趣的Matlab工具箱，我希望它能提高上面的矩阵乘法性能。经过几次尝试，我只能获得比Matlab本机矩阵乘法更小的增益。

浏览 4提问于2011-10-04得票数 3

1回答

QGraphicsItem -项转换失败

c++、qt、transformation、qgraphicsview、qgraphicsitem

我有一个在QGraphicsItem中使用的矩形。项目被添加到场景中，并且它有一个预定义的转换，这样它就可以正确地映射矩形。然而，对于另一个转换矩阵，我遇到了这个有趣的结果：QTransform(type=TxProject, 11=-4.51027我认为这一定是框架中的一个bug，因为它可以很好地转换一个点，但是它不能转换项目(我还有很多椭圆，rects不能被转换)。但是当我试图繁殖时，我做不到；它只发生在这个矩阵

浏览 0提问于2014-12-06得票数 0

回答已采纳

3回答

关于使用并行化解决问题的一般问题

c、openmp、mpi

让我们假设我们的任务是用MPI和/或OpenMP实现一些矩阵算法。在某些情况下，如OpenMP或MPI中的错误共享，通信依赖于矩阵维度(列在进程之间循环分布)，这会导致一些问题。例如，通过转置矩阵来解决这种情况，这是否是一种好的和常见的尝试，因为这将减少必要的通信，甚至避免虚假共享问题？在此之后，您将撤消转置。当然，假设这会带来更好的速度。

浏览 0提问于2011-03-27得票数 2

5回答

在C++中求解正规方程组

c++、boost、linear-regression、eigen、lapack++

我想解线性方程组：A是一个n x m矩阵(不是平方)，b和x都是n x 1向量。当A和b已知时，n为50-100阶，m约为2(换句话说，A可以是最大的100x2)。我知道x：$x = \inv(A^T A) A^T b$的解决方案我几乎没有找到解决这个问题的方法: uBLAS (Boost)、Lapack、Eigen等等，但我不知道使用这些包的x的计算时间有多快

浏览 17提问于2011-01-03得票数 7

回答已采纳

3回答

在python中存储大型对称稀疏矩阵的最有效方法

python、matrix、numpy、sparse-matrix

这将需要存储、操作、调整大小，并且(在某些情况下)可能对非常大的稀疏矩阵进行对角化。我希望能够有由0和几个(比方说<5) 1组成的行，并一次添加几个行(按照所使用的cpus数量的顺序)。

浏览 3提问于2011-08-29得票数 5

回答已采纳

3回答

就地重新排序简单2D矩阵

c++、algorithm

我有一个简单的2D (行，列)矩阵，我目前根据下面的算法重新排序，使用另一个数组作为最终容器来交换项目。{ } 注意: coef是一个通过下标进行双重访问的指针，_matrixCoef是一个矩阵助手类，它使用一个由运算符(row，col)进行双重访问的向量。编辑2：下面是上面通过操作符重载(row，col)访问的矩阵类： struct Matrix //&#x

浏览 40提问于2021-11-21得票数 1

回答已采纳

8回答

什么是用于矩阵运算的好的C++库

c++、optimization、matrix

我需要对矩阵做乘法运算。我正在寻找一个可以快速完成这项工作的库。我使用的是Visual C++ 2008编译器，并且我有一个核心的i7 860，所以如果这个库针对我的配置进行了优化，那么它就是完美的。

浏览 1提问于2010-12-30得票数 13

1回答

将数组(k，)或(k，n)乘以一维数组(k，)

python、arrays、numpy、matrix

我有时使用一维数组：或2D阵列(用scipy.io.wavfile读取的单或立体声信号)： A = np.array([[1, 2], [3, 4注意:当使用scipy.io.wavfile读取单声和立体声时，这个问题自然会出现。

浏览 1提问于2018-09-04得票数 3

回答已采纳

1回答

gsl、pygsl、lmsder、内存泄漏、输出不一致

python-2.7、memory-leaks、io、garbage-collection、pygsl

要使用pygsl进行参数拟合，我需要定义函数及其雅可比矩阵。然后，multifit_nlin的fitter lmsder在拟合过程中需要时调用这两个函数。当我调用雅可比矩阵时，它会产生一个数字矩阵。我可以将这个矩阵输出到屏幕上，我看到这个数字是正确的。接下来，我定义了一个lmsder类，并使用lmsder.set命令对其进行初始化。我使用lmsder.getJ()命令将雅可比矩阵输出到屏幕上，我看到的数字与之前相同。当然，这不是我想要对我的代码做的事情

浏览 1提问于2013-08-02得票数 1

1回答

阴影映射中的偏差矩阵混淆

opengl、shadow、shadow-mapping

我对阴影贴图中的偏差矩阵感到困惑。根据这个问题：，偏差矩阵用于缩小并转换为0..1x和0..1y。所以我想，如果我们不使用偏差矩阵，纹理将只填充1/4的场景大小？这是真的吗？或者这里面有什么魔法？

浏览 0提问于2013-01-30得票数 3

回答已采纳

1回答

法线不是在对象空间中(从正常语义中读取)吗？

matrix、unity3d、shader

当将法线从对象空间转换到世界空间时，预定义的矩阵_Object2World就像魅力:(mul(_Object2World，法线))。旋转和缩放(一致)通过行和列中的颜色检查元素给出了期望的矩阵，但是为什么有些代码使用法线左乘逆矩阵_World2Object进行对象空间转换，是一个统一版本的东西(我使用4.5.2)，或者法线存在于对象空间以外的其他空间，以及为什么这两个矩阵不响应非均匀缩放。

浏览 8提问于2014-09-01得票数 3

回答已采纳

3回答

转置二维阵列

beginner、ruby、matrix

您对下面的代码有什么改进建议吗？ # number of rows n = arr[0].count arr.each_with_index do |row, index1| row.each_with_index do

浏览 0提问于2015-04-16得票数 5

回答已采纳

3回答

两个矩阵的行间相关系数

matlab、statistics、distance

我有下面的代码，假设X是10x7，Y是4x7，D是一个矩阵，存储每对向量之间的相关性，我将如何使用函数来简化它，因为它目前运行得相当慢。如果解决方案是使用xcorr2函数，有人能告诉我是怎么做的吗？

浏览 1提问于2014-01-09得票数 3

点击加载更多