Anders Logg wrote: > 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. > > 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? > > Please test it and see if it works. Would be nice if Dag could run > the benchmarks again to confirm that I didn't mess anything up. > > Note that Dag Lindbo also contributed to the development of the code to the eval function ! So thanks to him as weel.
Kristen _______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
