Till,
I just took a look at the tests in deal.II:
http://www.dealii.org/svn/dealii/trunk/tests/codim_one/interpolation_3.cc
And there it seems that we don't have this problem. Are you sure you
are not trying to build a FESystem<0,1> object?
Luca.
On 21/set/09, at 11:53, Till Heinemann wrote:
The exception is due to the FESystem constructor's attempt to
dealii::FESystem<1,2>::initialize_unit_face_support_points(), as their
dimension is dim-1 referred to the FESystem's dimension, which is
already
dim-1=1, so we get a point with dim=0, which is not valid.
The names we use in deal.II are different from your notation.
dimension (dim) is the manifold dimension, and spacedimension
(spacedim) is the dimension of the Point<spacedim>, therefore in the
library the dim=1 case is with spacedim=2 (which is correct).
My guess is that you are using dim and dim-1 also when dim=1, and, of
course, this is not possible.
Luca.
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
