Re: [Libmesh-devel] Higher than second order elements.

2009-03-02 Thread Kirk, Benjamin (JSC-EG311)
> wrote: >> is that these tend to be more "industrial-type" applications, where >> higher-order elements are often not used for various reasons (non-smooth >> solutions, sharp complex geometry, etc...) > > In fact, hp-fem performs the best exactly with solutions that are both > non-smooth and sha

Re: [Libmesh-devel] Higher than second order elements.

2009-03-02 Thread Ondrej Certik
Hi, I CCed to our hpfem list, so that people in our group can have a look at this too. On Mon, Mar 2, 2009 at 5:35 PM, Roy Stogner wrote: > > > On Mon, 2 Mar 2009, Ondrej Certik wrote: > >> On Mon, Mar 2, 2009 at 12:16 PM, Roy Stogner >> wrote: >>> >>> >>> On Mon, 2 Mar 2009, Ondrej Certik wrot

Re: [Libmesh-devel] Higher than second order elements.

2009-03-02 Thread Roy Stogner
On Mon, 2 Mar 2009, Ondrej Certik wrote: > That's awesome that you are interested in this. At this moment, I am > quite busy with other things, but I am very interested in this > comparison. Does the isotropic hp-adaptivity work in 3d? It works for any elements we've got hierarchic basis functio

Re: [Libmesh-devel] Higher than second order elements.

2009-03-02 Thread Roy Stogner
On Mon, 2 Mar 2009, Ondrej Certik wrote: On Mon, Mar 2, 2009 at 12:16 PM, Roy Stogner wrote: On Mon, 2 Mar 2009, Ondrej Certik wrote: If I find time, I'd like to create Python interface to libmesh, because I used it in the past and I would like to compare libmesh h-adaptivity with our hp

Re: [Libmesh-devel] Higher than second order elements.

2009-03-02 Thread Ondrej Certik
On Mon, Mar 2, 2009 at 12:16 PM, Roy Stogner wrote: > > > On Mon, 2 Mar 2009, Ondrej Certik wrote: > >> If I find time, I'd like to create Python interface to libmesh, >> because I used it in the past and I would like to compare libmesh >> h-adaptivity with our hp-adaptivity. > > If you're already

Re: [Libmesh-devel] Higher than second order elements.

2009-03-02 Thread Roy Stogner
On Mon, 2 Mar 2009, Ondrej Certik wrote: > If I find time, I'd like to create Python interface to libmesh, > because I used it in the past and I would like to compare libmesh > h-adaptivity with our hp-adaptivity. If you're already doing that, you might want to go one step farther and compare l

Re: [Libmesh-devel] Higher than second order elements.

2009-03-02 Thread Ondrej Certik
On Sun, Mar 1, 2009 at 11:31 AM, Roy Stogner wrote: > > On Sun, 1 Mar 2009, Ondrej Certik wrote: > >> On Tue, Feb 24, 2009 at 12:26 PM, Roy Stogner >> wrote: >>> >>> On Tue, 24 Feb 2009, Vijay S. Mahadevan wrote: >>> I understand that all the developers are doing your PhDs >>> >>> It's gotte

Re: [Libmesh-devel] Higher than second order elements.

2009-03-01 Thread Kirk, Benjamin (JSC-EG311)
>>> I understand that all the developers are doing your PhDs >> >> It's gotten worse than that now that the primary developers are all >> graduated.  E.g. the adaptive p refinement was something I did on a >> lark as a break from my dissertation research.  It's harder to find >> time for things li

Re: [Libmesh-devel] Higher than second order elements.

2009-03-01 Thread Roy Stogner
On Sun, 1 Mar 2009, Ondrej Certik wrote: On Tue, Feb 24, 2009 at 12:26 PM, Roy Stogner wrote: On Tue, 24 Feb 2009, Vijay S. Mahadevan wrote: I understand that all the developers are doing your PhDs It's gotten worse than that now that the primary developers are all graduated.  E.g. the a

Re: [Libmesh-devel] Higher than second order elements.

2009-03-01 Thread Ondrej Certik
Hi Roy! On Tue, Feb 24, 2009 at 12:26 PM, Roy Stogner wrote: > > On Tue, 24 Feb 2009, Vijay S. Mahadevan wrote: > >> I understand that all the developers are doing your PhDs > > It's gotten worse than that now that the primary developers are all > graduated.  E.g. the adaptive p refinement was so

Re: [Libmesh-devel] Higher than second order elements.

2009-02-24 Thread Roy Stogner
On Tue, 24 Feb 2009, Vijay S. Mahadevan wrote: > I understand that all the developers are doing your PhDs It's gotten worse than that now that the primary developers are all graduated. E.g. the adaptive p refinement was something I did on a lark as a break from my dissertation research. It's h

Re: [Libmesh-devel] Higher than second order elements.

2009-02-24 Thread Vijay S. Mahadevan
> There are higher order elements, just not higher order Lagrange > elements; try the HIERARCHIC basis instead if you want C0 higher order > polynomials. You can even do adaptive p refinement with them. > (probably hp, too, but I wouldn't guarantee that's working) Thanks for the quick reply. I di

Re: [Libmesh-devel] Higher than second order elements.

2009-02-24 Thread Roy Stogner
On Tue, 24 Feb 2009, Vijay S. Mahadevan wrote: > Are there any higher than second order Finite Elements in Libmesh ? I > obviously do not see them in the class docs but just wondering if > there was some trick to apply say a 5th order Lagrange basis to a > QUAD4 elem. Then would Libmesh automatic

[Libmesh-devel] Higher than second order elements.

2009-02-24 Thread Vijay S. Mahadevan
Hi, Are there any higher than second order Finite Elements in Libmesh ? I obviously do not see them in the class docs but just wondering if there was some trick to apply say a 5th order Lagrange basis to a QUAD4 elem. Then would Libmesh automatically create the extra dofs needed to make this uniso