Quoting Anders Logg <[EMAIL PROTECTED]>: > On Wed, Feb 13, 2008 at 10:43:40PM +0100, Kristian Oelgaard wrote: > > Quoting Dag Lindbo <[EMAIL PROTECTED]>: > > > > > > > > > > Why do you need to evaluate functions at arbitrary points? > > > > > > > > I can understand it's nice to be able to do this, but in many cases > > > > there are other solutions (which are much more efficient). > > > > > > > > > > I think there is general interest in being able to compute line > integrals > > > over a computed solution. It is natural to use the finite elements as > > > interpolant to get values at quadrature points. Any thoughts? > > > > If you just want the values at quadrature points, you can project the > solution > > (Function) onto a QuadratureElement. > > > > Kristian > > I guess the quadrature points are given by some other entity (like an > edge) intersecting a cell so they wouldn't fit.
OK, that way around it wouldn't work. So using evaluate_basis(_all)() would definitely be the way to go. Kristian > -- > Anders > _______________________________________________ > DOLFIN-dev mailing list > [email protected] > http://www.fenics.org/mailman/listinfo/dolfin-dev > _______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
