dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

B1;2600;0cOn Fri, Jan 28, 2011 at 12:10:30PM +0100, Andre Massing wrote:
> Den 28.01.11 12.06, skrev Anders Logg:
> >On Fri, Jan 28, 2011 at 11:25:27AM +0100, Andre Massing wrote:
> >>Den 28.01.11 00.30, skrev Anders Logg:
> >>>On Thu, Jan 27, 2011 at 11:10:08PM +0100, Anders Logg wrote:
> >>>>I would like to add factory functions in UFC, in particular
> >>>>
> >>>>   finite_element* finite_element::create() const;
> >>>>   dof_map* dof_map::create() const;
> >>>>
> >>>>Are there any objections to this? It would be easy to add and easy to
> >>>>implement in FFC.
> >>>>
> >>>>I need it to implement refinement of FunctionSpaces in DOLFIN where
> >>>>fresh copies of ufc::finite_element and ufc::dof_map are needed given
> >>>>a ufc::finite_element or ufc::dof_map.
> >>>>
> >>>>An alternative is to hack create_sub_foo so that it returns a copy of
> >>>>the element/dofmap itself if the component argument 'i' is>= the
> >>>>number of sub elements/dofmaps.
> >>>
> >>>I suggest we add this to UFC 1.6 and make a release next week together
> >>>with DOLFIN 0.9.10, FFC 0.9.5 + other packages that need release. The
> >>>sooner the better: there will be time to make bug fix releases and
> >>>things will be tested before 1.0 in March.
> >>>
> >>>If there are any other things we want to add to UFC, now's a good
> >>>time.
> >>
> >>What about the tabulate_tensor function taking arguments for runtime
> >>quadrature rules, as we discussed a while ago? ATM I hacked UFC and
> >>my generated form files having an additional function
> >>
> >>virtual void tabulate_tensor(double* A,
> >>                              const double * const * w,
> >>                              const ufc::cell&  c,
> >>			     dolfin::uint num_quad_points,
> >>			     const double * const * coordinates,
> >>			     const double * weights) const
> >>
> >>Such an interface would add a flexible layer for those who wants to
> >>have some type of runtime quadrature.
> >>
> >>General opinons about that?
> >
> >The question is when this function will be implemented by form
> >compilers.
> >
> >We could add it and let FFC generate an empty function with throws an
> >exception.
>
> So I would like to delegate the following questions to the FFC
> experts: How difficult would it be to generate code for that?

I assume it would not be very difficult, but it would need to be
implemented by Kristian (as part of the quadrature code generator) and
he might have other things to do at the moment.

--
Anders