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 ------------------------------------------------------------------------------ LogMeIn Rescue: Anywhere, Anytime Remote support for IT. Free Trial Remotely access PCs and mobile devices and provide instant support Improve your efficiency, and focus on delivering more value-add services Discover what IT Professionals Know. Rescue delivers http://p.sf.net/sfu/logmein_12329d2d _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users