dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Mon, Mar 17, 2008 at 04:20:35PM +0100, Martin Sandve Alnæs wrote:
> 2008/3/17, Anders Logg <logg@xxxxxxxxx>:
> > It's now possible to evaluate functions at arbitrary points inside the
> >  mesh where the function is defined. If evaluated outside, you should
> >  get an informative exception.
> 
> Great!
> 
> >  Thanks to Kristen Kaasbjerg for providing the code.
> >
> >  There's a demo in demo/function/. The C++ version seems to work fine,
> >  but the Python version is broken. The problem is that FFC JIT does not
> >  generate the code for evaluate_basis (you should get an informative
> >  exception about this if you install the very latest FFC). Switching it
> >  on is a simple fix, but we might want to avoid that.
> >
> >  I don't know how to solve this. Maybe we should always generate code
> >  for evalute_basis, but that might slow down the overall performance
> >  of the JIT compiler. Suggestions?
> 
> Maybe you could cache the generated code for each finite element somehow in FFC?

That would be nice, but it sounds a little messy... FFC generates one
file with everything so I would need to split it into several files.
Something to think about...

-- 
Anders