передача массивов в функцию

#include <iostream>
#include <fstream>
#include "mg.h"
using namespace std;
const int FINE_MESH = 0;
 
void main()
{
    double ****phi, ***phi_exact, ****f, ****x, ****y, ****z, ****aux;
    int n_nodes = 1025;
    int n_levels;
    cout << "allocate arrays" << endl;
    n_levels = allocate_arrays(&phi, &phi_exact, &f, &x, &y, &z, &aux, n_nodes);
    cout << "generate fine mesh" << endl;
    generate_fine_mesh(x[FINE_MESH], y[FINE_MESH], z[FINE_MESH], n_nodes, true);
    cin >> n_levels;
}

#include "mg.h"
#include <iostream>
using namespace std;
 
int allocate_arrays(double *****phi, double ****phi_exact, double *****f, double *****x, double *****y, double *****z, double *****aux, int n_nodes)
{
    double ****x_temp, ****y_temp, ****z_temp, ****phi_temp, ****f_temp, ****aux_temp;
    double ***phi_exact_temp;
    bool coarsen = true; int n_levels = 1; int nodes = n_nodes;
    while (coarsen) {
        if (((nodes - 1) % 2 == 0) && ((nodes - 1) / 2 + 1 >= 3)) {
            nodes = (nodes - 1) / 2 + 1;
            n_levels++;
        }
        else {
            coarsen = false;
        }
    }
 
    x_temp = new double***[n_levels];
    y_temp = new double***[n_levels];
    z_temp = new double***[n_levels];
    phi_temp = new double***[n_levels];
    f_temp = new double***[n_levels];
    aux_temp = new double***[n_levels];
    phi_exact_temp = new double**[n_levels];
    nodes = n_nodes;
    for (int i_level = 0; i_level < n_levels; i_level++)
    {
        x_temp[i_level] = new double**[nodes];
        y_temp[i_level] = new double**[nodes];
        z_temp[i_level] = new double**[nodes];
        phi_temp[i_level] = new double**[nodes];
        f_temp[i_level] = new double**[nodes];
        aux_temp[i_level] = new double**[nodes];
        for (int i = 0; i < nodes; i++)
        {
            x_temp[i_level][i] = new double*[nodes];
            y_temp[i_level][i] = new double*[nodes];
            z_temp[i_level][i] = new double*[nodes];
            phi_temp[i_level][i] = new double*[nodes];
            f_temp[i_level][i] = new double*[nodes];
            aux_temp[i_level][i] = new double*[nodes];
            for (int j = 0; j < nodes; j++) {
                x_temp[i_level][i][j] = new double[nodes];
                y_temp[i_level][i][j] = new double[nodes];
                z_temp[i_level][i][j] = new double[nodes];
                phi_temp[i_level][i][j] = new double[nodes];
                f_temp[i_level][i][j] = new double[nodes];
                aux_temp[i_level][i][j] = new double[nodes];
            }
        }
        nodes = (nodes - 1) / 2 + 1;
    }
    nodes = n_nodes;
    for (int i = 0; i < nodes; i++)
    {
        phi_exact_temp[i] = new double*[nodes];
        for (int j = 0; j < nodes; j++) {
            phi_exact_temp[i][j] = new double[nodes];
        }
    }
    **x = *x_temp;
    **y = *y_temp;
    **z = *z_temp;
    **f = *f_temp;
    **aux = *aux_temp;
    **phi_exact = *phi_exact_temp;
    return n_levels;
}
void generate_fine_mesh(double ***x, double ***y, double ***z, int n_nodes, bool visualize)
{
    for (int i = 0; i < n_nodes; i++) {
        for (int j = 0; j < n_nodes; j++) {
            for (int k = 0; k < n_nodes; k++) {
                x[i][j][k] = (double)i / (double)(n_nodes - 1);
                y[i][j][k] = (double)j / (double)(n_nodes - 1);
                z[i][j][k] = (double)k / (double)(n_nodes - 1);
            }
        }
    }
    cout << "создана сетка(" << n_nodes << "x" << n_nodes << ")";
}