C++ 编译错误：error: invalid operands of types 'double***' and 'double' to binary 'operator/'？

Question

error: invalid operands of types 'double***' and 'double' to binary 'operator/'

源代码是这样的：
void Porous::Compute_Force(Flow_solver* Flow, Parallel_MPI* MPI_parallel) {

int X, Y, Z;
double F_epsilon;
if (processor_id != MASTER) {
    for (X = 2; X <= MPI_parallel->actual_rows_XYZ[0]+1; ++X) {
        for (Y = 2; Y <= MPI_parallel->actual_rows_XYZ[1]+1; ++Y) {
            for (Z = 2; Z <= MPI_parallel->end_XYZ2[2]+1; ++Z) {
                F_epsilon = 1.75 / sqrt( 150 * pow(Porosity[X][Y][Z],3) );
                Flow->force[X][Y][Z][0] = Porosity[X][Y][Z] * (Flow->viscosity/Permeability[X][Y][Z]
                                        - F_epsilon * Flow->velocity_magnitude[X][Y][Z]/sqrt(Permeability[X][Y][Z]) ) * Flow->velocity[X][Y][Z][0]
                                        + Porosity[X][Y][Z] * Flow->force[X][Y][Z][0];
                Flow->force[X][Y][Z][1] = Porosity[X][Y][Z] * ( Flow->viscosity/Permeability[X][Y][Z]
                                        - F_epsilon * Flow->velocity_magnitude[X][Y][Z] / sqrt(Permeability[X][Y][Z]) ) * Flow->velocity[X][Y][Z][1]
                                        + Porosity[X][Y][Z] * Flow->force[X][Y][Z][1];
                Flow->force[X][Y][Z][2] = Porosity[X][Y][Z] * ( Flow->viscosity / Permeability[X][Y][Z]
                                        - F_epsilon * Flow->velocity_magnitude[X][Y][Z] / sqrt(Permeability[X][Y][Z]) ) * Flow->velocity[X][Y][Z][2]
                                        + Porosity[X][Y][Z] * Flow->force[X][Y][Z][2];
                }
            }
        }
    }

}

Answer 1

https://blog.csdn.net/digent1/article/details/40538251