fiat team mailing list archive

[Marie Rognes <meg@xxxxxxxxxxxxxxxxxxxxx>]

Marie Rognes wrote:
> First, the following FIAT code:
>
>     from FIAT.RaviartThomas import *
>     from FIAT.shapes import *
>     RT0 = RaviartThomas(TRIANGLE, 0) 
>     print RT0.dual_basis().entity_ids
>
> gives:
>
>     {0: {0: [], 1: [], 2: []}, 1: {0: 0, 1: 1, 2: 2}, 2: []}
>
> Shouldn't this be
>
>     {0: {0: [], 1: [], 2: []}, 1: {0: [0], 1: [1], 2: [2]}, 2: []} ?
>   
I guess a more sensible suggestion would be

{0: {0: [], 1: [], 2: []}, 1: {0: [0], 1: [1], 2: [2]}, 2: {0: []}} 


--
Marie