#pragma omp parallel for
        for(i=0;i<n;i++){
            for(k=0;k< (i+1);k++){
                // int lim1 = 0;
                // if(lim1<(k-n/16 +1)){
                //     lim1 = k-n/16 +1;
                // }
                int lim1 = MAX(0,k-n/16+1);
                // int lim2 = n;
                // if(lim2>(k+n/16)){
                //     lim2 = k+n/16;
                // }
                int lim2 = MIN(n, k+n/16);
                for(j=lim1;j<lim2;j++){
                    A[n*i+j] += B[n*i+k]*C[n*k+j];
                }
            }
        }