[Ffc] evaluate_integrand

2010-04-12 Thread Anders Logg
We are doing some work where we need to do run-time quadrature over arbitrary polyhedra. For this reason, we would like to generate code for evaluating the integrand of a form at an arbitrary point within the cell. I would propose adding something like this to each of the *_integral classes in the

Re: [Ffc] evaluate_integrand

2010-04-12 Thread Garth N. Wells
On 12/04/10 20:47, Anders Logg wrote: We are doing some work where we need to do run-time quadrature over arbitrary polyhedra. We (Mehdi and I) do this already (using UFC), so I don't see why a new function is required. Can you explain why evaluate_tensor is not enough? Garth For this re

Re: [Ffc] evaluate_integrand

2010-04-12 Thread Garth N. Wells
On 12/04/10 21:19, Garth N. Wells wrote: On 12/04/10 20:47, Anders Logg wrote: We are doing some work where we need to do run-time quadrature over arbitrary polyhedra. We (Mehdi and I) do this already (using UFC), so I don't see why a new function is required. Can you explain why evaluate_

Re: [Ffc] evaluate_integrand

2010-04-12 Thread Anders Logg
On Mon, Apr 12, 2010 at 09:21:32PM +0800, Garth N. Wells wrote: > > > On 12/04/10 21:19, Garth N. Wells wrote: > > > > > >On 12/04/10 20:47, Anders Logg wrote: > >>We are doing some work where we need to do run-time quadrature over > >>arbitrary polyhedra. > > > >We (Mehdi and I) do this already (u

Re: [Ffc] evaluate_integrand

2010-04-12 Thread Garth N. Wells
On 12/04/10 21:29, Anders Logg wrote: On Mon, Apr 12, 2010 at 09:21:32PM +0800, Garth N. Wells wrote: On 12/04/10 21:19, Garth N. Wells wrote: On 12/04/10 20:47, Anders Logg wrote: We are doing some work where we need to do run-time quadrature over arbitrary polyhedra. We (Mehdi and I)

Re: [Ffc] evaluate_integrand

2010-04-12 Thread Anders Logg
On Mon, Apr 12, 2010 at 09:34:38PM +0800, Garth N. Wells wrote: > > > On 12/04/10 21:29, Anders Logg wrote: > >On Mon, Apr 12, 2010 at 09:21:32PM +0800, Garth N. Wells wrote: > >> > >> > >>On 12/04/10 21:19, Garth N. Wells wrote: > >>> > >>> > >>>On 12/04/10 20:47, Anders Logg wrote: > We are d

Re: [Ffc] evaluate_integrand

2010-04-12 Thread Garth N. Wells
On 12/04/10 21:49, Anders Logg wrote: On Mon, Apr 12, 2010 at 09:34:38PM +0800, Garth N. Wells wrote: On 12/04/10 21:29, Anders Logg wrote: On Mon, Apr 12, 2010 at 09:21:32PM +0800, Garth N. Wells wrote: On 12/04/10 21:19, Garth N. Wells wrote: On 12/04/10 20:47, Anders Logg wrote: W

Re: [Ffc] Bubble function

2010-04-12 Thread Andy Ray Terrel
On Mon, Apr 12, 2010 at 1:51 AM, Marie Rognes wrote: > Andy Ray Terrel wrote: >> >> Okay so Brezzi Fortin have just a bubble function with one node. >> Perhaps they meant P3 but their notation is ambiguous. >> ArnoldBrezziFortin1984 does call for the P3 in the original paper but >> then they menti

Re: [Ffc] evaluate_integrand

2010-04-12 Thread Anders Logg
On Mon, Apr 12, 2010 at 10:20:13PM +0800, Garth N. Wells wrote: > > > On 12/04/10 21:49, Anders Logg wrote: > >On Mon, Apr 12, 2010 at 09:34:38PM +0800, Garth N. Wells wrote: > >> > >> > >>On 12/04/10 21:29, Anders Logg wrote: > >>>On Mon, Apr 12, 2010 at 09:21:32PM +0800, Garth N. Wells wrote: > >

Re: [Ffc] evaluate_integrand

2010-04-12 Thread Andre Massing
Hi! Anders and I are working on an implementation for Nitsche's method, where the question of integration on cut cell and facets also (beside XFEM for you) came up. On Monday 12. April 2010 16.20.13 Garth N. Wells wrote: > On 12/04/10 21:49, Anders Logg wrote: > > On Mon, Apr 12, 2010 at 09:34

Re: [Ffc] evaluate_integrand

2010-04-12 Thread Garth N. Wells
On 12/04/10 23:35, Anders Logg wrote: On Mon, Apr 12, 2010 at 10:20:13PM +0800, Garth N. Wells wrote: On 12/04/10 21:49, Anders Logg wrote: On Mon, Apr 12, 2010 at 09:34:38PM +0800, Garth N. Wells wrote: On 12/04/10 21:29, Anders Logg wrote: On Mon, Apr 12, 2010 at 09:21:32PM +0800, Gar

Re: [Ffc] evaluate_integrand

2010-04-12 Thread Garth N. Wells
On 13/04/10 00:16, Andre Massing wrote: Hi! Anders and I are working on an implementation for Nitsche's method, where the question of integration on cut cell and facets also (beside XFEM for you) came up. On Monday 12. April 2010 16.20.13 Garth N. Wells wrote: On 12/04/10 21:49, Anders Logg