Reference stiffness and mass matrix

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

The stiffness matrix is assembled for homogeneous Neumann boundary conditions.

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

    mesh.m
    mesh.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