Reference stiffness and mass matrix

I have prepared a reference stiffness A and a reference mass matrix M on a reference grid. Comparing with these reference matrices should be the first step when you debug your own assembler.

The mass matrix is assembled for homogeneous Neumann boundary conditions.

The grid itself is available as a Matlab file, or as an XML file that can be used with DOLFIN:

    grid.m
    grid.xml

The following Matlab file contains the reference matrices. It also contains a few hints on how to check that your matrices are correct.

    reference.m