Towards an efficient meshfree solver