On Thu, May 6, 2010 at 7:22 PM, <[email protected]> wrote: > Hi James, > > actually you'd have to exchange the cases Wolfgang gave below, see the > documentation of GeometryInfo<dim>. > > Furthermore, you might also consider the following alternative approach: > a) Compute the volume/area of the cell by numerical quadrature (using > FEValues). > b) Compute the area/length of the faces of the cell by numerical quadrature > (using FEFaceValues). > c) Obtain an estimate of the length in one direction by dividing the cell > volume by the mean value of the area of the two faces that are orthogonal to > that direction (faces 0 and 1 for direction x(or 0)).
Thanks a lot to both of you! I'll have a go at both methods today and see what's easiest -- I'll post again with the result. Cheers, -- Med venlig hilsen, James Avery <[email protected]> _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
