Anders Logg wrote:
> On Fri, Jun 19, 2009 at 12:38:28PM +0100, Garth N. Wells wrote:
>>
>> Anders Logg wrote:
>>> Sorry about that. I'll check on it later (if not fixed by
>>> then). Today's midsummer here in Sweden. :-)
>>>
>> The problem is that cell facets have not been initialised. Where should 
>> this be done?
> 
> It should be done in the constructor of DofMap.
>

This could be a FFC bug. For the Poisson DG example, looking at the 
generated contents of

     bool ufc::dof_map::needs_mesh_entities(unsigned int d) const

I see

    switch ( d )
       {
       case 0:
         return false;
         break;
       case 1:
         return false;
         break;
       case 2:
         return true;
         break;
       }

Shouldn't 'case 1' also return true (that is, cell edges are required)?

Garth

> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@fenics.org
> http://www.fenics.org/mailman/listinfo/dolfin-dev


_______________________________________________
DOLFIN-dev mailing list
DOLFIN-dev@fenics.org
http://www.fenics.org/mailman/listinfo/dolfin-dev

Reply via email to