Re: [Libmesh-devel] Derivative evaluations on manifolds

2011-06-16 Thread Derek Gaston
Wow - we'll certainly run the full suite of tests on Monday (and maybe create a few new ones!). Thanks for catching this and fixing it! Derek On Jun 16, 2011, at 7:35 PM, Roy Stogner wrote: > > > On Thu, 16 Jun 2011, Derek Gaston wrote: > >> we use dphi and d2phi for 1D elements in 2D all t

Re: [Libmesh-devel] Derivative evaluations on manifolds

2011-06-16 Thread Roy Stogner
On Thu, 16 Jun 2011, Derek Gaston wrote: > we use dphi and d2phi for 1D elements in 2D all the time as well. > Are you sure it's wrong... Pretty sure. "dphi[i][p](1) = 0." doesn't leave a lot of room for uncertainty. > because we're definitely getting the > right answers... even for analytic

Re: [Libmesh-devel] Derivative evaluations on manifolds

2011-06-16 Thread Derek Gaston
Ok good! Because we use 2D in 3D all the time and actually now that I think about it we use dphi and d2phi for 1D elements in 2D all the time as well. Are you sure it's wrong... because we're definitely getting the right answers... even for analytic solutions. I don't have time to check f

Re: [Libmesh-devel] Derivative evaluations on manifolds

2011-06-16 Thread Roy Stogner
On Thu, 16 Jun 2011, Roy Stogner wrote: > Does libMesh not calculate dphi and d2phi properly for > lower-dimensional elements embedded at odd angles in > higher-dimensional spaces? Things might not be as bad as I thought - it looks like the 2D-embedded-in-3D code was correct, just the 1D-embedde

[Libmesh-devel] Derivative evaluations on manifolds

2011-06-16 Thread Roy Stogner
So here's an embarrassing question from me: Does libMesh not calculate dphi and d2phi properly for lower-dimensional elements embedded at odd angles in higher-dimensional spaces? I always assumed that worked (despite having edited the code myself...), and I even vaguely recall telling someone el