Вернуть двумерный массив из функции и сделать проверку

#include "stdafx.h"
#include <stdlib.h>
#include <ctime>
#include <iostream>
#include <stdio.h> //Здесь находится объявление функции fopen и объявлен тип структуры FILE
#include <time.h>
#include <sstream>
using namespace std;
 
#define N 20
 
double** INV_matrix(int n, double **A, double **D)
{
    int i = 0; int r; double eps;
    r = n;
    eps = 0.00001;
 
    while (i < r)
    {//Выбор ведущего элемента A[v][i]
        int v = i;
        for (int j = i + 1; j < n; j++)
            if (abs(A[j][i])>abs(A[v][i])) v = j;
        if (abs(A[v][i]) < eps) r = i - 1;
        else
 
        {    //перестановка строк
            if (v != i)
            {
                for (int j = i; j<n; j++)
                {
                    double z = A[i][j]; A[i][j] = A[v][j]; A[v][j] = z;
                }
                for (int j = 0; j<n; j++)
                {
                    double z = D[i][j]; D[i][j] = D[v][j]; D[v][j] = z;
                }
            }
 
            //деление i-й строки на диагональный элемент
            double aii = A[i][i];
            for (int j = i; j<n; j++)
                A[i][j] /= aii;
            for (int j = 0; j<n; j++)
                D[i][j] /= aii;
 
            //вычитание уравнений
            for (int k = 0; k<n; k++)
                if (k != i)
                {
                    double c = A[k][i] / A[i][i];
                    for (int j = i; j<n; j++)
                        A[k][j] -= c*A[i][j];
                    for (int j = 0; j<n; j++)
                        D[k][j] -= c*D[i][j];
                }
            i++;
        }
    }
 
    
 
            
if (r==n) return D;
else return NULL;
 
}
 
int main()
{
    FILE *f1;
    f1 = fopen("test_INV.txt", "r");
 
    setlocale(LC_CTYPE, "rus");
 
 
    //Выделение памяти для двумерного массива (указатель на массив указателей)
    int n;
    fscanf(f1, "%d", &n);
    cout << "Размерность матрицы - " << n << endl;;
 
    //Выделение памяти для исходной матрицы
    double **A = new double *[n];
    for (int i = 0; i < n; i++)
        A[i] = new double[n];
 
    //Выделение памяти для обратной матрицы
    double **D = new double *[n];
    for (int i = 0; i < n; i++)
        D[i] = new double[n];
 
 
    cout << "Исходная матрица" << endl;
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < n; j++)
        {
            // j - количество столбцов; i - количество строк
            fscanf(f1, "%lf", &A[i][j]);
            printf("%2g ", A[i][j]); //Чтобы не выводить незначащие нули, можно использовать %g; при исп. %f - появляются незначащие нули в конце
                                     //cout << A[i][j] << " ";
        }
        cout << endl;
    }
    cout << endl;
 
    // Теперь в массив A[][] записана исходная матрица из файла
 
    double **Q = new double *[n]; //В массиве Q - исходная матрица для того, чтобы сделать проверку
    for (int i = 0; i<n; i++) 
        {
            Q[i] = new double [n];
            for (int j = 0; j<n; j++) Q[i][j] = A[i][j];
            
        }
        
 
    //Обратную матрицу заполняем как единичную (на гл. диагонали стоят 1)
 
    for (int i = 0; i<n; i++)
        for (int j = 0; j<n; j++)
            if (i==j) D[i][i] = 1;
            else D[i][j] = 0;
    
        
    FILE *f2;
    f2 = fopen("result_INV.txt", "w");
 
 
    //Здесь надо сделать проверку на то, что A*D = E
 
if (INV_matrix != NULL) /*ВОТ ЭТО МЕСТО МНЕ НЕПОНЯТНО*/
{
    double **D1 = INV_matrix (n, A, D); //Вызов функции; возвращаемое значение - обратная матрица D;
 
    double **C = new double *[n];
    for (int i = 0; i<n; i++)
        C[i] = new double [n];
 
    //Создание и заполнение единичной матрицы
    double **E = new double *[n];
    for (int i = 0; i<n; i++)
        E[i] = new double [n];
 
    for (int i = 0; i<n; i++)
        for (int j = 0; j<n; j++)
            if (i==j) E[i][i] = 1;
            else E[i][j] = 0;
 
 
    //Перемножение исходной и обратной матриц
 
    for (int i = 0; i<n; i++)
        for (int j = 0; j<n; j++)
        {
            double d = 0; 
            for (int q = 0; q<n; q++)
                d += Q[i][q]*D1[q][j];
            C[i][j] = d;
        }
 
 
        for (int i = 0; i<n; i++)
        for (int j = 0; j<n; j++)
            if (C[i][j] == E[i][j]) bool H = 1;
            else bool H = 0;
            
 
            if (bool H = 1)
            {
                fprintf(f2, "Обратная матрица \n");
                for (int i = 0; i<n; i++)
                { 
                    for (int j = 0; j<n; j++)
                    fprintf (f2, "%8.3f", D[i][j]);
                    fprintf(f2, "\n");
                }
            }
 
            else fprintf (f2, "Произошла ошибка, проверьте алгоритм!");
 
 
        fclose(f1); fclose(f2);
 
    delete[] *A; delete[] A; //Исходная матрица
    delete[] *D; delete[] D; //Обратная матрица
    delete[] *Q; delete[] Q; //Для дублирования исходной матрицы
    delete[] *C; delete[] C; //Для сравнения с единичной матрицей (при проверке)
    delete[] *E; delete[] E; //Единичная матрица
    delete[] *D1; delete [] D1; //Для возврата обратной матрицы из функции
}
//else fprintf (f2, "Определитель равен 0");
 
    return 0;
}