Powerful mathematics.
Real-world solutions.

static Matrix Strass(Matrix A, Matrix B)
{
    if (A.Cols != B.Rows)
        throw new Exception("Matrices are not conformable for multiplication");
    if (A.Cols <= 2)
        return A * B;
    else
    {
        // get matrix size
        int N = A.Cols / 2;
        // Step 1: Divide matrices into quadrants
        Matrix A11 = A[..N, ..N], A12 = A[..N, N..],
               A21 = A[N.., ..N], A22 = A[N.., N..],
​
               B11 = B[..N, ..N], B12 = B[..N, N..],
               B21 = B[N.., ..N], B22 = B[N.., N..],
​
        // Step 2: Calculate the 7 Strassen products (M1 through M7)
        M1 = Strass(A11 + A22, B11 + B22),
        M2 = Strass(A21 + A22, B11),
        M3 = Strass(A11, B12 - B22),
        M4 = Strass(A22, B21 - B11),
        M5 = Strass(A11 + A12, B22),
        M6 = Strass(A21 - A11, B11 + B12),
        M7 = Strass(A12 - A22, B21 + B22),
​
        // Step 3: Combine products into the quadrants of C
        C11 = M1 + M4 - M5 + M7,
        C12 = M3 + M5,
        C21 = M2 + M4,
        C22 = M1 - M2 + M3 + M6,
​
        // Step 4: Assemble the final matrix
        C = new Matrix[,] 
        {
            { C11, C12 }, 
            { C21, C22 } 
        };
        return C;
    }
}
​
Matrix A = Rand(8, 8), B = Rand(8, 8), C = Strass(A, B), D = A * B;
Console.WriteLine($"A = \n{A}");
Console.WriteLine($"B = \n{B}");
Console.WriteLine($"C = \n{C}");
Console.WriteLine($"D = \n{D}");