I saw under http://docs.sympy.org/dev/tutorial/simplification.html#powsimp that it is impossible to combine radicals using powersimp:
"This means that it will be impossible to undo this identity with powsimp(), because even if powsimp() were to put the bases together, they would be automatically split apart again." I was wondering if it was possible to do this any other way. For a toy example I have import sympy L = sympy.symbols('L', real=True, finite=True, positive=True) sympy.sqrt(L) * sympy.sqrt(pi) and I would like to have it return sympy.sqrt(L * pi) Is there any way to do this? What I'd really like is if it combined these terms in this real example: import simplify import vtool as vt import sympy sigma, dist, L = sympy.symbols('sigma, distij, L', real=True, finite=True, positive=True) kernel = (1 / sympy.sqrt(sigma ** 2 * 2 * sympy.pi)) * sympy.exp((-dist ** 2) / (2 * sigma ** 2)) phi = (1 / L) * kernel logphi = sympy.simplify(sympy.log(phi)) logphi = sympy.logcombine(logphi) So I would get -distij**2/(2*sigma**2) - log(sqrt(2 * pi)*L*sigma) instead of -distij**2/(2*sigma**2) - log(sqrt(2)*sqrt(pi)*L*sigma) -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving emails from it, send an email to sympy+unsubscr...@googlegroups.com. To post to this group, send email to sympy@googlegroups.com. Visit this group at https://groups.google.com/group/sympy. To view this discussion on the web visit https://groups.google.com/d/msgid/sympy/c2d73b5e-60d0-4140-af8e-033ec7234890%40googlegroups.com. For more options, visit https://groups.google.com/d/optout.