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
