[dune-fem] Lumping Quadrature for Lagrange

[Claus-Justus Heine]

... my implementation of a lumping Dune::Fem quadrature can be found in
git on the branch acfem-playground (adaptive-conforming or continuous or
your favourite C-word FEM) in the proper sub-directory. Barely tested.

The stuff is modelled after the CachingQuadrature-class by replacing the
underlying traits-class and a supporting "low-level" lumping quadrature
class.

The intention is to be able to quickly exchange non-lumping quadrature
by lumping quadratures for testing purposes, without having to rewrite
more than 2 lines of code.

Cheers,

Claus

On 26.06.2013 21:32, Claus-Justus Heine wrote:
> Hi Andrea,
> 
> actually, cloning the code from Dune::Fem was harder than expected, but
> I have a first draft implementation which at least compiles; modelled
> after the CachingQuadrature. Still have to test that.
> 
> Thanks for your message; if it shows that the lumping quadrature
> actually work I can move it some branch in the git repo.
> 
> Best,
> 
> Claus
> 
> Am 26.06.13 18:26, schrieb Sacconi, Andrea:
>> Hi Claus,
>>
>> I think it hasn't been implemented yet. I have been using this
>> quadrature rule in my code for a while, implementing it myself. But
>> if you discover its existence, let me know :-)
>>
>> Cheers, Andrea
>> __________________________________________________________
>>
>> Andrea Sacconi PhD student, Applied Mathematics AMMP Section,
>> Department of Mathematics, Imperial College London, London SW7 2AZ,
>> UK a.sacconi11 at imperial.ac.uk
>>
>> ________________________________________ From:
>> dune-fem-bounces+a.sacconi11=imperial.ac.uk at dune-project.org
>> [dune-fem-bounces+a.sacconi11=imperial.ac.uk at dune-project.org] on
>> behalf of Claus-Justus Heine
>> [Claus-Justus.Heine at IANS.Uni-Stuttgart.DE] Sent: 26 June 2013 17:23
>> To: dune-fem at dune-project.org Subject: [dune-fem] Lumping Quadrature
>> for Lagrange
>>
>> Hi,
>>
>> another question: is there already a "lumping" quadrature for
>> Lagrange elements on simplex grids, i.e. a quad rule with weight
>> 1/dim and integration points at the vertices of a simplex? Its
>> extremely simple to implement, copying the code from
>> dune/fem/quadrature/ but if it is already there ... and yes, there
>> are much more efficient ways to implement mass-lumping, but ...
>>
>> Best,
>>
>> Claus
>>
>> -- Dr. Claus-Justus Heine Institut für Angewandte Analysis und
>> Simulation * Numerische Mathematik für Höchstleistungsrechner
>> Universität Stuttgart
>>
>> Fon: +49 (0) 711 - 685 65558 Fax: +49 (0) 711 - 685 65507
>>
>> _______________________________________________ dune-fem mailing
>> list dune-fem at dune-project.org
>> http://lists.dune-project.org/mailman/listinfo/dune-fem
>>
>> _______________________________________________ dune-fem mailing
>> list dune-fem at dune-project.org
>> http://lists.dune-project.org/mailman/listinfo/dune-fem
>>
> 
> 


-- 
Dr. Claus-Justus Heine
Institut für Angewandte Analysis und Simulation
* Numerische Mathematik für Höchstleistungsrechner
Universität Stuttgart

Fon: +49 (0) 711 - 685 65558
Fax: +49 (0) 711 - 685 65507