Is it possible to specify custom integration schemes in FEniCS

Question

Dear All

I am planning to develop a molecular finite element code. In some elements
I require custom integration schemes to be able to evaluate
finite element matrices of the Coulomb potential.
Googling for this only brought up a related question with the answer that
it cannot be done. That would be a pity and would force me to
use dealII instead which , because it works with C++ requires a lot more programming.

regards

Moritz

Answer 1

Actually you can do it. In FEniCS 1.3.0 you can merely tweak ffc/quadrature_schemes.py. It is pretty straightforward.

Comments

Dear Jan

Thanks a lot

I require this functionality urgently

regards

Moritz