C程序：如何优化矩阵幂运算

最近写了个C程序来计算矩阵的幂，代码如下：

 /*  矩阵幂运算，例如
        1 2        1 0
        3 4 ^ 0 => 0 1

        1 2        1 2
        3 4 ^ 1 => 3 4

        1 2        199 290
        3 4 ^ 4 => 435 634
    */
uint32_t* matrix_pow(const uint32_t* matrix, uint32_t exponent) {

    uint32_t* result = new_matrix();

    if (exponent == 0){
        result = identity_matrix();
    } else{
        result = cloned(matrix);
        for (int c = 1; c < exponent; c++){
            result = matrix_mul(result, matrix);
        }
    }

    return result;
}

其中, exponent是幂次，cloned()方法是复制矩阵，identity_matrix()返回一个单位矩阵，matrix_mul()用来计算两个计算的乘积。

不知有哪些优化的方法（可以完全改写上面的程序），可以使上面的函数效率更高，可以使用多线程，分治等。