I’ve found the new form compiler TSFC to be more robust with quadrature elements. It’s not included in FEniCS by default, but it can be installed manually using pip3
,
pip3 install git+https://github.com/blechta/tsfc.git@2018.1.0
pip3 install git+https://github.com/blechta/COFFEE.git@2018.1.0
pip3 install git+https://github.com/blechta/FInAT.git@2018.1.0
pip3 install singledispatch networkx pulp
then used in FEniCS Python scripts much like 'quadrature'
representation:
parameters["form_compiler"]["representation"] = 'tsfc'
One other trick is that, for especially complicated Form
s, I find that I sometimes need to manually raise Python’s default recursion limit, which can be done as follows:
import sys
sys.setrecursionlimit(10000)