제곱 및 큐브 항을 곱셈으로 변환

Question

sympy에서 ccode()를 사용하여 sage의 유효한 식을 유효한 C 코드로 변환하려고합니다. 그러나 제 표현에는 많은 제곱 및 큐브 용어가 있습니다. pow (x, 2)가 x * x보다 훨씬 느리기 때문에 변환하기 전에 표현식에있는 용어를 확장하려고합니다.

from sympy import Symbol, Mul, Pow, pprint, Matrix, symbols 
    from sympy.core import numbers 


    def pow_to_mul(expr): 
     """ 
     Convert integer powers in an expression to Muls, like a**2 => a*a. 
     """ 
     pows = list(expr.atoms(Pow)) 
     pows = [p for p in pows if p.as_base_exp()[1]>=0] 
     if any(not e.is_Integer for b, e in (i.as_base_exp() for i in pows)): 

       raise ValueError("A power contains a non-integer exponent") 
     repl = zip(pows, (Mul(*[b]*e,evaluate=False) for b,e in (i.as_base_exp() for i in pows))) 
     return expr.subs(repl) 

그것은 부분적으로 작동하지만, 한 전원이 곱셈의 인수와 같이 실패 :

>>>_=var('x') 
    >>>print pow_to_mul((x^3+2*x^2)._sympy_()) 
    2*x**2 + x*x*x 
    >>>print pow_to_mul((x^2/(1+x^2)+(1-x^2)/(1+x^2))._sympy_()) 
    x**2/(x*x + 1) - (x*x - 1)/(x*x + 1) 

왜 는 this 대화를 바탕으로, 나는 다음과 같은 코드를 작성? 어떻게 변경할 수 있습니까? 고맙습니다.

Answer 1

-ffast-math으로 컴파일하면 컴파일러에서이 최적화 작업을 수행합니다. 정확히 실패 무엇

>>> ccode(x**97 + 4*x**7 + 5*x**3 + 3**pi, user_functions={'Pow': [ 
... (lambda b, e: e.is_Integer and e < 42, lambda b, e: '*'.join([b]*int(e))), 
... (lambda b, e: not e.is_Integer, 'pow')]}) 
'pow(x, 97) + 4*x*x*x*x*x*x*x + 5*x*x*x + pow(3, M_PI)'