← Back to team overview

ffc team mailing list archive

Re: facet normal implementation

 

On Sun, Aug 23, 2009 at 03:08:28PM -0500, Peter Brune wrote:
> The higher-order normal is the reference normal transformed by the Piola
> transform.  Interestingly, it's actually easier to get the local facet length
> out of this than it is to compute it "manually".  This is of course:
>
> n_iso = (1 / |J_surface|) * J * n_ref (I might be missing an inverse... it's
> getting late)
>
> Anyways, you can normalize the normal computed here to get the local jacobian
> determinant for the surface, because you can just use your higher-order cell
> jacobian to get your normal going in the right "direction, and of course
> it's... normalized, making life easier.  Two birds with one stone.  This would
> mean that  facet integrals would need the corresponding cell Jacobian, but
> that's the case with both affine and higher-order normals.  Right now I'm doing
> this by appending a higher-order normal function as well as higher-order
> coordinates to the form with a transformer.
>
> Thoughts?

Sounds good to me but I don't have time to fix it.

--
Anders

Attachment: signature.asc
Description: Digital signature


References