Re: [Libmesh-devel] FENedelecOne

2012-09-05 Thread Roy Stogner
On Wed, 5 Sep 2012, John Peterson wrote: > On Wed, Sep 5, 2012 at 2:16 PM, Paul T. Bauman wrote: > >> 6. For higher-order tet elements, we'll need a tet14 for the face dofs. > > Interesting. 4 vertex + 6 midedge + 4 midface nodes? Right. We've needed a TET14 for a while; it would let us imple

Re: [Libmesh-devel] FENedelecOne

2012-09-05 Thread Paul T. Bauman
On Wed, Sep 5, 2012 at 3:45 PM, John Peterson wrote: > On Wed, Sep 5, 2012 at 2:16 PM, Paul T. Bauman wrote: > > > > 5. Just to be clear, the (excellent) report I cite in the code comes from > > the deal.ii website, but I came across it through a Google search. All > the > > code is my own, or c

Re: [Libmesh-devel] FENedelecOne

2012-09-05 Thread Roy Stogner
On Wed, 5 Sep 2012, Paul T. Bauman wrote: > The patch can be found here: > http://users.ices.utexas.edu/~pbauman/nedelec_one.patch > > Comments? OK for trunk? I'd give others some time to comment first, but all looks great to me. --- Roy ---

Re: [Libmesh-devel] FENedelecOne

2012-09-05 Thread John Peterson
On Wed, Sep 5, 2012 at 2:16 PM, Paul T. Bauman wrote: > > 5. Just to be clear, the (excellent) report I cite in the code comes from > the deal.ii website, but I came across it through a Google search. All the > code is my own, or copied from existing libMesh code, and nothing came from > deal.ii.

[Libmesh-devel] Handling edge/face orientation and element order

2012-09-05 Thread Paul T. Bauman
All, I wanted to solicit ideas and discussion for designs handling the edge/face orientation and element order. In particular, for scalar Lagrange elements of order >= 3 (and other scalar elements) or oriented vector elements (e.g. Nedelec), the orientation of the element matters and must be taken

[Libmesh-devel] FENedelecOne

2012-09-05 Thread Paul T. Bauman
All, I've created a patch against r5998 that implements first order Nedelec elements of the first type for triangles and quads, including a new vector_fe_ex3 which solves curl curl u + u = f for both the triangle and quads. The patch also includes updates to ExactSolution to compute HCurl and HDiv

Re: [Libmesh-devel] TensorTools

2012-09-05 Thread Roy Stogner
On Wed, 5 Sep 2012, Paul T. Bauman wrote: > On Wed, Sep 5, 2012 at 7:58 AM, Roy Stogner wrote: > > > I'm starting to waffle about whether or not we should include a > > libmesh_norm(a) { libmesh_deprecated(); return TensorTools::norm_sq(a); } > > shim. > > If this breaks, let me know and I'll

Re: [Libmesh-devel] TensorTools

2012-09-05 Thread Paul T. Bauman
On Wed, Sep 5, 2012 at 7:58 AM, Roy Stogner wrote: > > I'm starting to waffle about whether or not we should include a > libmesh_norm(a) { libmesh_deprecated(); return TensorTools::norm_sq(a); } > shim. > > That's probably just me. > > Go ahead and commit. > r5995 If this breaks, let me know and

Re: [Libmesh-devel] TensorTools

2012-09-05 Thread Roy Stogner
On Wed, 5 Sep 2012, Paul T. Bauman wrote: > I've made the suggested changes below in the updated patch here: > http://users.ices.utexas.edu/~pbauman/tensor_tools.patch > Tested in devel mode both --enable-everything and --enable-everything > --enable-complex > > OK for trunk? I'm starting to