Hello, I am working on developing a method where optimal discontinous test functions are computed for the global problem. The issue is that these optimal test function computations are decoupled from the global problem and are over each element individually. Has this been done or is there a way with par_loop() to perform a solve over each element of a mesh and have its solutions contribute to the global solve at the end? Each element problem has its own variational statement with its own test functions which we would define. Using these results we would construct the global stiffness matrix(function) as we solve this optimal test function problem over each element. We have implemented a very simple version of this with an in house code but your framework allows for so much flexibility we would love to be able to do this in Firedrake. Regards, -- Austin Kaul