A nice example of Julia’s expressiveness, applied to a finite element algorithm manipulating sparse matrices. Especially interesting is the use of custom indexing functions.