On 21/set/09, at 11:53, Till Heinemann wrote:
Looking at the fe_system.cc I understand the 1D case is to be caught
by a
preprocessor condition (line 1833), checking deal_II_dimension, which
appears to coincide with the spacial dimension. But it doesn't catch
the
case of 1D objects in 2D space.
So I wanted to somehow change that condition in fe_system.cc.
In fe_system.cc (last changed, rev. 18933), I read:
#if deal_II_dimension == 1
template <>
void
FESystem<1>::
initialize_unit_face_support_points ()
{
// no faces no work
}
template <>
void
FESystem<1,2>::
initialize_unit_face_support_points ()
{
// no faces no work
}
#endif
Which should do the job correctly...
What version of deal are you using? If you want to work with BEM, I
strongly suggest you switch to the svn version.
Luca.
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
