dolfin team mailing list archive

[Kristian Oelgaard <k.b.oelgaard@xxxxxxxxxx>]

Quoting Anders Logg <logg@xxxxxxxxx>:

> On Wed, Feb 13, 2008 at 10:43:40PM +0100, Kristian Oelgaard wrote:
> > Quoting Dag Lindbo <dag@xxxxxxxxxx>:
> > 
> > > >
> > > > 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
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev
>