Перемножение матриц (2 потока) help

for (i=0;i<N;i++)
    for (j=0;j<N;j++)
{
    AA_[i*N+j]=0;
    for (int k=0;k<N;k++)
       AA_[i*N+j]+=A[k*N+i]*A[k*N+j];
}

int N=1000;
    float *fA=new float[N*N]; // первая матрица
    float *fB=new float[N*N]; // вторая
    float *fC=new float[N*N]; //fC=fA*fB   
    float *fC__=new float[N*N]; // fC__=fA*fB - через потоки

unsigned int MyThread(LPVOID pParam)
{
    int *ind = (int*) pParam;
    int i=ind[0];
    int j=ind[1];

    fC__[i*N+j]=0;
    for (int k=0;k<N;k++)
        fC__[i*N+j]+=fA[i*N+k]*fB[k*N+j];

    return 0;
}

void CThreadDlg::OnOK()
{
    int i,j,k;
    int NThread=100; // количество потоков
    static int **ind;
             CWinThread **pThread=new CWinThread*[NThread];
    HANDLE *hThread=new HANDLE[NThread];
    float MinA=1000,MinB=10000;
    ind=new int*[NThread];
    for (i=0;i<NThread;i++)
    {
        ind=new int[2];
    }
           

    srand( (unsigned)time( NULL ) );
// Задаю матрицы
    for (i=0;i<N;i++)
    {
        for (j=0;j<N;j++)
        {
            fA[i*N+j]=fabs(rand()/(RAND_MAX+0.0)*(i-j)/(i+j+1.0));
            fB[i*N+j]=fabs(rand()/(RAND_MAX+0.0)*(i-j)/(i+j+1.0));
        }
    }

    for (i=0;i<N;i++)
    {
        for (j=0;j<N;j++)
        {
            fC[i*N+j]=0;
            for (k=0;k<N;k++)
                fC[i*N+j]+=fA[i*N+k]*fB[k*N+j];
        }
    }      

    for (i=0;i<N;i++)
    {
        for (j=0;j<N;)
        {
            for (k=0;k<NThread;k++)
                if (j<N)
                {
                    j++;
                    ind[k][0]=i;
                    ind[k][1]=j;
                    pThread[k]=AfxBeginThread(MyThread,(LPVOID)ind[k]);
                }
            for (k=0;k<NThread;k++)
            {
                hThread[k]= pThread[k]->m_hThread;
//              ::WaitForSingleObject(hThread[k],INFINITE);
            }
        }
    }

}