← Back to team overview

dolfin team mailing list archive

Re: Adding quadratic elements

 

On Fri, Dec 16, 2005 at 05:22:12PM +0100, Garth N. Wells wrote:
> On Fri, 2005-12-16 at 17:20 +0100, Garth N. Wells wrote:
> > On Fri, 2005-12-16 at 10:12 -0600, Anders Logg wrote:
> > > Yes, it would require some major work in quite a few places, FIAT,
> > > FFC, DOLFIN, and it's not on my TODO list for the next couple of
> > > months at least.
> > > 
> > 
> > Is the plan to support arbitrary order isoparametric mappings, or just a
> > limited number (linear and quadratic only, for example)? 
> > 
> > Garth 
> > 
> 
> Sorry Anders, I read over that inconspicuous little word "not" in front
> of "on my TODO list".
>
> Garth

Wishful thinking... :-)

/Anders


> > 
> > > /Anders
> > > 
> > > 
> > > On Fri, Dec 16, 2005 at 05:05:23PM +0100, Luca Antiga wrote:
> > > > Hi guys, thanks for the prompt reply.
> > > > In fact, I was referring to curved elements (sorry for not being 100% clear),
> > > > but nevertheless
> > > > I found the Poisson example very informative. A really nice way to go.
> > > > So, if I understand correctly, implementing curved tets would be on the same
> > > > level as implementing
> > > > bricks. Am I right?
> > > > 
> > > > Luca
> > > > 
> > > > On 12/16/05, Anders Logg <logg@xxxxxxxxx> wrote:
> > > > 
> > > >     Right.
> > > > 
> > > >     /Anders
> > > > 
> > > >     On Fri, Dec 16, 2005 at 09:49:25AM -0600, Robert C Kirby wrote:
> > > > 
> > > >     > No -- I think higher order geometry is what is meant - curved elements.
> > > >     >  We don't do those.
> > > >     >
> > > >     > On Dec 16, 2005, at 9:47 AM, Anders Logg wrote:
> > > >     >
> > > >     > >As Garth points out, we don't support prisms (yet).
> > > >     > >
> > > >     > >For higher order elements, take a look at the example in
> > > >     > >
> > > >     > >    src/demo/fem/convergence/
> > > >     > >
> > > >     > >which solves Poisson's equation in 2D and 3D with order ranging
> > > >     > >from q = 1 to q = 5.
> > > >     > >
> > > >     > >The only thing you need to change to do higher order elements is a
> > > >     > >number
> > > >     > >(and then recompile with FFC).
> > > >     > >
> > > >     > >/Anders
> > > >     > >
> > > >     > >
> > > >     > >On Fri, Dec 16, 2005 at 10:45:07AM +0100, Luca Antiga wrote:
> > > >     > >>Dear Dolfin developers,
> > > >     > >>  I've looked into Dolfin today and it's a really impressive piece of
> > > >     > >>software.
> > > >     > >>Congratulations for your work.
> > > >     > >>Just a curiosity: how hard it would be to code different element
> > > >     > >>types (e.g.
> > > >     > >>quadratic tetrahedra or prisms) into Dolfin?
> > > >     > >>Thanks a lot
> > > >     > >>
> > > >     > >>Luca
> > > >     > >>
> > > >     > >>
> > > >     > >>
> > > >     > >
> > > >     > >_______________________________________________
> > > >     > >DOLFIN-dev mailing list
> > > >     > > DOLFIN-dev@xxxxxxxxxx
> > > >     > >http://www.fenics.org/cgi-bin/mailman/listinfo/dolfin-dev
> > > >     >
> > > >     >
> > > >     > _______________________________________________
> > > >     > DOLFIN-dev mailing list
> > > >     > DOLFIN-dev@xxxxxxxxxx
> > > >     > http://www.fenics.org/cgi-bin/mailman/listinfo/dolfin-dev
> > > >     >
> > > > 
> > > > 
> > > >     _______________________________________________
> > > >     DOLFIN-dev mailing list
> > > >     DOLFIN-dev@xxxxxxxxxx
> > > >     http://www.fenics.org/cgi-bin/mailman/listinfo/dolfin-dev
> > > > 
> > > > 
> > > 
> > > > _______________________________________________
> > > > DOLFIN-dev mailing list
> > > > DOLFIN-dev@xxxxxxxxxx
> > > > http://www.fenics.org/cgi-bin/mailman/listinfo/dolfin-dev
> > > 
> > > 

-- 
Anders Logg
Research Assistant Professor
Toyota Technological Institute at Chicago
http://www.tti-c.org/logg/



References