PDEAssembles takes a set of partial differential equations and boundary conditions and creates a global assembly M, C, K, f in SparseMatrixCSC format ready for linear system solvers. - View it on GitHub