Anders Logg wrote: > No, it has not been implemented yet but it's something we desperately > need. If we define the L2 projection Pu of u onto V by > > (Pu, v) = (u, v) for all v in V > > then the problem is in evaluating the integral on the right-hand side: > > (u, v) > > where now u and v are defined on different meshes. > > My suggestion would be to evaluate this by quadrature at run-time, > using the generated UFC function > > evaluate_basis > > and IntersectionDetector/GTS which is very fast at finding the overlap > between cells.
Hmm, ok. I will try to compute the integral on the right-hands at runtime using quadrature and use it to determine Pu. _______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
