我想知道在c/c++中迭代矩阵的最快方法是什么。
到目前为止,我想到的最好方法是将矩阵映射到一维数组。
然后使用指针算术运算,是否有其他更快的方法?
维度在运行时已知但不在编译时已知,矩阵已完全填充。
到目前为止,我想到的最好方法是将矩阵映射到一维数组。
然后使用指针算术运算,是否有其他更快的方法?
维度在运行时已知但不在编译时已知,矩阵已完全填充。
#include <iostream>
#include <time.h>
#define XMAX 500
#define YMAX 400
#define ZMAX 300
int main()
{
srand(0);
register double sum = 0;
register int i;
register int j;
register int k;
double *arr_ptr;
arr_ptr = new double[XMAX*YMAX*ZMAX];
for (i=0; i<XMAX*YMAX*ZMAX; ++i)
{
*(arr_ptr+i) = rand()/double(RAND_MAX);
}
clock_t start, finish;
start = clock();
for (i=0; i<XMAX; ++i)
{
for (j=0; j<YMAX; ++j)
{
for (k=0; k<ZMAX; ++k)
{
sum += *(arr_ptr+i*YMAX*ZMAX+j*ZMAX+k);
}
}
}
finish = clock();
std::cout << "sum: " << sum << "\telapsed: " << finish - start << std::endl;
std::cin.get();
delete[] arr_ptr;
}
*(array+i)
比array[i]
快不了多少... 编译器会以相同的方式翻译它。 - Simonarray[i]
被定义为与*(array + i)
完全相同,因此这两种形式之间没有区别。 - caf