Re: [Libmesh-users] Gradient of C1 error

2012-12-06 Thread Ataollah Mesgarnejad
Ok done. You are right, when I write the output in Nemesis format it messes up the gradient. You can get my tweaked adaptivity_ex4 from: http://cl.ly/02131w2o2715 . So does this mean that this is a VisIt bug?! PS: I had to make it so it reads a mesh file otherwise Nemesis_IO wouldn't work and I

Re: [Libmesh-users] Gradient of C1 error

2012-12-04 Thread Derek Gaston
Ata, Are you sure you need Nemesis? This should really only be necessary once you get up over ~200 million DoFs or so. Anything smaller than that should work fine with Exodus. Not saying you haven't found a bug here... just trying to make sure you're not overly burdening yourself with extra wor

Re: [Libmesh-users] Gradient of C1 error

2012-12-04 Thread Ataollah Mesgarnejad
On Dec 4, 2012, at 4:53 PM, Roy Stogner wrote: > > On Tue, 4 Dec 2012, Ataollah Mesgarnejad wrote: > >> Ok done. You are right, when I write the output in Nemesis format it >> messes up the gradient. You can get my tweaked adaptivity_ex4 from: >> http://cl.ly/02131w2o2715 . So does this mean t

Re: [Libmesh-users] Gradient of C1 error

2012-12-04 Thread Roy Stogner
On Tue, 4 Dec 2012, Ataollah Mesgarnejad wrote: > Ok done. You are right, when I write the output in Nemesis format it > messes up the gradient. You can get my tweaked adaptivity_ex4 from: > http://cl.ly/02131w2o2715 . So does this mean that this is a VisIt > bug?! Possible, but not necessarily.

Re: [Libmesh-users] Gradient of C1 error

2012-12-04 Thread Roy Stogner
On Tue, 4 Dec 2012, Ataollah Mesgarnejad wrote: > I see what you mean. And no I didn't I did the gradient calculation in VisIt. > I don't use/calculate this gradient in my code. This is part of your (our) problem, I'm afraid: we don't output any shape functions more interesting than C0 quadrati

Re: [Libmesh-users] Gradient of C1 error

2012-12-04 Thread Ataollah Mesgarnejad
Roy, 1st as always thanks for the quick response. I see what you mean. And no I didn't I did the gradient calculation in VisIt. I don't use/calculate this gradient in my code. Here is the reason I suspected gradients: the problem I'm solving is complicated but in short what I solve for is a

Re: [Libmesh-users] Gradient of C1 error

2012-12-04 Thread Roy Stogner
On Tue, 4 Dec 2012, Roy Stogner wrote: > On Tue, 4 Dec 2012, Ataollah Mesgarnejad wrote: > >> I found a potential bug with Clough elements. I'm solving a 4th order >> problem using C1 Clough elements and when I plot magnitude of gradient of >> the solution I found out that it is discontinuous on

Re: [Libmesh-users] Gradient of C1 error

2012-12-04 Thread Roy Stogner
On Tue, 4 Dec 2012, Ataollah Mesgarnejad wrote: > I found a potential bug with Clough elements. I'm solving a 4th order > problem using C1 Clough elements and when I plot magnitude of gradient of > the solution I found out that it is discontinuous on mesh partition > boundaries. I would appreciat

[Libmesh-users] Gradient of C1 error

2012-12-04 Thread Ataollah Mesgarnejad
Dear all, I found a potential bug with Clough elements. I'm solving a 4th order problem using C1 Clough elements and when I plot magnitude of gradient of the solution I found out that it is discontinuous on mesh partition boundaries. I would appreciate if any comments. Thanks, Ata ---