Печать каждой итерации возведения в степень

#include "stdafx.h"
#include <iostream>
#include <ctime>
#include "windows.h"
#include "process.h"
 
using namespace std;
 
const int NUMBER_OF_THREADS = 1;
 
const int N = 3;
const int POW = 2;
 
int **matr_1;
int **matr_3;
 
int **CreateArray(int L){
 
    int i, j;
    int **arr = new int*[L];
    for (i = 0; i<L; i++)
        arr[i] = new int[L];
 
    for (i = 0; i<L; i++)
        for (j = 0; j<L; j++)
            arr[i][j] = 0;
 
    return arr;
}
 
void DeleteArray(int **arr, int L){
 
    for (int i = 0; i<L; i++)
        delete[] arr[i];
    delete[] arr;
}
 
void PrintArray(int **arr, int L)
{
    for (int i = 0; i < L; i++){
        for (int j = 0; j < L; j++)
            cout << arr[i][j] << " ";
        cout << endl;
    }
}
 
void Pow(int **matr_3, int **matr_1, int L, int pow, int id){
 
    int i, j, k, p = 1;
 
    int start = N * id / NUMBER_OF_THREADS;
    int end = N  * (id + 1) / NUMBER_OF_THREADS;
 
    /*cout << "id=" << id;
    cout << "start=" << start;
    cout << "end=" << end;
    cout << endl;
*/
    int **matr_2 = CreateArray(N);
 
    for (i = start; i < end; i++)
        for (j = 0; j < L; j++)
            matr_3[i][j] = matr_1[i][j];
 
    while (++p <= pow){
        for (i = start; i < end; i++)
            for (j = 0; j < L; j++)
            for (k = 0; k < L; k++)
                matr_2[i][j] += matr_3[i][k] * matr_1[k][j];
 
            for (i = start; i < end; i++)
            for (j = 0; j < L; j++){
                matr_3[i][j] = matr_2[i][j];
                matr_2[i][j] = 0;
            }
            //PrintArray(matr_3, N);
            //cout<<endl;
    }
    DeleteArray(matr_2, L);
}
 
unsigned __stdcall threadMult(void *param){
 
    int id = (int)param;
 
    Pow(matr_3, matr_1, N, POW, id);
 
    return 0;
}
 
int main()
{
    matr_1 = new int *[N];
    matr_3 = new int *[N];
 
    for (int i = 0; i<N; i++)
        matr_1[i] = new int[N];
 
    for (int i = 0; i<N; i++)
        matr_3[i] = new int[N];
 
    srand(time(0));
 
    for (int i = 0; i<N; i++)
        for (int j = 0; j<N; j++)
            matr_1[i][j] = rand() % 3;
    
    PrintArray(matr_1, N);
    cout << endl;
 
    clock_t start, end;
    start = clock();
 
    HANDLE hThreads[NUMBER_OF_THREADS];
    unsigned id;
    for (int i = 0; i < NUMBER_OF_THREADS; i++)
        hThreads[i] = (HANDLE)_beginthreadex(NULL, 0, threadMult, (void *)i, 0, &id);
    WaitForMultipleObjects(NUMBER_OF_THREADS, hThreads, TRUE, INFINITE);
    for (int i = 0; i < NUMBER_OF_THREADS; i++)
        CloseHandle(hThreads[i]);
 
    PrintArray(matr_3, N);
 
    end = clock() - start;
    cout << endl << "runtime: " << end / (double)CLOCKS_PER_SEC << endl;
    cout << endl;
 
    DeleteArray(matr_1, N);
    DeleteArray(matr_3, N);
 
    //system("pause");
 
    return 0;
}