Details about the examples can be found on the homepage of Prof. Sormann at http://itp.tugraz.at/LV/sormann/AKNumPhysik/
The first problem is solving the stationary heat equation with FEM.
Now how do we plot this? At first I tried using matplotlib which was not satisfactory because it cannot do color gradients across triangles in a mesh grid. Then I discovered the wonderful MayaVi package where you can do it as simple as
from enthought.mayavi import mlab
mlab.triangular_mesh(x, y, z, triangle_indices, temperature)
That's all! To add fancy stuff like a wireframe mesh, a colorbar and axes I did something like that before the final mlab.show()
mlab.triangular_mesh(x, y, z, triangle_indices, color=(0,0,0),
mlab.axes(extent=[0., 4., 0., 4., 0., 0.], nb_labels=5)
And here is the result:
If I have time, I'll keep you updated about the algorithm, the progress for the time-dependent heat conduction and more!