I think this projection is very closely related to the Clement interpolation operator. This is a well-defined local operator and it is ab operator that is used for a lot of theoretical error-analysis purposes.
I will have a closer look and get back to you. Kent On 31 October 2013 10:42, Mikael Mortensen <[email protected]> wrote: > Dear all, > > I have been inspired by the new LocalSolver class of Garth and I'm > currently trying out ideas to compute "projections" locally without going > through a global coefficient matrix and a global linear algebra solve. I > just wanted to ask if anyone has any experience with this? If so, I would > be very interested to hear about it. I'm particularly interested in > speeding up the following: > > V = FunctionSpace(mesh, "CG", 1) > p = Function(V) > dpdx = project(p.dx(0), V) > > My first efforts on using a local averaging operator has been put in > branch > https://bitbucket.org/fenics-project/dolfin/branch/mikael/local-average-operator. > For small problems I have found it to be 4-5 times faster than the regular > project shown above, and the accuracy is vey similar. > > Best regards > > Mikael > > > > ------------ > > Dr. Mikael Mortensen > Associate Professor > Department of Mathematics > Mechanics Division > University of Oslo > Phone: +47-22855866 > Mob: +47-41407201 > Email: [email protected] > > _______________________________________________ > fenics mailing list > [email protected] > http://fenicsproject.org/mailman/listinfo/fenics > >
_______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
