Матрица Якоби, решение не продвигается дальше второй итерации

#include "stdafx.h"
#include <iostream>
#include "stdio.h"
#include <stdlib.h>
#include <fstream>
#include <iomanip>
#include <ctime>
#include <omp.h>
 
 
using namespace std;
 
void main()
{
    setlocale(LC_ALL, "Russian");
    omp_set_num_threads(100);
    int n, iteration, m = 0;
    fstream f1;
    double e = 0;
    bool p = 0;
    cout << "Размер пространства ";
    cin >> n;
    cout << "\nКол-во итераций ";
    cin >> iteration;
    f1.open("C:\\123\\matrix.txt", ios::in);
    /*f1.open("C:\\123\\matrix.txt", ios::out);
    for (int i = 0; i < 40000; i++)
    {
    m = rand() % 1000 + 1;
    f1 << m << " ";
    }*/
    double **array, **array2;
    array = new double *[n];
    array2 = new double *[n];
    for (int i = 0; i < n; i++)
    {
        array[i] = new double[n];
        array2[i] = new double[n];
    }
    for (int i = 0; i < n; i++) //Заполнение массива элементами по периметру
    {
        for (int j = 0; j < n; j++)
        {
            if ((i == 0) || (i == (n - 1)) || (j == 0) || (j == (n - 1)))
            {
                f1 >> m;
                array[i][j] = m;
            }
            else
                array[i][j] = 0;
            cout << array[i][j] << "\t";
        }
        cout << "\n";
    }
    cout << "\n";
 
    for (int k = 0; k < iteration; k++) // Основной цикл, все что происходит - происходит в нем
    {
        //#pragma omp parallel for 
        for (int i = 0; i < n; i++) //Циклы всех расчетов
        {
            for (int j = 0; j < n; j++)
            {
                if ((i == 0) || (i == (n - 1)) || (j == 0) || (j == (n - 1)))
                    array2[i][j] = array[i][j];
                else
                    array2[i][j] = (array[i - 1][j] + array[i + 1][j] + array[i][j - 1] + array[i][j + 1]) / 4;
                e = array2 - array;
                if ((e < 0.01) && (k != 0))
                    p = 1;
            }
        }
 
        for (int i = 0; i < n; i++) //Вывод в консоль
        {
            for (int j = 0; j < n; j++)
                cout << array2[i][j] << "\t";
            cout << "\n";
        }
        memcpy(array, array2, n*n);
        cout << "\n";
        if ((k != 0) && (p = 1))
            break;
    }
    system("pause");
}