Re: [Libmesh-users] Gradient of C1 error

[Roy S. <roy...@ic...>]

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 quadratics in any formats
that viz software can read.  So basically the C1 3-split cubic
functions you get with CLOUGH in libMesh get reinterpolated onto C0
1-split quadratic functions.  I wouldn't expect the gradients to be
continous *anywhere*.

There shouldn't be a larger discontinuity than average along partition
boundaries, though.  That still might indicate a localization bug,
either in your code or ours.

If it is a problem in library code... I don't suppose there's any
chance you could reproduce it with a tweaked adaptivity_ex4 ?
---
Roy