Fix reuse of grid function in skeleton case

python/dune/perftool/pdelab/localoperator.py

@@ -158,6 +158,9 @@ def localoperator_basename(form_ident):
 
 
 def name_gridfunction_member(coeff, restriction, diffOrder=0):
+    # We reuse the grid function for volume integrals in skeleton integrals
+    if restriction == Restriction.NEGATIVE:
+        restriction = Restriction.NONE
     restr = "_n" if restriction == Restriction.POSITIVE else ""
     name = "local_gridfunction_coeff{}_diff{}{}".format(coeff.count(), diffOrder, restr)
     define_gridfunction_member(name, coeff, restriction, diffOrder)