On 10/08/12 16:15, Christian Jacobs wrote:
> Public bug reported:
...
> After changing the following line in assemble/Momentum_DG.F90:
> assemble_element = .not.dg.or.element_neighbour_owned(U, ele)
> to:
> assemble_element = .not.dg.or.element_owned(U, ele)
I was convinced at the time that element_neighbour_owned(U, ele) would
return .true. if ele was owned. i.e. it's a superset of
element_owned(U, ele). However, it's possible that for some corner
cases this may not be true.
On 10/08/12 16:15, Christian Jacobs wrote:
> Public bug reported:
...
> After changing the following line in assemble/ Momentum_ DG.F90: or.element_ neighbour_ owned(U, ele) or.element_ owned(U, ele)
> assemble_element = .not.dg.
> to:
> assemble_element = .not.dg.
I was convinced at the time that element_ neighbour_ owned(U, ele) would
return .true. if ele was owned. i.e. it's a superset of
element_owned(U, ele). However, it's possible that for some corner
cases this may not be true.
Does it work if you change the test to:
.not.dg .or. (element_owned(U, ele) .or. element_ neighbour_ owned(U, ele))
?
Lawrence