def matrix_from_expressions(row_exprs, column_exprs):
row_exprs = [fold_constants(expand(expr)) for expr in row_exprs]
result = numpy.zeros((len(row_exprs), len(column_exprs)), dtype=object)
for i, row_expr in enumerate(row_exprs):
for j, col_expr in enumerate(column_exprs):
result[i, j] = differentiate(row_expr, col_expr)
return result
def test_expand():
from pymbolic import var, expand
x = var("x")
u = (x + 1)**5
expand(u)
def test_expand():
from pymbolic import var, expand
x = var("x")
u = (x+1)**5
expand(u)
def s2s_rhs(t, yf, ys):
return fold_constants(expand(mat[1, 1] * ys()))
def f2s_rhs(t, yf, ys):
return fold_constants(expand(mat[1, 0] * yf()))
