Point-based implementation of multilevel fast multipole algorithm for higher-order Galerkin's method