dolfin team mailing list archive

[Kristian Ølgaard <k.b.oelgaard@xxxxxxxxx>]

On 22 November 2011 21:07, Anders Logg <logg@xxxxxxxxx> wrote:
> On Tue, Nov 22, 2011 at 07:34:38PM +0000, Garth N. Wells wrote:
>> On 22 November 2011 19:12, Kristian Ølgaard <k.b.oelgaard@xxxxxxxxx> wrote:
>> > Hi,
>> >
>> > The IntervalCell::facet_area in DOLFIN currently returns 0.0, while
>> > the facet determinant (a scale factor) 'det = 1.0;'
>> > in the generated code for the tabulate_tensor function in 1D (facet integrals).
>> >
>> > Is 0.0 the correct return value, or should it return 1.0?
>> >
>>
>> I would think 1.0.

That's what I'm leaning towards too.

> Facet is the thing of dimension one lower than the cell dimension so
> in this case 1 - 1 so it would mean the area of an endpoint. So I
> believe 0 is correct.

I see your 'point' :), but is it consistent with how we evaluate
integrals on endpoints?

Kristian

> --
> Anders
>
>
>> Garth
>>
>> > Kristian
>> >
>> > _______________________________________________
>> > Mailing list: https://launchpad.net/~dolfin
>> > Post to     : dolfin@xxxxxxxxxxxxxxxxxxx
>> > Unsubscribe : https://launchpad.net/~dolfin
>> > More help   : https://help.launchpad.net/ListHelp
>> >
>>
>>
>>
>