← Back to team overview

dolfin team mailing list archive

Re: [HG DOLFIN] Change test for cell-based data in VTKFile.

 

On Mon, Dec 22, 2008 at 01:28:26PM +0100, Martin Sandve Alnæs wrote:
> On Mon, Dec 22, 2008 at 1:14 PM, Garth N. Wells <gnw20@xxxxxxxxx> wrote:
> >
> >
> > Martin Sandve Alnæs wrote:
> >>
> >> On Mon, Dec 22, 2008 at 1:02 PM, Garth N. Wells <gnw20@xxxxxxxxx> wrote:
> >>>
> >>> Martin Sandve Alnæs wrote:
> >>>>
> >>>> This check isn't general either.
> >>>>
> >>> Can you give me an example of when it fails?
> >>
> >> It will fail for symmetric tensor elements.
> >>
> >
> > Which we don't yet support :).
> 
> Which SFC with UFL will support quite soon,
> and which our shared interface UFC allows :)
> 
> 
> >>>>  // Test for cell-based element type
> >>>> -  const uint tdim = mesh.topology().dim();
> >>>> -  bool only_cell_dofs = dofmap.needs_mesh_entities(tdim);
> >>>> -  for (uint i = 0; i < tdim; i++)
> >>>> -  {
> >>>> -    if (dofmap.needs_mesh_entities(i))
> >>>> -      only_cell_dofs = false;
> >>>> -  }
> >>>> -  if (only_cell_dofs)
> >>>> +  if( dofmap.local_dimension() == std::pow(mesh.topology().dim(), rank)
> >>>> )
> >>>>    data_type = "cell";
> >>>>
> >>>>
> >>>> This is the check I suggested to Anders earlier:
> >>>>
> >>>> const uint tdim = mesh.topology().dim();
> >>>> if(dofmap.num_entity_dofs(tdim) == dofmap.local_dimension())
> >>>>
> >>> What does num_entity_dofs(tdim) return?
> >>>
> >>> Garth
> >>
> >> The number of dofs associated with a cell (mesh entity of dimension tdim).
> >> dofmap.num_entity_dofs(0) would return the number of dofs associated
> >> with the vertices.
> >
> > OK. That would make the test easy.
> >
> >>
> >>>> Or hasn't FFC implemented UFC 1.1 yet?
> >
> > FFC hasn't implemented num_entity_dofs(tdim)
> >
> > Garth
> 
> It's more than ten months since UFC 1.1 was
> released, it would be nice if FFC implemented it.
> 
> Martin

It should be fairly simple to fix, but a lot of work is piled up for
FFC that I plan to look at after the release of DOLFIN 0.9.0.

Will the current fix work for a while?

-- 
Anders

Attachment: signature.asc
Description: Digital signature


Follow ups

References