Zero guess the RieszMap solver

jrmaddison · jrmaddison · commit ed5213acefce · 2025-09-18T14:01:50.000+01:00
diff --git a/firedrake/cofunction.py b/firedrake/cofunction.py
@@ -485,6 +485,7 @@ def __call__(self, value):
             else:
                 solve, rhs, soln = self._solver
                 rhs.assign(value)
+                soln.zero()
                 solve()
                 output.assign(soln)
         elif ufl.duals.is_primal(value):
