dolfin team mailing list archive

["Martin Sandve Alnæs" <martinal@xxxxxxxxx>]

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?

-- 
Martin