将 Python 风格的矩阵转换为 C 风格的矩阵

Question

我正在使用 Python(x,y) / Spyder / Sympy 进行一些符号矩阵计算，然后我想将其写入我的 C 程序中。

是否有一种工具（Python 与否）可以将这个以 Python 风格编写的结果矩阵作为一行或多行（有效）C 代码输出？

例如如果我定义

A = [[a, 0] [b, 0]]
B = [[c, d] [0, 0]]

然后让Python计算并打印A*B然后我会得到

Matrix([[a*c, a*d], [b*c d*c]])

然后我想将其转换为：

result[0][0] = a*c;
result[0][1] = a*d;
result[1][0] = b*c;
result[1][1] = d*c;

我可以在 Python 中手动完成它的核心，因为基本算术运算符在 Python 和 C 中是相同的。

但是一旦我们开始使用幂和三角函数，情况就会发生一些变化......而且可能还有其他一些我现在无法想到的东西，它们会在我的矩阵中出现变得更复杂并使用更高级的标准数学函数（相信我，他们会的）。

所以一些可以为我自动化的代码或工具会很棒。

Answer 1

在关于任意维矩阵的问题之后：Arbitrary matrix，我使用np.ndenumerate()对矩阵进行迭代。

import numpy as np

arr = np.random.random([4,5,2,6])
print(arr)
c_code = ''

for idx, value in np.ndenumerate(arr):
    print(idx, value)
    index_c = ""
    for index in idx:
        #print('print index')
        #print(index)
        index_c = index_c + '[' + str(index) + ']'

    cell_content = arr[idx]

    #print('printing content')
    #print(cell_content)

    #print('print index_c')
    #print(str(index_c))

    single_entry = ''
    single_entry = 'arr'+str(index_c)+'='+ str(cell_content)

    #print('printing single entry')
    #print(single_entry)

    c_code = c_code + str(single_entry) + ';' + '\n'

print('Printing C_code')
print(c_code)

我留下了一些测试照片。

我使用arr = np.random.random([4,5,2,6])生成了一个[4,5,2,6]随机值的维度矩阵（只是为了测试，你可以使用自己的矩阵）。

矩阵的 C 代码包含在 c_code 变量中，类似于：

arr[3][3][1][2]=0.15905407213419076;
arr[3][3][1][3]=0.11357812551781832;
arr[3][3][1][4]=0.8884023323331601;

例如。