def main(): print('Differentiation:') x, y = map(Symbol, 'x y'.split()) expr = (x - acos(y))*exp(x + y) Dexpr = expr.diff(y) print(Dexpr) print(Dexpr._print_ccode())
def _get_array(): X, Y, Z = inp = array.array('d', [1, 2, 3]) args = x, y, z = se.symbols('x y z') exprs = [x+y+z, se.sin(x)*se.log(y)*se.exp(z)] ref = [X+Y+Z, math.sin(X)*math.log(Y)*math.exp(Z)] def check(arr): assert all([abs(x1-x2) < 1e-13 for x1, x2 in zip(ref, arr)]) return args, exprs, inp, check