Matrix Multiply

Naive algorithm: triple nested loop, O(n³). For each element C[i,j], compute the dot product of row i of A and column j of B.

Cilia

func multiply(Matrix A, B) -> Matrix {
    const m = A.rows()
    const n = A.columns()
    const o = B.rows()
    const p = B.columns()

    assert(n == o, "A.columns() must equal B.rows()")

    Matrix C(m, p)

    for i in 0..<m {
        for k in 0..<n {
            for j in 0..<p {
                C[i, j] += A[i, k] * B[k, j]
            }
        }
    }

    return C
}

C++

auto multiply(const Matrix& A, const Matrix& B) -> Matrix {
    const int m = A.rows();
    const int n = A.columns();
    const int o = B.rows();
    const int p = B.columns();

    assert(n == o && "A.columns() must equal B.rows()");

    Matrix C(m, p);

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

    return C;
}