Julia implementation of 2D Galerkin P1 finite element methods with parallel Karhunen Loeve computation and domain decomposition preconditioners. - View it on GitHub