Перемножение матриц. Оптимальный алгоритм

for (int i = 0; i < nSize; i++)
    {
        for (int k = 0; k < nSize; k++)
        {
            for (int j = 0; j < nSize; j++)
            {
                C[i][j] += A[i][k] * B[k][j];
            }
        }
    }

for (int i = 0; i < nSize; i++)
    {
        for (int j = 0; j < nSize; j++)
        {
            for (int k = 0; k < nSize; k++)
            {
                C[i][j] += A[i][k] * B[k][j];
            }
        }
    }

#include <iostream>
#include <time.h>
 
void mult(int nSize, double** A, double** B, double** C)
{
    for (int i = 0; i < nSize; i++)
    {
        for (int j = 0; j < nSize; j++)
        {
            for (int k = 0; k < nSize; k++)
            {
                C[i][j] += A[i][k] * B[k][j];
            }
            std::cout << C[i][j] << " ";
        }
        std::cout << std::endl;
    }
}
 
void newMult(int nSize, double** A, double** B, double** C)
{
    for (int i = 0; i < nSize; i++)
    {
        for (int k = 0; k < nSize; k++)
        {
            for (int j = 0; j < nSize; j++)
            {
                C[i][j] += A[i][k] * B[k][j];
            }
        }
    }
    /*std::cout << "Vot tut ne verno:" << std::endl;
    for (int i = 0; i < nSize; i++)
    {
        for (int j = 0; j < nSize; j++)
        {
            std::cout << C[i][j] << " ";
        }
        std::cout << std::endl;
    }*/
}
 
void writeMatrix(int sSize, double** T)
{
    for (int i = 0; i < sSize; i++)
    {
        for (int j = 0; j < sSize; j++)
        {
            std::cout << T[i][j] << " ";
        }
        std::cout << std::endl;
    }
}
 
int main(int argc, char* argv[])
{
    int size = 3;
    double** oneMatrix = new double* [size];
    double** twoMatrix = new double* [size];
    double** product = new double* [size];
    
    clock_t Time, nTime;
 
 
    for (int i = 0; i < size; i++)
    {
        oneMatrix[i] = new double[size];
        twoMatrix[i] = new double[size];
        product[i] = new double[size];
 
        for (int j = 0; j < size; j++)
        {
            oneMatrix[i][j] = rand()%2;
            twoMatrix[i][j] = rand()%2;
        }
    }
    //std::cout << "oneMatrix:" << std::endl;
    //writeMatrix(size, oneMatrix);
 
    //std::cout << "twoMatrix:" << std::endl;
    //writeMatrix(size, twoMatrix);
 
    std::cout << "Product:" << std::endl;
    
    Time = clock();
    mult(size, oneMatrix, twoMatrix, product);
    std::cout << "Time = " << clock() - Time << std::endl;
 
    nTime = clock();
    newMult(size, oneMatrix, twoMatrix, product);   
    std::cout << "nTime = " << clock() - nTime << std::endl;
 
 
    std::system("pause");
}