Потоки Транспонирование матрицы

#include <iostream>
#include <thread>
template<typename T> bool matrix_alloc(T**& mat, int size);
template<typename T> void matrix_free(T** mat, int size);
void transponse(int** mat, int first, int last, int size);
 
int main(void){
    int  size = 26;
    int** mat;
    if(! matrix_alloc(mat, size))
        return 1;
 
    //заполняем матрицу
    for(int i = 0; i < size; ++i){
        for(int j = 0; j < size; ++j){
            mat[i][j] = 'A' + i;
            std::cout << char(mat[i][j]) << ' ';
        }
        std::cout << std::endl;
    }
    std::cout << std::endl;
 
    const int TN   = 3;
    const int step = size / TN + 1;
    std::thread ts[TN];
    for(int i = 0; i < TN; ++i){
        int p = i*step;
        ts[i] = std::thread(transponse, std::ref(mat), p, std::min(p + step, size), size);
    }
 
    for(int i = 0; i < TN; ++i)
        ts[i].join();
 
    //вывод транспонированной матрицы
    for(int i = 0; i < size; ++i){
        for(int j = 0; j < size; ++j)
            std::cout << char(mat[i][j]) << ' ';
        std::cout << std::endl;
    }
    matrix_free(mat, size);
    std::cin.get();
    return 0;
}
 
//транспонирование только квадратных матриц
void transponse(int** mat, int first, int last, int size){
    for(int i = first; i < last; ++i){
        for(int x = i; x < size; ++x)
            std::swap(mat[i][x], mat[x][i]);
    }   
}
 
template<typename T>
bool matrix_alloc(T**& mat, int size){
    if((mat = new (std::nothrow) T*[size]) == nullptr)
        return false;
 
    for(int i = 0; i < size; ++i){
        if((mat[i] = new (std::nothrow) T[size]) == nullptr){
            matrix_free(mat, i);
            return false;
        }
    }
    return true;
}
 
template<typename T>
void matrix_free(T** mat, int size){
    for(int i = 0; i < size; ++i)
        delete[] mat[i];
    delete[] mat;
}