← Back to team overview

dolfin team mailing list archive

Re: Evaluation of functions

 

I have experienced some strange behavior with the IntersectionDetector.
When I carry out the cell searching as the first thing after the mesh
variable has been created, the searching is successful. However, if I
solve a pde and use the eval function, or a locally created intersection
detector, some points are simply not found.
Has anyone else experienced this (using the python wrappings) ?

Kristen



> On Mon, Mar 17, 2008 at 09:19:02PM +0100, cosby@xxxxxxxxx wrote:
>> Would it work if I transfered the updated files in dolfin/function/ to
>> version 0.7.2 or are the any new changes that would break dependencies
>> on
>> these files ?
>>
>> Kristen
>
> Don't know, you'll have to try and see what happens...
>
> --
> Anders
>
>
>>
>> > On Mon, Mar 17, 2008 at 05:04:43PM +0100, Kristen Kaasbjerg wrote:
>> >> 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.
>> >
>> > Yes! My fault.
>> >
>> >
>>
>> _______________________________________________
>> DOLFIN-dev mailing list
>> DOLFIN-dev@xxxxxxxxxx
>> http://www.fenics.org/mailman/listinfo/dolfin-dev
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev
>



Follow ups

References