Many things...
There is finally an example of a stationary simulation with Dirichlet boundaries only.
Global matrices seem to work and the classes for solution of electromagnetic problems have been stated.
write_vtk function was upgraded to handle elements instead of polygons by using the datased unstructured grid.

Please look at the file called "small_post_proc_test.py" for an example on how to perform the simulation.

The vtk file produced by this simulation is called capa.vtk