我在模拟软件中工作,在阵列上完成的许多操作之一是按数字缩放矢量.
我有这样的代码：
//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(),但这个操作非常慢.
