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

Avoid direct vector acces for gradvec

At the moment loopy tries to acces a vector class variable and throws
a not implemented error. Probably needs some changes in the horizontal
add code.
parent d6d3c69c
No related branches found
No related tags found
No related merge requests found
......@@ -245,7 +245,11 @@ def generate_accumulation_instruction(visitor, accterm, measure, subdomain_id):
)
# In the case of FastDGGridOperator we can write directly into the resiudal/jacobi
if get_option('fastdg'):
#
# TODO: At the moment this only works if we do not vectorize
# (over gradients) because loopy tries to acces a vectorclass
# variable.
if get_option('fastdg') and not index:
ft = get_global_context_value("form_type")
if ft=='residual':
accum = accum + ".data()"
......@@ -255,7 +259,7 @@ def generate_accumulation_instruction(visitor, accterm, measure, subdomain_id):
expression = prim.Sum((assignee,result))
instruction(assignee=assignee,
expression=expression,
forced_iname_deps=frozenset(inames + visitor.inames + vecinames),
forced_iname_deps=frozenset(inames),
forced_iname_deps_is_final=True,
depends_on=insn_dep,
)
......@@ -268,7 +272,7 @@ def generate_accumulation_instruction(visitor, accterm, measure, subdomain_id):
expression = prim.Sum((assignee,result))
instruction(assignee=assignee,
expression=expression,
forced_iname_deps=frozenset(inames + visitor.inames + vecinames),
forced_iname_deps=frozenset(inames + visitor.inames),
forced_iname_deps_is_final=True,
depends_on=insn_dep,
)
......
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