Sympy Codegen для MatrixSymbol Resp. сумма произвольного числа элементов ⇐ C++

[Anonymous]

Я попытался генерировать код C для матриц произвольного размера с помощью Sympy, но код недействителен. Корльная причина, по -видимому, заключается в том, что функция «суммы» симпей не является должным образом переводится.
Вот пример: < /p>
from sympy import symbols, Sum
from sympy.utilities.codegen import codegen

n = symbols('n', integer=True)
i = symbols('i')
my_sum = Sum(i, (i,0,n))
pprint(codegen(("f", my_sum), "C", "test", header=False, empty=False)[0][1])
< /code>
output: < /p>
#include "test.h"
#include
double f(int n) {
double f_result;
f_result = Sum(i, (i, 0, n));
return f_result;
}
< /code>
Есть ли способ генерировать правильный код C для этого простого выражения? < /p>
Правильный код C ++ может прочитать < /p>
#include

double f(int n) {
double f_result;

auto sympy_sum_expr [](int i)->i {return i;};
auto sympy_sum = [](std::function expr, int from, int to)->double
{
double sum = 0.0;
std::foreach([&](int i)->void {sum+=expr(i);}, std::iota_view{from,to+1}
return sum;
};

f_result = sympy_sum(sympy_sum_expr, 0, n);
return f_result;
}


Подробнее здесь: https://stackoverflow.com/questions/794 ... f-elements