Вопрос по работе с "CALLOC()"

void SomeFunc(double **Arg, int R, int C)
{
    for(int i=0;i<R;i++) {
         for(int j=0;j<C;j++) {
               // Че-то-там
         }
    }
}

void MatrixAlloc(double **T,int R,int C)
{
    T = (double**) calloc(R, sizeof(double));
    for(int i = 0; i < R; i++){
      T[i] = (double*) calloc(C, sizeof(double));
      for(int j = 0; j < C; j++) T[i][j] = 0;
    }
}

__fastcall TKalman::TKalman()
{
//  double F[16][16],dF[16][16],Xmk[16],Xk_1[16],Xk[16]/*,Zk[8],Uk_1[3]*/,D[8][8],dD[8][8],C[8][8],dC[8][8],
//         R[16][3],dR[16][3],Rpi[3][16],Pk[16][16],Pmk[16][16],Pk_1[16][16],K[16][8],M[8][16],Mt[16][8],Uk_1[3],Uk[3];
    F = MatrixAlloc(16,16);
    dF = MatrixAlloc(16,16);
    FXk_1 = MatrixAlloc(16,1);
    Xmk = MatrixAlloc(16,1);
    Xk_1 = MatrixAlloc(16,1);
    Xk = MatrixAlloc(16,1);
    Zk = MatrixAlloc(16,1);
    D = MatrixAlloc(8,8);
    dD = MatrixAlloc(8,8);
    C = MatrixAlloc(8,8);
    dC = MatrixAlloc(8,8);
    R = MatrixAlloc(16,3);
    dR = MatrixAlloc(16,3);
    Rpi = MatrixAlloc(3,16);
    Pk = MatrixAlloc(16,16);
    Pmk = MatrixAlloc(16,16);
    Pk_1 = MatrixAlloc(16,16);
    K = MatrixAlloc(16,8);
    M = MatrixAlloc(8,16);
    Mt = MatrixAlloc(16,8);
    dM = MatrixAlloc(8,16);
    Uk = MatrixAlloc(3,1);
    Uk_1 = MatrixAlloc(3,1);
}
__fastcall TKalman::~TKalman()
{
//  double F[16][16],dF[16][16],Xmk[16],Xk_1[16],Xk[16]/*,Zk[8],Uk_1[3]*/,D[8][8],dD[8][8],C[8][8],dC[8][8],
//         R[16][3],dR[16][3],Rpi[3][16],Pk[16][16],Pmk[16][16],Pk_1[16][16],K[16][8],M[8][16],Mt[16][8],Uk_1[3],Uk[3];
    free(F);
    free(dF);
    free(FXk_1);
    free(Xmk);
    free(Xk_1);
    free(Xk);
    free(Zk);
    free(D);
    free(dD);
    free(C);
    free(dC);
    free(R);
    free(dR);
    free(Rpi);
    free(Pk);
    free(Pmk);
    free(Pk_1);
    free(K);
    free(M);
    free(Mt);
    free(dM);
    free(Uk);
    free(Uk_1);
}

double **TKalman::MatrixAlloc(int R,int C)
{
    double **T = (double**) calloc(R, sizeof(double));
    for(int i = 0; i < R; i++){
      T[i] = (double*) calloc(C, sizeof(double));
      for(int j = 0; j < C; j++) T[i][j] = 0;
    }
}