программировал умножение матрицы на вектор, и заметил что умножение с использованием простых динамических массивов выполняется в два раза быстрее в среднем чем с использованием std::vector, с чем это связано? ведь проверки на выход за пределы массив нет, вроде бы и оверхеда нет никакого почти значительного по сравнению с массивами, почему все-таки столь различное время работы
#include <iostream>
#include <vector>
#include <ctime>
//#include <omp.h>
void fillVector(std::vector<int>& vector)
{
for (int i = 0;i < vector.size();++i) {
vector[i] = rand() % 10;
}
}
void fillMatrix(std::vector<std::vector<int>>& matrix)
{
for (int i = 0;i < matrix.size();++i) {
for (int j = 0;j < matrix[i].size();++j) {
matrix[i][j] = rand() % 10;
}
}
}
void matmul(const std::vector<std::vector<int>>& matrix,
const std::vector<int>& vector, std::vector<int>& result)
{
//#pragma omp parallel for
for (int i = 0;i < matrix.size();++i) {
result[i] = 0;
for (int j = 0; j < matrix[i].size(); ++j) {
result[i] += matrix[i][j] * vector[j];
}
}
}
void fillVector(int* vector, const int n)
{
for (int i = 0;i < n;++i) {
vector[i] = 1 + rand() % 10;
}
}
void fillMatrix(int** matrix, const int n, const int m)
{
for (int i = 0;i < n;++i) {
for (int j = 0;j < m;++j) {
matrix[i][j] = 1 + rand() % 10;
}
}
}
void matmul(int** matrix, int* vector, int* result, const int n, const int m)
{
//#pragma omp parallel for
for (int i = 0;i < n;++i) {
result[i] = 0;
for (int j = 0; j < m; ++j) {
result[i] += matrix[i][j] * vector[j];
}
}
}
int main()
{
const int n = 15000;
const int m = 15000;
std::vector<std::vector<int>> matrix(n, std::vector<int>(m));
std::vector<int> vector(m);
std::vector<int> result(n);
fillVector(vector);
fillMatrix(matrix);
clock_t startTime = clock();
matmul(matrix, vector, result);
clock_t endTime = clock();
std::cout << endTime - startTime << "msec\n";
int* vector1 = new int[m];
int** matrix1 = new int* [n];
for (int i = 0;i < n;++i) {
matrix1[i] = new int[m];
}
int* result1 = new int[n];
fillVector(vector1, m);
fillMatrix(matrix1, n, m);
startTime = clock();
matmul(matrix1, vector1, result1, n, m);
endTime = clock();
std::cout << endTime - startTime << "msec\n";
delete[] result1;
delete[] vector1;
for (int i = 0;i < n;++i) {
delete[] matrix1[i];
}
delete[] matrix1;
return 0;
}
matrix[i][j]
для обычных массивов это просто пара умножений и сложений. В случае вектора там немножко больше. Плюс студия любит вставлять внутрь даже[]
дополнительные проверки. я попробовал бы вынестиmatrix[i]
за пределы цикла также.int sum=0;..sum+= matrix[i][j] * vector[j];..result[i]=sum;
Так как чтение, запись в одно и тоже место это лишняя работа для доступа к памяти.