c – MKL或BLAS例程将矢量乘以标量不合适的位置

我在模拟软件中工作,在阵列上完成的许多操作之一是按数字缩放矢量.

我有这样的代码

//Just some initialization code,don't bother about this part
int n = 10000;
std::vector<double> input(n,42.0);
std::vector<double> output(input.size());

double alpha = 69.0;

//the actual calculation:
for (size_t i = 0; i < n; ++i) {
    output[i] = input[i] * alpha;
}

我有MKL库,所以如果我的计算是“就地”完成的,可以写下面的内容

cblas_dscal(n,alpha,&input[0],1);

但是,这会改变输入变量,这不是我想要的.

我尝试使用mkl_domatcopy(),但这个操作非常慢.

解决方法

我想出的解决方案是调用cblas_dcopy()然后调用cblas_dscal().

它不是所有世界中最好的,但它仍然比原始循环更快.

相关文章

本程序的编译和运行环境如下(如果有运行方面的问题欢迎在评...
水了一学期的院选修,万万没想到期末考试还有比较硬核的编程...
补充一下,先前文章末尾给出的下载链接的完整代码含有部分C&...
思路如标题所说采用模N取余法,难点是这个除法过程如何实现。...
本篇博客有更新!!!更新后效果图如下: 文章末尾的完整代码...
刚开始学习模块化程序设计时,估计大家都被形参和实参搞迷糊...