Исключение переполнения стека, возникающее при вызове метода

Я работаю над проблемой в C ++, но получаю исключение переполнения стека, и я не могу сказать, почему. Основной метод вызывает problem28 (), но первая строка должна вывести «check» в мой вывод, чего не происходит. Если я определяю размер сетки как 501 или меньше, он работает нормально, но что-то еще, и это вызывает исключение переполнения стека.

Любая помощь будет оценена.

#define right 0
#define down 1
#define left 2
#define up 3
#define gridsize 1001

int* next(int row, int col, int dir) {
int* newPos = new int[2];
newPos[0] = row;
newPos[1] = col;
switch(dir) {
case right:
newPos[1] += 1;
break;
case down:
newPos[0] += 1;
break;
case left:
newPos[1] -= 1;
break;
case up:
newPos[0] -= 1;
break;
}
return newPos;
}

int problem28() {
cout << "check" << endl;
int grid[gridsize][gridsize];
for (int i = 0; i < gridsize; i++)
for (int j = 0; j < gridsize; j++)
grid[i][j] = 0;
int* pos = new int[2];
pos[0] = pos[1] = gridsize / 2;
int dir = right;for (int i = 1; i <= 1001; i++) {
grid[pos[0]][pos[1]] = i;
pos = next(pos[0], pos[1], dir);
int* npos;

npos = next(pos[0], pos[1], (dir + 1) % 4);
if (grid[npos[0]][npos[1]] == 0)
dir = (dir + 1) % 4;
}
cout << "generated grid" << endl;

int total = 0;
for (int i = 0; i < gridsize; i++) {
total += grid[i][i];
total += grid[i][gridsize - i - 1];
}
total -= grid[gridsize / 2][gridsize / 2];

return 0;
}

int main() {
problem28();

system("pause");
return EXIT_SUCCESS;
}