fiat team mailing list archive

["Robert Kirby" <robert.c.kirby@xxxxxxxxx>]

Well, FIAT really does simplices because they're hard.  Hexes and rectangles
have tensor-product bases for H(div), H(curl), and H1, so you should really
be doing spectral-type techniques.  See the thesis of Heintsch (sp?) at NYU
for some interesting work on high order H(div) and H(curl) a few years ago.
 Not sure if it got published, but it was very good
I now have a user at Sandia who wants the other shapes, so I will be looking
at doing this, either in Python or else with straight tensor-products in
C++.  Of course, pyramids and prisms present wrinkles of being somewhere
in-between.



On Jan 23, 2008 9:07 PM, Harish Narayanan <harish@xxxxxxxxxxx> wrote:

> Greetings,
>
> I understand from the documentation in FIAT that all its reference
> elements treat spaces as simplicial complexes. Is there an obvious
> difficulty in extending this to non-simplicial elements (assuming the
> mesh is suitably generated elsewhere)?
>
> I would be curious to read some literature on this topic if someone were
> willing to point me to it.
>
> Thanks,
> Harish
> _______________________________________________
> FIAT-dev mailing list
> FIAT-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/fiat-dev
>