Skip to content
Snippets Groups Projects
Commit d5a9cbaa authored by René Heß's avatar René Heß
Browse files

[bugfix] in jacobian accterm splitting

parent 32212687
No related branches found
No related tags found
No related merge requests found
......@@ -153,10 +153,10 @@ def split_into_accumulation_terms(expr):
if ma.restriction != restriction else ma.expr
for ma in indexed_jac_args}
jac_expr = replace_expression(jac_expr, replacemap=replacement)
jac_accum_expr = replace_expression(jac_expr, replacemap=replacement)
if not isinstance(jac_expr, Zero):
ret.append(AccumulationTerm(jac_expr, test_arg, indexmap, newi))
ret.append(AccumulationTerm(jac_accum_expr, test_arg, indexmap, newi))
else:
if not isinstance(replace_expr, Zero):
ret.append(AccumulationTerm(replace_expr, test_arg, indexmap, newi))
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment