FYI: it is now fixed, so you can try the current git master.
On 04/18/2018 02:56 PM, Robert Cimrman wrote:
it is a bug that manifests when using a single quadrature point per element. It will be fixed soon, but in the meantime, try increasing the order to 2 in integrals. After that your script runs, but I am getting some nan values - you will need some more debugging of the code you inherited.