Re: [Libmesh-devel] Ghosted vectors

2009-03-02 Thread Roy Stogner
On Mon, 2 Mar 2009, Tim Kroeger wrote: Update: It doesn't crash any more. With the patch that I sent you in the previous mail, it seems to work. That's interesting. Did you ever track down the source of the problem? The patch you provided looks like it could have fixed some inaccuracy bugs

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] Ghosted vectors

2009-03-02 Thread Tim Kroeger
Dear Roy, On Wed, 25 Feb 2009, Tim Kroeger wrote: > On Wed, 25 Feb 2009, Tim Kroeger wrote: > >> I'm currently trying to test the ghosted vectors on my application, >> but the queue on the cluster is full, so I have to wait. > > Update: My application crashes with the ghosted vectors, but I'm not

Re: [Libmesh-devel] Ghosted vectors

2009-03-02 Thread Tim Kroeger
Dear Roy, On Wed, 25 Feb 2009, Tim Kroeger wrote: On Wed, 25 Feb 2009, Tim Kroeger wrote: Besides, this can't be the right fix, can it? Why should any parallel communication be required to copy one consistent ghosted vector to another? It seems as if VecCopy is simply failing to copy the wh

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