allow radial basis constructor to default to comm world when appropriate
extension of elem->on_boundary() for a particular side: elem->on_boundary(s)
avoid a race condition when preparing multiple meshfree objects
create a level set for a particular boundary id
adding interface to compute gradients for MeshfreeInterpolation objects